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ПРЕДИСЛОВИЕ РЕДАКТОРА 


Дорогие читатели журнала, уважаемые авторы и коллеги! 

В конце 2004 года произошло важнейшее событие в жизни нашего журнала «Хи- 
мическая физика и мезоскопия». Решением Президиума Уральского отделения Россий- 
ской Академии наук утвержден новый состав редколлегии и журнал признан регио- 
нальным академическим журналом. 

Состав редколлегии существенно увеличен, что продиктовано потребностями 
быстрого и качественного рецензирования поступающих статей для ускорения выпус- 
ка номеров журнала. 

В журнале расширяется круг помещаемых материалов. Кроме заявленных разде- 
лов по тематике журнала появится уже в этом номере раздел «Персоналии», в кото- 
ром будут сведения о выдающихся ученых и специалистах, работающих в области хи- 
мической физики и мезоскопии. Планируется раздел, отражающий достижения моло- 
дых ученых и студентов, а также информация о достижениях научных коллективов, 
касающихся тематики журнала, и рекламные объявления по новым наноматериалам, 
процессам нанотехнологии, наноразмерным методам исследования, соответствующей 
приборной и технологической технике, наряду со сведениями по новым программным 
продуктам. 

Обращаем внимание спонсоров и рекламодателей на то, что журнал становится 
подписным изданием и имеет выход в Интернет через сайты Удмуртского научного 
центра Уральского отделения Российской Академии наук, Института прикладной ме- 
ханики УрО Р АН, Физико-технического института УрО РАН, Удмуртского государ- 
ственного университета, Ижевского государственного технического университета, 
Министерства образования и науки Р оссийской Федерации. Ваша информация через 
Интернет станет доступной для научной общественности не только Уральского от- 
деления РАН и сотрудников Вузовской науки Урала и Поволжья, но и ученых и специа- 
листов всей России. 

В данном номере журнала опубликованы работы ученых Москвы, Екатеринбурга 
и Ижевска, помещены Персоналии двух выдающихся ученых, юбилеи которых отмече- 
ны в начале 2005 года. Речь пойдет о нашем главном редакторе, академике РАН Алек- 
сее Матвеевиче ЛИПА НОВ Е (70 лет, 3 марта 2005 года) и о члене редколлегии, про- 
фессоре, докторе химических наук Геннадии Ефремовиче ЗАИКОВЕ (70 лет, 7 января 
2005 года). 

Редколлегия журнала сердечно поздравляет юбиляров и желает им крепкого здо- 
ровья и творческого долголетия. 

Заместитель главного редактора 

журнала «Химическая физика и мезоскопия» В. И. КОДОНОВ 
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МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ФИЗИКО-ХИМИЧЕСКИХ 
ПРОЦЕССОВ 

УДК 536.7+541.12 

ПРИМЕНЕНИЕ ФУНКЦИЙ ЛАГРАНЖА И ПРОСТРАНСТВЕННО - 
ЭНЕРГЕТИЧЕСКИХ ПАРАМЕТРОВ ДЛЯ ОПРЕДЕЛЕНИЯ 
ВЕРОЯТНОСТИ ПРОЦЕССОВ 

А.М. ЛИПАНОВ, В.И. КОДОЛОВ, Г.А. КОРАБЛЕВ 

Научно-образовательный центр химической физики и мезоскопии УдНЦ 
УрО РАН, Ижевск, Россия 


АННОТАЦИЯ. Предложен метод оценки вероятности протекания процессов с ис- 
пользованием функций Лагранжа. Обосновано введение пространственно - энерге- 
тического параметра, полученного на основе анализа поверхностной и объемной 
энергий систем и вывода приведенной энергии. Приведены примеры применения 
предложенного метода для расчета энергий активации химических реакций. 

ВВЕДЕНИЕ 

Для наносистем и процессов с их участием вычислительный аппарат на основе кванто- 
вой химии или на основе термодинамических представлений может дать неадекватную 
картину реальных процессов и его использование вызывает сомнения в достоверности 
получаемых результатов. К сожалению обоснованного перехода от квантово- 
химических расчетов, выполненных для пикосистем, к термодинамическим, примени- 
мым для макросистем, в настоящее время нет. Единый вычислительный аппарат, кото- 
рый можно применить к системам любого размера, только создается как единая сеть 
превращений пикосистем через наносистемы в макросистемы с учетом путей перехода, 
обусловленных направленностью взаимодействий компонентов систем. 

Здесь может использоваться для характеристики энергии систем принцип адди- 
тивности, а для оценки энергий взаимодействий и активности систем понятия «свобод- 
ной энергии» и « энергии активации» . 

Нобелевский лауреат И. Пригожин предложил для оценки условий максимального 
выхода реакции использовать неопределенные множители Лагранжа [17]. Эти множи- 
тели или функции использованы при применении решетчатой модели раствора, образо- 
ванного молекулами двух или более сортов [18] полимеров. 

Таким образом, возникает возможность, используя функции Лагранжа, провести 
оценку процессов образования наночастиц и наносистем из пикосистем, микрочастиц и 
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микросистем из соответствующих наноразмерных структур и, наконец, макрочастиц и 
макросистем из организованных микросистем. В отличие от квантово-химических рас- 
четов, основанных на анализе всех электронных состояний, т.е. полной энергии систе- 
мы, когда для упрощения расчетов нужно подбирать потенциалы полей для атомных 
остовов, использование Лагранжиана сокращает количество взаимодействующих элек- 
тронных облаков до того минимума, при котором происходит наиболее выгодные с 
энергетической точки зрения процессы. 

ОБОСНОВАНИЕ МЕТОДОЛОГИИ 


Вероятность протекания процессов в химических системах обусловлена изменением 
электронных плотностей в реагирующих связях и перераспределением электронной 
плотности в переходном состоянии (активированном комплексе). Оценку таких изме- 
нений можно провести с помощью квантово-химических вычислительных эксперимен- 
тов, которые в основном базируются на волновом уравнении: 


Ну/ = Еу/ (1 ) 

При этом учитывается полная энергия, включающая кинетическую и потенциаль- 
ную энергии системы. Используемый в уравнении Гамильтониан записывают как 


где Т - кинетическая энергия, 11 - потенциальная энергия. При этом обычно кинетиче- 
ская энергия включает вклады энергий поступательного, вращательного и колебатель- 
ного движений, а потенциальная энергия соответствует энергиям притяжения и оттал- 
кивания. Сложные системы можно мысленно разбить на составляющие фрагменты, 
>читывая, что каждый фрагмент остается в системе, давая определенные вклады в ее 
кинетическую и потенциальную энергии. Чем активнее фрагмент, тем больше разница 
между его вкладами в кинетическую и потенциальную энергии. 

Отметим, что разность между кинетической и потенциальной энергиями не что 
иное, как Лагранжиан - 


Ь = Т-П 

Вместе с тем функции Гамильтона и Лагранжа взаимосвязаны: 


(3) 


Н = ^РЧ~ 1 

где 


( 4 ) 

( 5 ) 
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А.М. ЛИПАНОВ, В.И. КОДОЛОВ, Г.А. КОРАБЛЕВ 


По смыслу анализируемых функций представляется целесообразным рассматривать 
реакционность отдельных фрагментов системы, атомов, молекул, ионов, радикалов и 
иных химических частиц с применением функций Лагранжа, поскольку активность 
частиц реализуется в основном через соответствующие им движения и мерой активно- 
сти может служить разность кинетической и потенциальной энергий. Если рассматри- 
вать химические связи в химических частицах, в которых ограничены поступательные 
и, в меньшей степени вращательные движения, то, по-видимому, все движения можно 
свести к колебаниям. Параметрами колебательных движений атомов в связях, во мно- 
гом определяемых распределением электронной плотности в связи, являются частоты и 
амплитуды колебаний (частоты и длины волн связей). Увеличение амплитуды колеба- 
ний приводит в конце концов к разрушению связи. В свою очередь разрушение связей 
обусловлено в большинстве случаев возбуждением электронных оболочек, вплоть до 
ионизации отдельных атомов, находящихся в связи до её разрушения. Разность между 
энергией колебаний и энергией притяжения отражает по смыслу функцию Лагранжа 
(Ь). Если последовательно рассматривать устойчивость или реакционную способность, 
двигаясь от атомов к более сложным многоатомным образованиям, то функции Ла- 
гранжа и его уравнения возможно применить к любым пикоразмерным частицам. 
Можно отождествить энергию колебаний с поверхностной энергией, а энергию притя- 
жения - с объемной энергией. Рассматривая их как системы, состоящие из взаимодей- 
ствующих между собой внутренней и поверхностной частей или «ядра» и «оболочки», 
обладающих своими энергиями и массами, можно записать общие приведённые массу 
и энергию в виде 


а 


П) 


Е \ Е 1 

Е \ +е 2 


(6) 


т х т 2 
т х +т 2 ’ 


( 7 ) 


где е І и т, - энергия и масса «ядра», е 2 и пц- энергия и масса «оболочки». Здесь можно 
определить энергию «оболочки» как поверхностную энергию, а энергию «ядра» как 
объёмную энергию. Для атома объёмная энергия е, может выражать его потенциаль- 
ную энергию, определяемую энергией притяжения электронной оболочки к атомному 


ядру согласно уравнению 



( 8 ) 


где ^ - — , а 2 и п 
п 


эффективный заряд ядра атома и эффективное главное квантовое 


число, Гі - орбитальный радиус і-той орбитали химической частицы (для атома - орби- 
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тальный радиус). Энергия е 2 определяется подвижностью электронного облака внеш- 
них орбиталей, что может быть выражено через орбитальную энергию, которую можно 
представить или через энергию ионизации атома (Еі), или через энергию связи электро- 
нов на данной орбитали (\ѵ,): 

е' 2 = Еі; в 2 " = ѵѵіПі (9,9а) 

где Еі — энергия ионизации атома, ѵѵ, — орбитальная энергия, н, — число электронов на і- 
той орбитали. Подставив эти значения в уравнение приведённой энергии, получим па- 
раметр, который в работе [3] назван эффективным пространственно-энергетическим 
параметром Р э . 

Р э ' = (ЕіЯі 2 /гі)/(Еі + Чі 2 /гі); Р э " = (\ѵ;ПіЯі 2 /гі)/(яі 2 /гі + щщ) (10,10а) 

Это выражение можно преобразовать в уравнение 


р_ Я 2 Е,г, 

* п 


+ Е,Гі 


( 9 ) 


При этом Р 0 = Р э хгі (Па) 

В данном параметре находят отражение два противоположных, точнее противодейст- 
вующих, фактора. С одной стороны, в атоме отрицательно заряженные электронные 
оболочки «прижимаются» к ядру (и), с другой стороны, более подвижные валентные 
электроны на поверхности атома стремятся «скомпенсировать» накопленную энергию с 
окружающей средой (е 2 ). Поскольку направления их действий противоположны, то чем 

меньше Е„ тем активнее частица и наоборот, чем больше, , тем она более стабиль- 
на. По физическому смыслу эти понятия близки к кинетической и потенциальной энер- 
гии в функции Лагранжа. Вместе с тем полученный параметр характеризует состояние 
атома или другой химической частицы или компонента в окружающей среде при опре- 
делённых внешних условиях. 

Таким образом, следует, что пространственно-энергетический параметр Р 0 являет- 
ся своеобразным «паспортом» химической частицы и близок по смыслу к химическому 
потенциалу д, 


Мі = 


г дС] (дГ^ 

у Вп ‘кр,»« У 8п ^ту,п,„ 


{ дн"' 






У дП ЧВ,Ѵ,п 


( 10 ) 
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О - свободная энергия Гиббса, Р - свободная энергия Гельмгольца, Н - энтальпия, И - 
внутренняя энергия, п, - число молекул і-го компонента, 8 - энтропия, р - давление, V - 
объем. 

Такое представление открывает возможность применения Р-параметра для оценки 
вероятных взаимодействий по аналогии с известными способами использования хими- 
ческого потенциала. Поэтому вводится безразмерный параметр а, равный относитель- 
ной разности Р-параметров взаимодействующих атомов-компонентов 

си = ДРіАТРі/х), (и) 

где х — количество участвующих во взаимодействии атомов. 

Для парных взаимодействий - 


а = 


АР 

(Щ 

2 


(12) 


В работе [3] параметр а используется в качестве коэффициента, характеризующего 
изоморфизм компонентов, по аналогии с использованием химического сродства - 

А== -^ Ѵ Л (13) 


где д, химические потенциалы компонентов, ѵ, - количество частиц, участвующих во 
взаимодействии компонентов. При изоморфизме или изоструктурных превращениях 
чем меньше А и а, тем более вероятны превращения. Если эта аналогия верна, то с 
уменьшением а должна возрастать скорость взаимодействия. Вероятность протекаю- 
щих процессов обычно оценивают по константе скорости процесса, которая определя- 
ется по уравнению Аррениуса - 


к = к 0 ехр 



(14) 


где ко предэкспоненциальный множитель, зависящий от изменения энтропии процес- 
са, Е а - энергия активации процесса. Естественно, в большей степени влияет на ско- 
рость процесса именно энергия активации. Поэтому необходимо оценить, как про- 
странственно-энергетические параметры влияют на энергию активации и предэкспо- 
ненциальный множитель. При термодинамическом подходе к анализу уравнения (16) 
можно записать - 


/ КТ ( 

к = Х^т~ ех Р 
п 


АС 

КТ 




(15) 
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ИЛИ 



(16) 


где трансмиссионный множитель, кь — константа Больцмана, Т - температура, Ь — 
постоянная Планка, ДО* - изобарно-изотермный потенциал, Д8 # - изменение энтропии в 
процессе, ДН - изменение энтальпии при переходе в активированное состояние. 

В свою очередь 


Е а = ДН* + КТ 


(17) 


Если вышеизложенный подход верен, то при введении относительных величин, т.е 
проведении операции обезразмеривания, поскольку а является безразмерным 
параметром, можно представить - 


с1Е а /Е а = па 

После интегрирования формулу для Е а можно преобразовать к виду: 
Е а =аехр(псс) = аещ){^Ру , 


(18) 


(19) 


где а - переводной коэффициент в систему СИ (4,184 кДж/моль), п - число 
взаимодействующих компонентов, / - коэффициент, приводящий значения Р к а. 

ПРИМЕНЕНИЕ ПРОСТРАНСТВЕННО-ЭНЕРГЕТИЧЕСКИХ ПАРАМЕТРОВ 
ДЛЯ ОЦЕНКИ ЭНЕРГИЙ АКТИВАЦИИ РЕАКЦИЙ И ВЕРОЯТНОСТИ ИХ 
ПРОТЕКАНИЯ 

Для подтверждения приведённых аргументов в пользу применения функций Лагранжа 
и выведенных пространственно-энергетических параметров для оценки энергий 
активации взаимодействия атомов и отдельных фрагментов («блоков») в газовой фазе. 
При этом были внесены уточнения и дополнения в предполагаемую схему расчётов Р- 
параметров. Например, для ионов в формуле для Р и учитывается радиус иона - 



Для сложных систем расчёт Р-параметров проводился по формулам: 



( 21 ) 
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_І_ 

Я 


1 


1 


N Р 

1 ѵ у- п 


N Р" 

ІѴ 2 1 п 


( 22 ) 


Здесь N 1 и Ы 2 - числа однородных атомов в молекуле (например для Н 2 N=2). 

Введение множителей N 1 и N 2 означает, что в данной модели используется и 
принцип аддитивности Р-параметров. Например, система СН 2 СОСНз рассматривалась, 
как состоящая из трех таких блоков СН 2 , СО и СН3. Р-параметр каждого из них 
находился в соответствии с уравнениями (10, 11), затем алгебраическое сложение Р і? Р и 
— параметров этих блоков дает суммарный Рс — параметр всей структуры. При этом 
радиус атома водорода значительно зависит от характера химической связи. Для атома 
водорода в основном состоянии в качестве размерной характеристики использовался 
боровский радиус (0,529 А) и рассчитывался Р|° -параметр (для электронейтрального 
атома). Эта же величина Р-параметра водорода использовалась и для структур с 
ковалентной связью, то есть предполагалось, что Р^ = Р и = 9,0624 эВ. Для реакций, где 
возможно взаимодействие с ионом Н’ использовался ионный радиус водорода (1,36 А), 
что дает в расчете Р и " = 3,525 эВ для иона водорода. 

Для расчета энергий активаций свободных атомов с компонентами 
использовались рассчитанные значения Р-параметров, приведенные в табл. 1 , 2. 

Следует отметить, что при а -> 0, Е а -» 4,184 (^і). Тогда как значения Е а 

моль 

(в эксперименте) находятся для этих случаев (табл. 3) обычно в пределах от 0 до 4 
КДж 


). Иначе говоря уравнение (19) менее точно работает при значениях а — >0. В 


( 

моль 

табл. 4 приведены расчеты энергии активации газофазных реакций с участием 
галогеноводородных молекул, ионный характер которых учитывался через 

Р 


коэффициент у = 


Р 


НС1 


где ц - дипольный момент данного соединения; цнсі - 


дипольный момент молекул НС1. В данном случае уравнение (8) принимает вид: 


Е а = ауе па 


(23) 


где п = 1,2, 3... 

В табл. 5 приведены выборочные расчеты Е а для молекулярных реакций по 
уравнениям (19 и 23) и по уравнению: 


а 


І.И 


Р 1 - Р 11 
г |,и г с 

(Р‘и+Рс)/2 


(24) 


Анализ табл. 5 приводит к выводу о том, что уравнение (19) и в этих случаях 
работает достаточно удовлетворительно. Точность расчетов обычно находится в 
пределах 10-15%, что не превышает относительной погрешности большинства 
экс периментов. 


12 


ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. Том 7, № 1 


Таблица 1. Р - параметры атомов 


ПРИМЕНЕНИЕ ФУНКЦИЙ ЛАГРАНЖА И ПРОСТРАНСТВЕННО - ЭНЕРГЕТИЧЕСКИХ 
ПАРАМЕТРОВ ДЛЯ ОПРЕДЕЛЕНИЯ ВЕРОЯТНОСТИ ПРОЦЕССОВ 


Приме- 

чание 

для Н' 




о 

§ 

для О 1 



для К 3 ' 


для N 

для 0' 

для (У' 






о 

си 

н 

1 

с 

03 

о 

Н 

3,525 




19,900 

258,7 



22,746 


550,9 

3,8419 

12,724 

4,3774 

О 

оч 

ОО 

4,5199 

4,3486 

4,7898 

Ги (А) 

1,36 




2,60 

0,20 



ОО 


ЦП 

о 4 

1,36 

1,36 

1,345 

г—Ч 

ОО 

1,96 

2,20 

0,98 

о 

Он 

н 

1 

с 

Ін" / Ч 

03 

о 

н 

9,0624 


26,750 


86,810 

12,822 

36,111 

68,984 

108,69 

158,62 

12,621 

41,797 

16,389 

І9І‘П 

10,410 

9,1638 

2,7402 

■< 

оэ 

т 

о 

Он 

н 

4,794 


15,943 


51,739 

г-~ 

сч 

чо 

17,586 

33,664 

56,63 

82,642 

5,225 

17,304 

5,882 

8,125 

8,859 

9,567 

4,694 

Ро(эВА) 

4,794 

5,641 

10,302 

16,515 

19,281 

6,257 

11,329 

ОО 

о 

<о 

40 

22,966 

26,012 

5,225 

12,079 

04 

ОО 

оо^ 

8,125 

8,859 

9,567 

4,694 

Ч 2 (эВА) 

14,398 

35,395 

35,395 

37,243 

37,243 

52,912 

52,912 

52,912 

53,283 

53,283 

71,380 

71,380 

94,641 

59,842 

73,346 

77,65 

10,058 

2 

Іч 

0,529 

0,596 

0,596 

0,620 

0,620 

0,488 

0,487 

0,487 

0,521 

0,521 

0,414 

0,414 

0,360 

0,728 

0,851 

1,044 

1,713 

Е (эВ) 

13,599 

11,260 

24,383 

47,860 

64,480 

14,54 

29,60 

47,426 

77,472 

97,89 

13,618 

35,118 

17,423 

12,268 

11,84 

10,451 

ОО 

СП 

ЦП 

Валентные 

орбитали 

сл 

"а 

04 

Ъ 

СЧ 

СО 

04 

С/Э 

04 

04 

04 

“си 

сч 

"с/э 

04 

04 

“си 

гч 

"Ь, 

СЧ 

~а 

ГЧ 

"а, 

СП 


“а 

ІЛ 

СО 

СП 

Атом 

I 

о 

2 

О 

Си 

г— < 

и 

Вт 

НН 

аз 

г 


ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. Том 7, № 1 


13 


Таблица 2. Структурные Р- параметры 
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Таблица 4. Энергия активации газофазных реакций галогеноводородных систем ( 
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ПРИМЕНЕНИЕ ФУНКЦИИ ЛАГРАНЖА И ПРОСТРАНСТВЕННО - ЭНЕРГЕТИЧЕСКИХ 
ПАРАМЕТРОВ ДЛЯ ОПРЕДЕЛЕНИЯ ВЕРОЯТНОСТИ ПРОЦЕССОВ 
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При этом п = 1, 2, 3, 4, 5, 6. 

Отметим, что основному "разрешенному" состоянию порогового значения 
энергии (энергии активации) соответствует п = 3 — для систем, содержащих радикалы и 
молекулы (как и для реакций с участием свободных атомов). 

Таким образом, данный подход дает достаточно рациональный, практически 
удобный и научно-обоснованный метод расчета энергии активации химических 
реакций. 
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МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ТЕЧЕНИЯ ЖИДКОСТИ С 
ПЕРЕМЕННОЙ СТРУКТУРОЙ 

А.М. ЛИПАНОВ, А. А. ВАХРУШЕВ, А.В. ВАХРУШЕВ 


Институт прикладной механики УрО РАН, Ижевск, Россия 


АННОТАЦИЯ. Сформулирована задача расчета течения жидкости, в которой при 
изменении макропараметров течения изменяются её микроструктурные свойства и 
реологическое поведение. Выполнены численные исследования течения жидкостей с 
температурным пределом изменения вязкости. 


ВВЕДЕНИЕ 

В связи необходимостью расчета микро - и нанотечений жидкостей, весьма 
актуальным стало исследование потоков, в которых происходят структурные измене- 
ния [1-3]. В таких потоках при варьировании макропараметров течения (давление, 
температура и др.) изменяются микроструктурные параметры жидкости, а, 
следовательно, и реологическое поведение. При этом происходит дискретное 
(скачкообразное) увеличение или уменьшение коэффициента вязкости. Понимание 
процессов, сопровождающих течение такого типа жидкостей, играет важную роль при 
выявление фундаментальных закономерностей влияния изменения структурных 
свойств на параметры течения, а также при проектировании и оптимизации различных 
промышленных и технологических процессов. 

В настоящей работе, в развитие ранее выполненных исследований [4.5], 
приведены результаты расчетов течения жидкостей с температурным пределом изме- 
нения её структуры и вязкости. 

ПОСТАНОВКА ЗАДАЧИ 

Постановка задачи расчета течения [6,7] со структурным изменением вязких свойств 
включает в себя, наряду с гидродинамической моделью, модель реологического 
поведения жидкости с дискретным изменением вязкости, а также температурную 
зависимость структуры от характеристик потока. При этом граница области изменения 
вязкости не имеет заранее заданного положения, а передвигается в зависимости от эво- 
люции температуры и других параметров потока. Вязкость на границе раздела областей 
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с различной структурой жидкости меняется “ступенчато”. Таким образом, имеем за- 
дачу течения “двухфазовой” жидкости с переменной границей раздела фаз. 

Система уравнений, описывающая данную задачу, включает: обобщенную 
систему уравнений гидродинамики (1), (2) и уравнение энергии (3) в безразмерных 
переменных: 


да ду/ да ди/ да 

— + = А(п- а) + сь , 

ді ду дх дх ду 


д'у/ д'у/ 

— т + — т = -&>, 
дх 2 ду~ 


дт д у/ дт д у/ дт і т г 

— - — | = АТ А Ф т, 

ді ду дх дх ду Ре Ке 


Ф =-2 


д 2 г/ д 2 у/ 

дх 2 ду 2 дхду дхду ду 2 дх 2 


- 2 - 


д'т] д 2 у/ + д 2 ?і д 2 у/ 


(1) 

( 2 ) 

( 3 ) 

(4) 


Ф т =4 


д 2 


V 


дхду 


+ 


( я 2 
д у/ 


ду 2 дх 2 


О I// 


(5) 



Рис-1. Продольное сечение прямого канала: а-Ь - участок нагрева или охлаждения 
стенки канала 
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где со - завихренность потока; у/ - функция тока; Т - температура потока; 
7 = т/(х,у,...) - вязкость жидкости; I, х, у - пространственно-временные координаты; 

л 5 2 52 

А = — у 3 7 - оператор Лапласа; Ке - число Рейнольдса; Ре - число Пекле; Е - 

дх ду 

число Эккерта; Ф щ - диссипативная функция в уравнении переноса “вихря”; Ф г - 
диссипативная функция в уравнении энергии. 

Система уравнений ( 1 )-( 5 ) дополняется дискретной зависимостью вязкости от 
температуры 


V = г] 0 ,Т <Т 0 ; 
П = г]\,Т >Т 0 , 


(6) 


где 7 0 , 7, - дискретные значения вязкости; Т 0 - температура структурного перехода. 

Сформулированная задача, в общем случае, не имеет аналитического решения, 
поэтому были использованы конечно-разностные численные методы, реализованные в 
пакете программ расчета для однопроцессорных и мультипроцессорных вычислитель- 
ных комплексов [8]. 

РЕЗУЛЬТАТЫ РАСЧЕТОВ 

В качестве примера, приведем расчет течения, в котором, в основном потоке, коэффи- 
циент вязкости равен т ] шіп , а при охлаждении среды до определенной температуры про- 
исходит его дискретное увеличение от г/ шіп до г / тах , то есть изменяется реологическое 
поведение жидкости. Характерное число Рейнольдса, соответствующее невозмущен- 
ному потоку в канале, — Ке = 10 при вязкости потока г] тт . 

Т аким образом, в ходе вычислений рассматривалась “двухфазовая” среда; исполь- 
зовался температурный предел изменения вязкости, равный 30 % от среднего значения 
температуры Т в невозмущенном потоке. Результаты расчетов приведены на рис. 2 
(линии тока), рис. 3 (поле температур), рис . 4 (поле вектора скорости потока). 

Процессы, происходящие в потоке при возникновении локальной области струк- 
турного изменения вязких свойств среды, аналогичны обнаруженным в [4, 5] при моде- 
лировании неподвижного скачка вязкости. Происходит замедление потока в зоне охла- 
ждения и его ускорение (в силу условия постоянства расхода) вблизи противоположной 
стенки канала. 

Видно (рис. 2 -Ь), что при коэффициенте увеличения вязкости т] тах Іт] тііі = 12.5 
происходит образование циркуляционного течения. Наличие обратного основному по- 
току движения жидкости приводит к отклонению температурного фронта (рис. 3-Ъ) 
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Рис.2. Линии тока при течении жидкости со структурным изменением вязкости 
при коэффициенте увеличения вязкости т) тах /т ] шш , равном: а - 10; 
Ь - 12,5; с - 15; сі - 17,5; е - 20 

вблизи охлажденной стенки. Вмести с этим, изменяется положение границы области 
повышенной вязкости. 

Описанные явления приводят к нестационарности процесса. Поле скоростей по- 
стоянно перестраивается под влиянием смещения границы структурного увеличения 
вязкости, что, в свою очередь, приводит к изменению распределения температуры в по- 
токе. Возникают колебательные движения, и течение не выходит на стационарный ре- 
жим. С ростом коэффициента увеличения вязкости г] тах Іт] шіп амплитуда колебания зна- 


26 


ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. Том 7, № 1 


МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ТЕЧЕНИЕ! ЖИДКОСТИ С ПЕРЕМЕННОЙ СТРУКТУРОЙ 



8 


Рис.З. Распределение поля температур при течении жидкости со структурным из- 
менением вязкости 


чений в поле скоростей увеличивается, возникают дополнительные циркуляционные 
зоны, хорошо заметные на рис. 4 при изменении направления векторов скорости в раз- 
личных точках расчетной области. 

Число дополнительных зон ооратного течения с ростом вязкости увеличивается 
от одной (рис. 4-6) до двух (при росте вязкости в 20 раз, соответствующем рис. 4-е): 
образуются основная циркуляционная зона в области охлаждения и две вторичные 
вблизи противоположной стенки. 
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а 



Ь 




е 


Рис.4. Направление вектора скорости ѵ при росте вязкости в 10,12.5,15,17.5,20 
раз (течение жидкости со структурным изменением вязкости) 
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Рис.5. Положения границы повышенной вязкости в различные моменты времени 

При дальнейшем увеличении максимального значения вязкости т] тах в области 
охлаждения поток теряет устойчивость, положение границы «фазового» изменения 
вязких свойств постоянно изменяется во времени (рис.5). 

Управляющим параметром рассмотренной реологической модели структурного 
изменения вязкости может являться значение температурного предела. Так, при 
уменьшении степени охлаждения жидкости, необходимой для изменения вязкости, по- 
явление циркуляционных зон происходит при меньших значения коэффициента увели- 
чения вязкости г] тах / Т) тіп и наоборот. 

Таким образом, выполненные численные показывают, что изменение структуры 
жидкости существенно изменяет параметры течения, вызывая появление циркуляций в 
потоке и различные нестационарные явления. 
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АННОТАЦИЯ. Приведена методика расчета энергии активации самодиффузии и 
объемной диффузии атомов в твердых телах, основанная на эффективной энергии 
межатомного взаимодействия. Результаты проведенных многочисленных расчетов 
удовлетворительно согласуются с экспериментальными данными по диффузии эле- 
ментов в твердых телах. 


ВВЕДЕНИЕ 

В таких технологических методах производства как сварка, наплавка, термическая об- 
раоотка и ряда других основной вклад в образование соединений между материалами, 
формирование их структуры и получаемые при этом свойства вносит процесс диффу- 
зии. Его основной характеристикой является энергия активации направленного пере- 
мещения атомов химических элементов в кристаллической решетке твердого тела. 

Теоретическим расчетам и экспериментальным методам определения энергии ак- 
тивации посвящено достаточно много исследований [1-3]. Однако сопоставление лите- 
ратурных данных выявило наличие весьма существенных противоречий [4]. Например, 
для хорошо изученного процесса самодиффузии кремния расхождение теоретических и 
экспериментальных значений энергии активации составляет 1,56 эВ, а значения для 
диффузии железа в сплаве Ре+1%Мп, полученные на основе метода радиоактивных 
изотопов, отличаются от значений полученных с помощью металлографического мето- 
да [5] на 2,37 эВ. Кроме того, анализ приведенных в работе [6] результатов показал, что 
диффузионная подвижность атомов компонента В в окружающей его среде атомов 
компонента А существенно меняется при наличии в ней атомов третьего компонента С, 
как в сторону увеличения, так и уменьшения в зависимости от свойств последнего. 

Согласно современным представлениям о строении твердых тел диффузия атомов 
происходит в основном по вакансионному механизму, так как для этого требуются 
наименьшие затраты энергии по сравнению с другими возможными механизмами [1-3, 
7,8]. Миграция атома в соседний с ним узел или междоузлие возможна только при со- 
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общении ему дополнительной энергии, которая равная по величине энергии того по- 
тенциального барьера, который он должен преодолеть. Установлено, что данная вели- 
чина общей энергии атома равна сумме энтальпий образования и перемещения вакан- 
сий [2]. Поскольку атомы твердых тел состоят из ионов и относительно свободных 
электронов, движущихся в определенных пространственных зонах, можно предполо- 
жить взаимосвязь между энергией активации диффузии и электронным взаимодействи- 
ем атомов в кристаллической решетке тела. Исходя из этого, энергию активации (Е а ) 
рассмотрим как энергию межатомного (парного) взаимодействия атомов А окружаю- 
щей среды и диффундирующего в ней атома компонента В. Если энергия поступатель- 
ного движения атома В больше энергии межатомного взаимодействия с атомами А, то 
он диффундирует за пределы окружающих его атомов А. Если же атом В будет иметь 
энергию движения меньше, чем энергия его взаимодействия с атомами А, то он не 
сможет мигрировать за пределы этих атомов. В этом случае не будет происходить об- 
менного процесса между атомами А и В, то есть диффузия невозможна. Наконец, если 
энергии атомов равны, то создаются условия для образования пары: междоузельный 
атом - вакансия. Таким образом, задача оценки энергии активации сводится к опреде- 
лению реальной энергии парного или иного вида взаимодействия диффундирующего 
атома и атомов окружающей его среды. 

ОБОСНОВАНИЕ МЕТОДИКИ 


Теоретический расчет энергии активации в такой постановке достаточно затрудните- 
лен. Поэтому различные авторы связывают ее с некоторыми другими величинами, ко- 
торые можно рассчитать или экспериментально установить. Согласно данным работы 
[9], при самодиффузии энергия активации должна быть близка по величине к теплоте 
испарения кристалла. Однако на опыте наблюдаются более низкие ее значения. В ос- 
новном энергию активации связывают с температурой плавления химических элемен- 
тов [8]. Это пригодно лишь для грубой оценки, поскольку отсутствует критерий выбора 
эмпирической константы в приводимых уравнениях. Для более сложных систем леги- 
рования, чем кристаллы, расчетных методов вообще не существует. 

В работе [4] показано, что условие образования стабильных фаз для многих сис- 
тем, в том числе и сложных [9], сводится к условию равенства пространственно- 
энергетических параметров атомов. Для парного взаимодействия двух разнородных 
атомов это условие запишется: 


р'эіК\=Р," /К 2 , ( 1 ) 

где Р э и Р э - пространственно-энергетические параметры взаимодействующих атомов; 
Кі и Кг - их координационные числа. 
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При таком подходе равенство энергии этих атомов есть энергетическое условие 
образования системы: атом А - атом В - вакансия атома В , поскольку атом В имеет 
двойное энергетическое взаимодействие с атомом В за счет энергии валентных элек- 
тронов. Это как раз и дает условие образования вакансии. 

МЕТОДИКА РАСЧЕТА 


Рассмотрим бинарную систему (М -М ). С помощью Р-параметра по методике ра- 
боты [4] определяем усредненную эффективную энергию атома по его валентным элек- 
тронам ( Р 0 /2гі). Разделив это значение на число эффективных валентных электронов п, 
получим некоторую эффективную энергию атома, приходящуюся на один валентный 
электрон: 


^ ° 2 П п ' 


( 2 ) 


Эта величина является прямой характеристикой энергии активации процессов 


диффузии и самодиффузии в системах типа М — М 

Используя ранее установленные в [4,9] зависимости с использованием Р- 
параметров, с учетом (2) получим результирующую величину эффективной энергии 
парного взаимодействия атомов А и В тождественную энергии активации процесса 
диффузии: 



= 2 


[ г,п^ 


+ 


{ Г іП^ 


у Роу А \Р°УВ 


( 3 ) 


где Е а энергия активации процесса диффузии атома В в среде атомов А. 

При диффузии атома В в однородной по составу и подобной ему среде с атомами 
В значение Е а приобретает смысл энергии активации самодиффузии атома В. Исходя из 
(3) она определится как: 


где Е ас - энергия активации самодиффузии. 

По уравнениям (2) - (4) рассчитаны значения энергии активации самодиффузии 
ряда химических элементов периодической системы Менделеева Д.И. (табл.1) и их 
парною взаимодействия (табл. 2). Во всех случаях за величину п принимали или число 
всех валентных электронов, или число наиболее удаленных от ядра валентных электро- 
нов конкретного атома. Какой из этих двух вариантов, и при каких условиях реализует- 
ся в экспериментах предсказать для всех случаев пока затруднительно. Тем не менее, 
сопоставление полученных результатов с экспериментальными данными работ [1-3, 8, 
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Таблица 1 . Энергии активации объемной самодиффузии атомов 


Атом 

<2 0 (эВА) 
по (2) 

'гі(А) 

Кол-во ва- 
лентных 

электронов, 

п 

Еас (эВ) 
расчет 
по (4) 

Е ас (ЭВ) 

эксперимент 

и 

3,487 

1,586 

1 

0,55 

0,57 

Ве 

13,347 

1,04 

2 

1,60 

1,7-1,63 

с 

51,739 

0,596 

4 

7,23 

7,07 

М§ 

15,436 

1,279 

2 

1,51 

1,40 

А1 

31,624 

1,312 

3 

1,47 

1,48 

8і 

54,394 

1,068 

4 

4,24 

4,76 

2п 

7,623 

1,065 

2 

0,896 

0,885-0,95 

2п 

12,798 

1,065 

2 

1,50 

1,34 

С(3 

8,349 

1,184 

2 

0,881 

0,83 

Ое 

61,175 

1,090 

4 

3,508 

3,15 

а2г 

58,955 

1,593 

4 

2,31 

2,25 

Р2г 

46,883 

1,093 

4 

1,84 

1,65 

РТІ 

39,114 

1,477 

4 

1,66 

1,52 

V 

71,579 

1,401 

5 

4,26 

4Д98 

Те 

50,542 

1,111 

6 

1,896 

17,5-2,03 

Ш 

21,616 

1,476 

2 

1,83 

1,68 

8п 

20,761 

1,240 

4 

1,05 

1,01-1,06 

8Ь 

41,87 

1,193 

5 

1,76 

1,55-2,08 

Ре 

29,026 

1,227 

3 

2,957 

2,704-3,202 

№ 

46,954 

1,139 

4 

2,599 

2,90* 

Си 

20,353 

1,191 

3 

2,204 

2,1-2,33 

А§ 

18,044 

1,285 

2 

1,754 

1,91* 

Сг 

70,674 

1,453 

5 

2,424 

2,27* 

Мо 

123,44 

1,520 

4 

5,076 

5,18* 

Та 

88,95 

1,413 

3 

5,245 

4,748* 

\Ѵ 

126,77 

1,360 

6 

3,883 

3,088-4,616* 

N6 

85,555 

1,589 

3 

4,4868 

4,532* 

2г 

58,773 

1,593 

3 

3,0829 

3,487* 

- - по данным райоты р] 
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Таблица 2. Энергия активации объемной диффузии элементов 


в кремнии (8і): Т пл =1688 К; орбиталь ЗР 1 с Р 0 =6,716 эВА; п=1; 
орбиталь ЗР 2 +38 2 Р 0 =54,394 эВА; п=4 

Параметры 

диффундирующего элемента 

диффузионной среды 

Е а , эВ 

Е а , эВ 

Атом 

ГК 

Р 0 , эВА 

п в 

Р т , эВА 

Р 0 , эВА 

п А 

Р т , эВА 

расчет 

эксп-нт 

Іл 

998 

3.487 

1 

- 

6.716 

1 

5.55 

0.77 

0.62-0.79 

А1 

1523 

3 1 .624 

1 

31.624 

54.394 

4 

52.868 

4.09 

3-3.77 

С 

1688 

51.739 

4 

35.627 

54.394 

4 

- 

3.44 

2.9 

8Ь 

1688 

103.85 

5 

121.94 

54.394 

4 

- 

3.92 

3.92 

в германии: орбиталь 4Р 1 , Р 0 = 7,128; п=1; ц = 1,090 А; орбиталь 4Р" -Р 0 = 19,361 эВА; г, 
= 1,090 А; п=2; орбитали: 48 2 + 48 2 - Р 0 =61,17 эВА; п= 4; г тах 1,090 А . 

Іл 

- 

3,487 

1 

- 

7,128 

1 

- 

0,469 

0,46 

2п 


19,599 

12,798 

2 

2 


61,176 

61,176 

4 

4 

— 

2,78- 

2,10 

2,80-2,16 

А1 


31,624 

3 

- 

61,176 

4 

- 

2,55 

2,70 

Іп 


40,749 

3 

- 

61,176 

4 

- 

2,96 

3,2 


в никеле (№): " 

Г пл = 1713 К; орбиталь 48 1 +48 1 с 

Ро= 18,838 эВА; п=2 

Си 

683 

13,272 

2 

11,06 

18,838 

2 

13,605 

2,233 

1,539* 

в вольфраме (\Ѵ): Т пл = 3473 К; орбиталь 5с1 +50 с Р 0 =61,959 эВА; п= 

2 

Мо 

1806 

30,848 

4 

- 

61,959 

2 

- 

3,135 

3,474* 

2г 

2000 

18,547 

2 

- 

61,959 

2 

- 

3,55 

3,367* 


в алюмі 

ании (А1): 1 пл - 933 К; орбитали ЗР 1 ;38* ;38‘ с Р 0 = 31,624 эВА; п=2 

!Аё__ 

“ 

7,717 

1 

- 

17,451 

2 

- 

1,550 

1,413** 

Си 

“ 

13,272 

2 

- 

31,624 

3 

- 

1,610 

1,413** 

8і 

- 

6,716 

1 

- 

31,624 

3 

- 

1,740 

1,322** 


в 

меди (Си) 

I пл = 1356 К; орбиталь 36' 

с Р 0 = 13,272 эВА; п=2 


А§ 

“ 

23,558 

2 

- 

13,272 

2 

- 

1,691 

1,612** 

Мп 

“ 

18,025 

2 

- 

7,081 

1 

- 

1,597 

1,005** 

8п 

” 

20,761 

2 

- 

13,272 

2 

- 

1,672 

1,352** 


В уге: 1 ПЛ (Ее)- 1805 К; орбитали 48 ;48 с Р 0 = 18,462 эВ 

А; п=2 


С 


51,739 

4 

- 

29,026 

3 

- 

2,899 

1,586** 

Си 

“ 

13,272 

1 

- 

29,026 

3 

- 

2,309 

2,644** 

Мп 

“ 

74,043 

5 

- 

29,026 

3 

- 

2,799 

2,860** 

* - по данным работы [5 

; ** - по данным работы [15] 
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10-13] (табл. 1 и 2) показало их удовлетворительное совпадение и достоверно можно 
утверждать, что при парном взаимодействии атомов чаще всего п равно максимально- 
му числу валентных электронов. 

Диффундирующий элемент В и окружающая его среда атомов А могут находить- 
ся в одном фазовом состоянии (например в твердом), либо в разных. Учет фазового со- 
стояния элементов в процессе диффузии осуществляется через эмпирическое уравне- 
ния: 

Р т =Цх± к ~ Т '^\ ( 5 ) 

г Д е Т — температура, при которой протекает процесс диффузии, К; Т т — температура 
плавления. К; к - коэффициент равный 1 0 для твердых тел. 

Если атомы среды находятся в жидкой фазе, а диффундирующий элемент имеет 
температуру плавления выше, то, а = 0,9 и после единицы ставится знак «-». Т т - тем- 
пература плавления диффундирующего элемента, а Т — температура плавления окру- 
жающей его среды. Если в твердое тело диффундирует элемент, находящийся в жидкой 
фазе, то а= 1,1; знак «+» после единицы. 

С учетом (5) уравнение (3) запишется: 



= 2 


Г г, п л 


+ 


Г г 


\Рт)а \Рт) в 


(6) 


Результаты расчетов по уравнению (6) на примере объемной диффузии в кремнии 
приведены в табл. 2, которые достаточно хорошо согласуются с экспериментальными 
данными. 

Рассмотрим оксидные и интерметаллидные системы. На основании правил, выве- 
денных в работе [4], определим энергию межатомного взаимодействия в системе ато- 
мов А и В по уравнению: 

1 _ 1 1 

в А + О В ’ (7) 


где <2 а = 




\ г і'А 


; Я в = 




У г і1в 


<2а <2в - эффективные энергии атомов элементов А и 


В, входящих в систему; - энергия системы. 

Далее находим общую эффективную энергию активации диффузии атома В в 
системе “диффундирующий атом В - диффузионная среда атомов А ” по уравнению: 
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Таблица 3. Энергия активации объемной самодиффузии и диффузии в оксидах 


Диффузионная 

среда 

Диффундирующий 

элемент 

Энергия активации 
диффузии 

Систе- 

0. 

Атом 

Ро 

Яв 

Ов 

п 

Еаэ (ЭВ) 

Е аъ (эВ) 

ма 

(эВ) 


(эВ) 

(эВ) 

(эВ) 


расчет 

эксперим. 

СаО 

7,6415 

Са 

15,803 

4,676 

9,3509 

2 

2,90 

2,71 

СаО 

7,6415 

О' 2 

17,304 

6,362 

9,3509 

2 

3,47 

3,50 

ВеО 

9,819 

Ве 

13,347 

6,417 

12,834 

2 

3,88 

4,0 

М§0 

9,365 

м ё 

15,436 

6,035 

12,068 

2 

3,67 

3,43 

8гО 

7,7135 

8г 

17,367 

4,730 

9,459 

2 

2,93 

2,76 

ЗпО 

11,954 

5п 

20,761 

8,372 

16,743 

2 

4,92 

5,48 

РЬО 

14,292 

РЬ 

20,936 

3,573 

21,719 

4 

2,69 

2,86 

СоО 

9,425 

Со(ІІ) 

14,372 

1,521 

12,169 

8 

3,8 

4,1 

СоО 

9,425 

О' 2 

17,304 

6,362 

12,169 

2 

3,8 

4,1 

№0 

11,850 

№(ІІ) 

18,838 

2,067 

16,539 

8 

1,76 

1,92 

№0 

11,850 

О' 2 

17,304 

3,181 

16,539 

4 

2,51 

2,50 

РеО 

9,063 

Ре(ІІ) 

18,462 

1,447 

11,573 

8 

1,25 

1,29 

СГ2О3 

13,903 

Сг(ІІІ) 

30,27 

6,944 

20,833 

3 

4,63 

4,34 

СГ2О3 

13,903 

О' 2 

17,304 

6,362 

20,833 

2 

4,36 

4,38 

Ре 2 0 3 

15,351 

Ре(ІІІ) 

29,026 

8,087 

24,262 

3 

5,3 

4,86 

Ре 2 0 3 

15,351 

О' 2 

17,304 

6,362 

24,262 

2 

4,50 

4,81 

ЫцАІз 

7,110 

А1(И) 

17,469 

13,099 

6,549 

2 

2,508 

2,028* 

№А1 

2,577 

N1(1) 

6,708 

5,889 

5,889 

1 

1,793 

1,769* 

№ 3 А1 

4,722 

N1(1) 

6,708 

5,889 

1,963 

3 

1,390 

0,5179* 

РезАІ 

4,659 

Ре(І) 

7,098 

5,784 

1,928 

3 

1,363 

1,725** 

РеАІ 

12,809 

Ре(І) 

7,098 

5,784 

5,784 

1 

2,556 

2,080** 

РеАЦ 

3,461 

Ре(ІІ) 

18,462 

15,046 

5,015 

3 

2,047 

1,864** 

РеА1 3 

4,362 

Ре(ІІІ) 

29,026 

23,656 

7,885 

3 

2,758 

1,634** 

РезАІ 

4,659 

АІ(ІІІ) 

31,624 

23,939 

7,979 

3 

2,922 

2,189** 

РеАІ 

12,809 

А1(І) 

6,055 

4,583 

2,292 

1 

1,943 

1,088** 

РеАЬ 

3,461 

А1(І) 

6,055 

4,583 

2,292 

1 

1,836 

2,388** 

РеА1 3 

4,362 

А1(І) 

6,055 

4,583 

2,292 

1 

1,502 

1,214** 
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Пример расчета энергии активации объемной диффузии Ре в оксиде Ре 3 0 4 
Ре(ІІ) (}і = 1Р573; (}Г= Р447; Ре(ІІІ) (^ = 24.262; <5Г= 8.073; О(ІІ) (^= 4Р797 
Оксид Ре 3 0 4 есть сложная система двух компонентов РеО и Рег 0 3 . Поэтому рассчи- 
таем СЬ как двухатомную связь каждого компонента: Ре(ІІ)-0 С> 5 = 9.063; Ре(ІІІ)-0 
0 5 = 24.262 и всей системы: 0 5 = 5,699 

Рассчитаем СЬ при диффузии Ре(ІІІ) в системе Ре 3 0 4 С> 5 = 3,341 
Рассчитаем С> 5 при диффузии Ре(ІІ) в системе Ре 3 0 4 (} 3 = 1,154 

При одинаковом энергетическом вкладе этих атомов железа энергия его активации 
определится как средняя величина: 3,341 + 1,154=2,25 эВ; экспериментальная величина 
составляет по данным [14] 2.12-2.38 эВ 
* - по данным работы [16]; ** - по данным работы [17] 


1 = 1 1 

Е аэ ^ ^ ^ 


где Е аэ - эффективная энергия активации диффузии диффундирующего в системе ато- 
ма; бв - минимальное значение эффективной энергии диффундирующего атома В, оп- 
ределяемое по уравнению: 




\ 2 г і п 


7 В 


■ІѴ, 


(9) 


где М- число элементов в системе. Для бинарных оксидов и интерметаллидов N=2. На 
примере диффузии атомов железа в сложном оксиде Ре 3 0 4 в табл. 3 приведен расчет 
эффективной энергии активации данного диффузионного процесса. При этом у ионов 
кислорода использовалось значение ионного радиуса, а не орбитального. Величина п 
для атомов окружающей среды принимается равной 1 . Сопоставление эксперименталь- 
ных [14] и расчетных значений энергии активации удовлетворительно согласуются ме- 
жду собой (в пределах погрешности эксперимента), а учет температурных факторов 
может еще более улучшить это согласование. 

Таким образом, расчет энергии активации атома в сложной системе в принципе 

аналогичен расчету простой системы М — М . 
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ВЫВОДЫ 

1. Основываясь на данной методике, можно проводить расчеты и коэффициентов 
диффузии атомов отдельного элемента в сложных системах легирования. Особенно 
это актуально при диффузионной сварке разнородных металлов и сплавов, так как 
от величины их взаимного диффузионного проникновения друг в друга зависят 
свойства получаемых сварных соединений. 

2. Поскольку результирующая величина эффективной энергии любой сложной струк- 
туры может быть сравнительно просто рассчитана по принципам применения Р- 
параметров компонентов ее составляющих, то данный разработанный метод может 
найти применение для прогнозирования и расчета энергии активации процессов 
диффузии атомов не только в простых, но и в сложных системах твердых тел и ма- 
териалов в различных технологических процессах. 
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НА ПАРАМЕТРЫ ПАРОГАЗОВОЙ СМЕСИ 
ПЛАЗМОГАЗОДИНАМИЧЕСКОЙ УСТАНОВКИ 

А.П. БЕСОГОНОВ, Ю.Ф. КИСАРОВ, Л.В. ШИШКИНА 


Институт прикладной механики Уральского отделения Российской Акаде- 
мии наук, Ижевск, Россия 


АННОТАЦИЯ. Предложена одномерная математическая модель расчета процесса 
неравновесной конденсации металла в сверхзвуковом сопле Лаваля. Установлено, 
что давление в камере плазмогазодинамической установки на количество конденсата 
влияет очень слабо, а увеличение содержания паров металла в исходной парогазовой 
смеси приводит к уменьшению концентрации конденсата в фиксированных сечениях 
сопла. 


Как выяснилось в процессе исследований, наночастицы металлов обладают чрезвычай- 
но интересными с практической точки зрения и отличными от субмикронных частиц 
свойствами. Получение веществ в нанодисперсном виде становится важной технологи- 
ческой задачей. Данная работа посвящена численному моделированию процесса обра- 
зования металлических наночастиц в сопле плазмогазодинамической установки, соз- 
данной для получения кластеров и нанесения нанопокрытий на различные поверхности 
[ 11 - 

УРАВНЕНИЯ ГАЗОВОЙ ДИНАМИКИ 

Расчет начинается с подготовки граничных условий к системе уравнений газовой ди- 
намики. Находим безразмерную температуру в критическом сечении. 

Г кг 

Х кг=-у- ( 1 ) 

і 0 

Здесь индекс кг относится к параметрам в критическом сечении сопла, а 0 - к 
параметрам торможения потока. 

По следующим формулам находим остальные параметры смеси 
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_к_ 

Ркг ~ Ро ^ к г к-\ ’ 


Ркг - 


V 


кг 


К 0 Т кг 


5 




и расход по формуле 

0. = Ркг Ѵ кг Р кг ’ 


( 2 ) 


( 3 ) 


где Р - площадь сечения сопла, р - плотность смеси, ѵ - скорость, к - показатель адиа- 
оаты, р 8ит - суммарный молекулярный вес, Ко - универсальная газовая постоянная. 

На протяжении всех дальнейших расчетов расход остается постоянным. По за- 
висимости 


20 = 


Р (рѴ 2 + р 8ит ) 

0 


( 4 ) 


находим начальное условие для уравнения импульсов. Здесь р 5ит - давление смеси. За- 
тем решается система уравнений газовой динамики [2] 


+р 5ит ) 



СІ X 


а=рѵр, 


Р кит ^ ^ 
^ СІ X 


,2 ке + 1 8 


= - + I I к и (Т) , 

2 7=1 /=1 


Р вит 


-РКоТ 


ке 


ОС: 


5іі т 


Рвит ’ 

( = 1 Рі 


Р,= а ‘ рК ° Т и , і“1*кг, 
/ Рі 

Р = Р(Х) , 


( 5 ) 


где рі - давление і -го компонента смеси, Н о- полная энтальпия, р Г суммарный моле- 
кулярный вес /-го компонента, а массовая концентрация /-го компонента смеси, 
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К) - коэффициенты в полиномах для расчета полной энтальпии, ке- число газообразных 
компонентов смеси. 

Конденсация моделировалась с помощью соотношений для мономолекулярной 
теории, которые могут быть записаны следующим образом: 

А в +Л і— ( 6 ) 


А м + А, <-$- - А 

где Л е - частица, содержащая # мономеров; К/, - константы равновесия прямых и 

обратных реакций. 

Тогда, как и в химической кинетике, уравнения для роста-уменьшения кластеров 
запишутся в виде [3]: 


д 1± 

ді 


4 4+і 


( 7 ) 


4 = К п \ 4-і - к ~* 4 = К 1 п \ Ѵі 


4-і 


ѵѴі 


/, 


л 


8 

П 8 ) 


( 8 ) 


где -скорость образования #-ой фракции в единицу времени,^ - функция распреде- 
ления частиц ^-фракции, п & - равновесная функция распределения кластеров по разме- 
рам. 

Система уравнений (7) неудобна при численном анализе, поэтому сделаем заме- 
ну переменных: 


§ ]+1 =§ ] +А к , А к = А к _ х -тоё, Д,=1, к = 

Чтобы свести бесконечно большое число уравнений в системе (7) к ограничен- 
ному все пространство размеров делилось на группы со своим шагом Л к , который уве- 
личивался по мере возрастания номера группы в несколько раз. Коэффициент увеличе- 
ния шага тосі выбирался в зависимости от затрат времени на расчет и соответствия 
экспериментальным данным. При такой замене происходит сжатие по переменной §: 
сначала при у— 1 10 учитывается каждая фракция, затем они объединяются по две, да- 
лее по четыре и т.д. до у- 100. Количество объединяемых фракций можно изменить. Да- 
лее предполагаем, что между расчетными точками у и у'+У функция распределения по- 
стоянна и равна значению в точке у+У. 
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Краевые условия к системе (7) ставятся следующим образом: 



( 10 ) 


где а р - суммарная доля неконденсирующихся составляющих смеси, а соп - массовая 
доля конденсата. 

Массовая доля конденсата рассчитывалась нами по формуле: 



(П) 


І= 2 к=8і 


Уравнения газовой динамики (5) для сверхзвуковой части сопла решались мар- 
шевым методом, уравнения (7) - методом прогонки. 

РАСПРЕДЕЛЕНИЕ ПАРАМЕТРОВ В ЗАВИСИМОСТИ ОТ СОСТАВА 
РАБОЧЕГО ТЕЛА 

Были проведены расчеты для различных соотношений аргона и железа в рабочем теле 
[4]. Массовая доля железа ос р е менялась от 0,2 до 0,8. Температура и давление в камере 
были Т 0 = 7000 К , ро = 0,2 МПа. Представленные кривые показывают изменение 
параметров в трех сечениях: 1 - х— 12 мм, 2 - х = 16 мм, 3 - х— 20 мм. Все сечения попа- 
дают в зону конденсации. Наиболее интересные результаты расчетов представлены на 
рис. 1 - рис.З. 

Температура в любом из сечений растет с увеличением содержания железа в ра- 
бочем теле (рис. Іа). Это, по-видимому, объясняется особенностями расширения и теп- 
лофизическими характеристиками компонентов смеси, т.к. то же можно наблюдать и 
для течения с равновесной конденсацией (рис. За), и для "замороженного" течения (рис. 
- 5В )- Распределение температуры по сечениям ( Тз >Т / ) объясняется особенностями рас- 
пределения этого параметра по длине сопла. Наиболее высокая температура наблюда- 
ется при равновесной конденсации, т.к. здесь выделяется больше тепла из-за образова- 
ния максимально возможного количества к-фазы. Затем по количеству выделяющегося 
тепла идет неравновесный процесс. Характеры графиков равновесного и неравновесно- 
го процессов совпадают, а в "замороженном" течении имеют другую кривизну, и тем- 
пература там ниже. Изменение температуры в сечении для неравновесного и равновес- 
ного течений около 54 и 250 К соответственно; оно немного растет по мере удаления 
сечения от критического, а для "замороженного" 650 К. 

Массовая доля конденсата при увеличении содержания железа в смеси падает 
(рис. 16). На первый взгляд это неправильно, однако содержание конденсата в каком 
либо сечении сопла в большей степени зависит от удаленности его от точки росы, а не 


44 


ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. Том 7, № 1 


ВЛИЯНИЕ КОНЦЕНТРАЦИИ ЖЕЛЕЗА И ДАВЛЕНИЯ В КАМЕРЕ НА ПАРАМЕТРЫ ПАРОГАЗОВОЙ 
СМЕСИ ПЛАЗМОГАЗОДИНАМИЧЕСКОЙ УСТАНОВКИ 


Т-10 3 , К « р.. ю 2 




Н-10%,1 3 н^іо", М 




Рис. 1. Зависимость температуры смеси (а), массовой доли конденсированной фа- 
зы (Ь), числа частиц в единице объема (с) и среднемассового размера час- 
тиц (б) от начального состава смеси в потоке с неравновесной конденсаци- 
ей для сечений: 1 - х=12 мм , 2 - х = 16мм, 3 - х = 20 мм 

от расстояния до критического сечения. Точка росы при увеличении массового содер- 
жания железа смещается к срезу сопла (рис. 2а), а значит, количество конденсата в вы- 
бранных сечениях должно уменьшаться. Эту особенность подтверждают и расчеты, 
проведенные для равновесного течения (рис. 36). Массовая доля к-фазы для равновес- 
ного процесса выше, чем для неравновесного. Поведение этих кривых для обоих про- 
цессов примерно одинаково. Изменение содержания конденсата в сечении для нерав- 
новесного течения около 0,02, и оно уменьшается с удалением сечения от критическо- 
го, а для равновесного - 0,0367 и примерно постоянно во всех сечениях. 
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X 1 0 ; 1,1 

6.500 


4.500 


2.259 


2.465 


2.064 



Хс _ 

Р 10І Па 


0.980 


0.580 


0.175 




Ь> 


8.000 


Рис. 2. Зависимость положения точки росы и температуры смеси в этой точке (а), 
давления в точке росы и парциального давления паров железа (Ь) от на- 
чального состава смеси 

Обратимся теперь к графикам, характеризующим конденсированную фазу. 

Число частиц в единице объема (рис. 1в) растет с увеличением содержания же- 
леза в смеси. Это может происходить потому, что, во-первых, увеличивается число 
атомов железа в раоочем теле и, во-вторых, они не успевают конденсироваться. Из- 
менение этого параметра составляет около 7,0 -ІО 21 1/м 3 и уменьшается по мере удале- 
ния от критического сечения. 
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т-юѴк 



Рис. 3. Зависимость температуры смеси (а) и массовой доли конденсированной 
фазы (Ь) для равновесного течения, температуры для «замороженного» те- 
чения (с) от начального состава смеси в потоке для сечений: 1 - х= 12 мм , 
2 - х = 1 6 мм , 3 - х = 20 мм 

Среднемассовый размер частиц нелинейно уменьшается с увеличением варьи- 
руемого параметра (рис. 1г). При этом его изменение составило 0,475 нм, и оно умень- 
шается с удалением сечения от критического. 

Точка росы, как было отмечено выше, сдвигается к срезу сопла с увеличением 
варьируемого параметра (рис. 2а). При этом ее координата в случае равновесного тече- 
ния находится несколько дальше от среза сопла из-за неточностей в определении дав- 
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ления насыщенных паров над плоской поверхностью жидкого железа. Смещение коор- 
динаты точки росы происходит на 4,9 мм. 

То же можно сказать и о температуре в точке росы, ее изменение - около 81,6 К. 

Давление смеси в точке росы уменьшается с увеличением массовой доли железа, 
потому что, во-первых, происходит расширение газа и, во-вторых, меняется характер 
расширения из-за изменения состава смеси (см. рис.1 а, рис. За, рис.Зв) 

Давление паров железа в точке росы, наоборот, растет при увеличении варьи- 
руемого параметра. Это может происходить по двум причинам: 

- увеличивается количество атомов железа, 

запаздывание конденсации в неравновесном процессе по сравнению с равно- 
весным. 

Разрыв между давлением смеси и давлением паров для равновесного процесса, 
благодаря более активному образованию к-фазы, значительно выше, чем для неравно- 
весного. Изменения давления смеси и давления паров в точке росы 383 и 162 Па соот- 
ветственно. 

Характер изменения всех графиков для точки росы неравновесного и равновес- 
ного процессов совпадает. 

Максимальное перенасыщение и минимальное расстояние от точки Вильсона до 
точки росы для этой группы расчетов наблюдалось при меньших а Ре . Количество кон- 
денсата, число частиц в единице объема, среднемассовый размер частиц в точке Виль- 
сона практически не менялись. 

РАСПРЕДЕЛЕНИЕ ПАРАМЕТРОВ В ЗАВИСИМОСТИ ОТ ДАВЛЕНИЯ В 
КАМЕРЕ 

Для выяснения вопроса о влиянии давления торможения на параметры смеси были 
произведены расчеты с различными значениями р 0 = (1 ,0-5,0)- 1 0 5 Па и одинаковыми 
температурой в камере Т 0 = 7000 К, составом смеси а Ре = 0,375. Исследуемые сечения 
были те же, что и в предыдущем случае, и они также находились в зоне конденсации. 
Некоторые результаты представлены на рис. 4 - рис. 6. 

Зависимость давления в сечениях от давления в камере здесь не помещена. По 
характеру она практически не отличается от равновесной. Но для той и другой сущест- 
вует очень слабо выраженная особенность в зоне конденсации, происходящая из-за из- 
менения состава смеси. 

Температура при увеличении давления в камере нелинейно растет (рис. 4а, 
рис. 6а) за счет выделения тепла при конденсации, и за счет изменения состава смеси. 
Для замороженного течения она не зависит от давления в камере, остается постоянной. 
Температура в разных сечениях сопла равновесного процесса больше отличается от 
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т- 1б 3 , к *-ю 2 




Н-ІО'Ѵ К-ю", м 




Рис. 4. Зависимость температуры смеси (а), массовой доли конденсированной фа- 
зы (Ь), числа частиц в единице объема (с) и среднемассового размера час- 
тиц (ё) от давления в камере установки в потоке с неравновесной конден- 
сацией для сечений: 1 - х=1 2 мм, 2 - х = 1 6 мм , 3 - х = 20 мм 

этих же параметров неравновесного. Изменения температуры для равновесного и не- 
равновесного процессов 190 и 83 К соответственно. Характер зависимостей температу- 
ры от давления в камере одинаков для равновесного и неравновесного процессов. 

Количество конденсата растет при увеличении давления в камере (рис.4б, 
рис.бо), вероятно из-за увеличения числа столкновений. Но это происходит нелинейно, 
возможно потому, что при конденсации повышается температура, в результате чего 
уменьшается число эффективных столкновений. Количество образующегося в равно- 
весном процессе конденсата более значительно, чем в неравновесном. Изменение мас- 
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т, -ІО' 3 , К ХИ0 3 ,М 

■і й 7 



Рис. 5. Зависимость положения точки росы и температуры смеси в этой точке от 
давления в камере установки 


т- іо' 3 , к 



Р- КГ 1 , Па 


а й -10 : 



1.000 3.000 

Р-ЮЧПл 


Рис. 6. Зависимость температуры смеси (а) и массовой доли конденсированной 
фазы (Ъ) для равновесного течения от давления в камере установки в пото- 
ке для сечений: 1 - х=12 мм, 2 - х = 16 мм, 3 - х = 20 мм 
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совой доли конденсата в смеси для равновесного и неравновесного процессов состав- 
ляет 1,8-10 и 6,1 • 10 соответственно. Вниз по течению вдоль сопла оно очень слабо 
растет. 

Число частиц в единице объема возрастает прямо пропорционально давлению в 
камере (рис. 4в). При увеличении давления начинают конкурировать между собой два 
процесса: 

- рост количества конденсата и как следствие увеличение температуры и паде- 
ние давления и /V, 

- увеличение давления в камере, рост давления по длине сопла и N. 

И как видим, второй процесс преобладает. Самое большое изменение этого па- 
раметра наблюдается в сечении, ближайшем к точке росы (Л N = 1 1,72 • ІО 21 1/м 3 ), за- 
тем в более удаленных сечениях оно постепенно сглаживается. 

Среднемассовый размер частиц К 43 (рис. 4г) при увеличении варьируемого па- 
раметра возрастает. Характер графиков подобен характеру кривых массовой доли кон- 
денсата и, видимо, объясняется теми же причинами. Изменение К 43 составляет около 
2,9 ■ 10 12 м и очень слабо растет с приближением сечения к срезу сопла. 

1 очка росы при увеличении давления в камере сдвигается к критическому сече- 
нию сопла (рис. 5). Изменение координаты точки росы на исследуемом участке варьи- 
руемого параметра составило 7,2 • 10 3 м. Некоторая нелинейность в возрастании тем- 
пературы в точке росы обусловлена изменением состава смеси при конденсации. Изме- 
нение температуры в точке росы 187,5 К . Для равновесного процесса кривые ведут се- 
бя аналогично. 

Максимум на кривой перенасыщения, расстояние от точки росы до точки Виль- 
сона, а также массовая доля конденсата, число частиц в единице объема, среднемассо- 
вый размер частиц возрастали с увеличением давления в камере. 

ВЫВОДЫ 

В работе исследовалось влияние на процесс конденсации начального состава смеси 
и давления в камере. Выявлено, что при увеличении ос р е количество конденсата в 
потоке уменьшается. Изменение ро не приводит к существенному изменению пара- 
метров потока в сопле. 

В заключении следует отметить, что полученные результаты достаточно хо- 
рошо согласуются как с экспериментальными данными, так и с данными, получен- 
ными другими авторами [5,6]. 
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А. В. ВАХРУШЕВ, А.М.ЛИПАНОВ 

Институт прикладной механики УрО РАН, Ижевск, Россия 

АННОТАЦИЯ. Предложен метод расчета потенциала парного взаимодействия на- 

ночастиц на основе аппроксимации результатов численных расчетов метода моле- 
кулярной динамики. На основе потенциала парного взаимодействия наночастиц 
рассчитана зависимость предела прочности монодисперсного порошкового нано- 
композита от размеров составляющих его наночастиц. 


ВВЕДЕНИЕ 

Развитие нанотехнологий обусловливает актуальность моделирования процессов ста- 
тического и динамического взаимодействия наночастиц. При этом используют методы, 
основанные как на положениях классической, так и квантовой механики в зависимости 
от целей и необходимой точности моделирования: метод эмпирических потенциалов, 
квантово - механическое моделирование из первых принципов, полуэмпирический 
подход, молекулярную динамику, метод Монте-Карло и т.д. Перечисленные методы 
позволяют рассчитывать параметры наночастиц исходя из информации на атомарном 
уровне, однако их применение к системам наночастиц требует больших вычислитель- 
ных ресурсов и времени. Поэтому необходимо развитие экономичных методов вычис- 
лений. 

Задачей данной работы являлось создание такой методики. С этой целью выпол- 
нен расчет сил парного взаимодействия наночастиц в зависимости от типа атомов, их 
составляющих, и построен потенциал взаимодействия наночастиц путем аппроксима- 
ции результатов численных расчетов. 

ПОСТАНОВКА ЗАДАЧИ 

Задача расчета парного взаимодействия наночастиц имеет два этапа: первый - расчет 
внутренней структуры и равновесной конфигурации (формы) каждой отдельной, сво- 
бодной о г взаимодействия наночастицы; второй - расчет парного взаимодействия двух 
наночастиц. 
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На первом этапе движение атомов, образующих в каждую наночастицу, опреде- 
ляется системой дифференциальных уравнений Ланжевена [1] 


ё 2 х 


N. 


ш • = Е р у + " р . (О - «іШі 

а і ]=і ш 

при краевых условиях 

х і =Хі 0 ,Ѵ і =Ѵ і0 л = 0,х 1 еО к , 


і = 1,2,.., N.. , 


(1) 

( 2 ) 


где ]У к - число атомов, составляющих каждую наночастицу; т г масса і-го атома; 
Х і 0 > х і - начальные и текущие координаты і-го атома, соответственно; К - силы меж- 
атомного взаимодействия; О к - область, занимаемая соответствующей наночастицей: 
Уо> У - начальная и текущая скорости і-го атома, соответственно; ОС. - коэффициент 

“трения” в атомарной системе; Р(і)- случайный набор сил при заданной температуре, 
задаваемый распределением Гаусса. 

Силы межатомного взаимодействия, обычно, являются потенциальными и опре- 
деляются из соотношения 




дф Ш 

д ъ> 


>і = 1,2,. 


,Ы к ,] = 1,2, 


к. 


( 3 ) 


где Ру - радиус-вектор, определяющий положение і-го атома относительно і-го атома: 

ф (Ріі) ' потенциал, зависящий от взаимного расположения всех атомов, п - число ти- 
пов взаимодействия между атомами. 

Потенциал Ф(р у ), в общем случае, задается в виде суммы нескольких состав- 
ляющих, соответствующих различным типам взаимодействия: 

Ф(Ру) = Ф сЬ + Ф ѵа + ф [а + Ф р8 + Ф ѵѵ + ф е5 + Ф ьь . (4) 

Здесь имеются в виду потенциалы: Ф сЬ - химических связей; Ф ѵа - валентных углов; 

Ф 1а - торсионных углов; Ф р8 - плоских групп; Ф ѵѵ - Ван-дер-ваальсовых контактов; 

Ф е$ - электростатики; Ф ьь - водородных связей. 

Указанные слагаемые имеют различный функциональный вид, а значения пара- 
метров для потенциалов взаимодействия определяются на основе экспериментов (кри- 
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сталлографических, спектральных, калориметрических и др.) и квантомеханических 
расчетов [1]. 

Задав начальные координаты (а значит, и силы взаимодействия атомов) и скоро- 
сти всех атомов каждой наночастицы в начальный момент времени, согласно уравне- 
ниям (2), находим из решения уравнений движения (1 ) изменение координат и скоро- 
стей атомов каждой наночастицы во времени. Так как на наночастицы не действуют 
внешние силы, то с течением времени они принимают некую равновесную конфигура- 
цию атомов, которую будем использовать для следующего этапа расчетов. 

На данном этапе решения задачи рассматриваем две взаимодействующие нано- 
частицы, расположенные друг от друга на расстоянии 8 (рис.1). В этом случае уравне- 
ние (1) примет вид: 


ГПІ 


сГХ| 

ёі 2 


Ѵ , 2 


= I г + цц 


1=1 


СІХ, 

«іПІі ТГ’ 

СЦ 


і = 1 , 2 ,..,(ІЧ 1 + М 2 ), 


( 5 ) 



Рис.1. Схема взаимодействия наночастиц; А - увеличенное изображение фраг- 
мента наночастицы 
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при краевых условиях 



( 6 ) 


где N,,N 3 - число атомов в первой и второй наночастицах, соответственно; О, О г 
области, занимаемые первой и второй наночастицами, соответственно. 

Решение (5) при краевых условиях ( 6 ) позволяет рассчитать траектории движения 
атомов каждой наночастицы, а, следовательно, и наночастиц в целом. При этом сум- 
марные силы взаимодействия между частицами будут определяться соотношением 


РЫ=-РЬ2=]Г ^Ру, 


(7) 


где і, .і - атомы первой и второй наночастицы, соответственно. 

МЕТОДЫ РЕШЕНИЯ 

Сформулированная в предыдущем разделе задача, в общем случае, не имеет аналитиче- 
ского решения, поэтому, как правило, используют численные методы [1]. В данной ра- 
боте численное интегрирование уравнений движения атомов наночастицы в процессе 
релаксации проводилось по схеме Верлета [2]: 



( 8 ) 



(9) 



(Ю) 


где Х| , V, - координаты и скорость /-го атома на п- ом шаге по времени; Ді - шаг по 
времени; индекс “ 0 ” обозначает параметры в момент времени і= 0 . 
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РЕЗУЛЬТАТЫ РАСЧЕТОВ 

Реализацию предложенной выше методики рассмотрим на примере расчета взаимодей- 
ствия частиц алюминия. 

В расчетах использовался потенциал взаимодействия атомов Морзе[ 1]: 

Фщ =О т (ехр(-2Х т (р у -Ро ))-2ехр(-А, т (Ру -р 0 ))) (12) 

В соответствии с (12), силы межатомного взаимодействия вычисляются как 
Ру = 2Е) т^т (~ехр(-2А, т (ру -р 0 )) + ехр(-А. т (ру - р 0 ))), (13) 

где Э т Д т ,р 0 - константы материала. 

На первом этапа задачи в качестве начальных принимались координаты атомов, 
размещенных в узлах кристаллической решетки макроматериала (рис.2(1)). В процессе 
релаксации, согласно расчетам по уравнениям (6)-(9), исходная система атомов пере- 
страивается в новую “равновесную” конфигурацию (рис. 2(2)), удовлетворяющую ус- 
ловию приближения к минимуму потенциальной энергии системы (рис. 2, график). 



Рис.2. Исходная кристаллическая (1) и кластерная (2), после релаксации, струк- 
туры наночастицы из 1331 атома; график изменения потенциальной 
энергии этой системы атомов в процессе релаксации 
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При вычислении сил взаимодействия двух наночастиц одинакового размера ис- 
пользовались их параметры после свободной релаксации. 

На рис. л представлены результаты расчетов, показывающие влияние размеров на- 
ночастиц на силу их взаимодействия. Из графика видно, что более крупные наночасти- 
цы притягиваются сильнее, т.е. максимальная сила взаимодействия увеличивается с 
ростом размера частицы. Поделим силу взаимодействия наночастиц на максимальное 
ее значение для каждого размера наночастицы, соответственно. Полученный при этом 
график “относительной “ (безразмерной) силы (рис. 4) показывает, что данная величина 
практически не зависит от размера наночастицы, так как все кривые сближаются и мо- 
гут быть аппроксимированы одной линией. 

На рис. 5 приведена зависимость максимальной силы притяжения между наноча- 
стицами от их диаметра, которая характеризуется нелинейностью и общей тенденцией 
роста максимальной силы с увеличением размера наночастиц. 

Суммарная сила взаимодействия наночастиц определяется произведением двѵх 
приведенных графиков (рис. 4 и рис. 5). 



6.СЕ-08 

4.СЕ-08 

2.СЕ-08 

О.0Е+О0 

-2Д1Е-08 

-4.СЕ-0Й 

-8.СЕ-08 

-8.СЕ-08 
-1.СЕ-07 -1 


Рис. 3. Зависимость силы взаимодействия Р наночастиц от расстояния 8 между 

ними и размера частицы: 1-6=2,04; 2-6=2,40; 3-6=3,05; 4-6=3,69; 5-6=4,09 
(нм) 


58 


ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. Том 7, № 1 


РАСЧЕТ ПОТЕНЦИАЛА ПАРНОГО ВЗАИМОДЕЙСТВИЯ НАНОЧАСТИЦ 



Рис. 4. Зависимость “относительной ’ силы взаимодействия Р наночастиц от 
расстояниями 8 между ними 
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Рис.5. Зависимость максимальной силы притяжения Р тах наночастиц от диаметра 
СІ наночастицы 
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Используя полиноминальную аппроксимацию кривой на рис. 4 и степенную - кри- 
вой на рис. 5, получим 


Р = (-1 . 133 6 + 3.088' - 3.418 4 - 0.588'' + 0.828 - 0.00335)1 о\ (14) 


= 0.5 -ІО -9 -О 14 ". 


(15) 


где сі, 5 -диаметр наночастиц и расстояние между ними в нанометрах (нм), соответст- 
венно; Р шах - максимальная сила взаимодействия наночастиц в ньютонах (н). 

Следует отметить, что предлагаемые аппроксимации ограничены размерами на- 
ночастиц от 2 до 10 нм. Это объясняется тем, что для частиц менее 2нм возможно слия- 
ние частиц в одну, а верхний предел аппроксимации обусловлен изменением формы 
наночастицы при увеличении ее размера, более 10 нм [3-5]. 

С использованием зависимостей (14)-(16) был рассчитан предел прочности нано- 
композита при различных схемах "‘упаковки” наночастиц в композите (рис. 6). 

На рис. 7 представлена зависимость предела прочности нанокомпозита, сформиро- 
ванного из монодисперсных наночастиц, от их размеров. Видно, что с уменьшением 


размера наночастиц прочность наноматериала увеличивается и наоборот. Расчеты по- 
казали, что на прочностные характеристики нанокомпозита существенно влияет тип 
"упаковки” наночастиц в материале: прочность материала растет при увеличении плот- 
ности упаковки наночастиц. Особо отметим, что прочность материала изменяется об- 


ратно пропорционально диаметру наночастицы в степени 0.5, что соответствует экспе- 
риментально установленному закону изменения прочности наноматериалов (закон 
Холла - Петча)[6]: 



Рис.6. Различные типы “упаковок” наночастиц в композите 
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РАСЧЕТ ПОТЕНЦИАЛА ПАРНОГО ВЗАИМО ДЕЙСТВИЯ НАНОЧАСТИЦ 
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Рис. 7. Зависимость предела прочности нанокомпозита, сформированного из мо- 
нодисперсных наночастиц, от размера наночастиц 

ст = С-с! (17) 

где С=С тах =2.17 'ІО 4 - максимальная плотность упаковки; С= С т і п =6.4 'ІО 3 - мини- 
мальная плотность упаковки. 

Работа поддержана Российским фондом фундаментальных исследований. Про- 
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УДК 541.11 


МЕТАСТАБИЛЬНЫЕ САМОАССОЦИАТЫ (МАЛЫЕ КЛАСТЕРЫ) 

В СИСТЕМЕ МЕДБ+ АРГОН 

Г.К. МОИСЕЕВ 

Российская академия наук, Уральское отделение, ГУ Институт металлур- 
гии, Екатеринбург, Россия 


АННОТАЦИЯ. Рассчитаны термодинамические свойства и функции метастабиль- 
ных самоассоциатов Сиг, Сиз, Сиб. Методами термодинамического моделирования с 
использованием модели идеальных растворов продуктов взаимодействия рассчита- 
ны составы конденсированной и газовой фаз при Р = ІО 5 Па и Т = 400-3000 К с уче- 
том существования малых кластеров. 


1. ВВЕДЕНИЕ 

Возможность существования летучих кластеров Си 2 и Си 3 показана в работах [1-3]. Ве- 
личины АН 29 8 конденсированного кластера Си 2 (СЭО(к*Си 2 )) оценены в [4]. Согласно 
исследований [5], выполненных методами молекулярной динамики, наиболее вероят- 
ными малыми кластерами меди являются к*Си 2 , Си 3 и Сиб. Для этих частиц при ис- 
пользовании разных геометрических представлений о структуре кластеров и методов 
молекулярной динамики рассчитаны энергии связи между атомами в изолированных 
кластерах (Еь(к*Си п ), эВ/атом при п = 2, 3 и 6). Данных о других метастабильных кла- 
стерах нами не обнаружено. Для летучих атомов (Си,) и димеров (Си 2 ) полные термо- 
динамические характеристики имеются в [1], также как и для конденсированных ато- 
мов меди (к*Сиі). Для остальных упомянутых частиц полная термодинамическая ин- 
формация отсутствует. 

Цели настоящей работы были следующими: - рассчитать термодинамические 
свойства и функции летучего кластера Си 3 , конденсированных кластеров Си 2 , Си 3 и 
Сиб; - с участием атомов и перечисленных кластеров выполнить равновесные компью- 
терные эксперименты по определению состава меди и газовой фазы в широком интер- 
вале температур. 
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2. МЕТОДЫ ИССЛЕДОВАНИЯ 

2Л. Определение термодинамических характеристик 

Летучий кластер Сщ. С использованием СЭО(Сиі и С112), равных соответственно 337.6 
и 478.96 кДж/моль [1], и методики [6] построена зависимость 

1§(СЭО(Си п ))~ 2.3765 + 0.1519-П, п=1-3, (1) 

1 Де СЭО(Сііп) в кДж/моль и п — число атомов в кластере. Подобным образом построены 
зависимости 


298 (СиД) ~ 2.05853 + 0.1 62335'П, п= 1-^-3, , (2) 

1§(Н 298 "Н д (Си„)) ~ 3.58662 + 0.20558-п, п = КЗ (3) 

где 8 29 8 и Н 2 98 -Нр в Дж/(К-моль) и Дж/моль, соответственно. Для оценки С Р (Т) ис- 
пользована аддитивная схема [7]. Проверочные расчеты показали, что значения 
Ср(Т)(Си п ) при 298 К и 2000 К изменяются закономерно: 


т,к 

п= 1 

2 

3 

Ср(Т) 298 

20.74 

36.58 

54.87 

Дж/(К-моль) 2000 

21.53 

38.85 

58.27 


Результаты оценки термодинамических свойств Сиз, температурные зависимо- 
сти приведенной энергии Гиббса приведены в табл. 1. 

Конденсированные кластеры. Оценку СЭО(к*Си„) при п = 2. 3 и 6 выполнили с исполь- 
зованием уравнения [4] 

АН 298 ~ 5940. 1 •ггЕ ь (к*Си п ), Дж/моль, (4) 

где величины Е ь (к*Си п ) (в эВ/атом) взяты как среднеарифметические из максимальных 
значений для каждого п по данным [5], где варьировалась геометрия кластера и мето- 
дика расчета. При п = 2, 3 и 6 число вариантов соответствовало соответственно 4, 8 и 
20. Выбор максимальных значений Еь(к*Си п ) обусловлен выбором наиболее прочных 
связей между атомами в кластерах с конкретными значениями п. 

Для расчета 8 298 » Н 298 ~Н о , АН плавления , Ср(Т) использованы методы, описанные 
в [6-8]. За температуру плавления кластеров принята Т пл меди = 1357 К [3]. Принятые 
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Таблица 1. Принятые за достоверные термохимические свойства метастабильных кластеров меди 


МЕТАСТАБИЛЬНЫЕ САМОАССОЦИАТЫ (МАЛЫЕ КЛАСТЕРЫ) В СИСТЕМЕ МЕДБ+ АРГОН 
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за достоверные термодинамические свойства конденсированных кластеров меди и их 
температурные зависимости приведенной энергии Гиббса даны в табл. 1. 

2.2. Методика компьютерных экспериментов 

Использованы методология термодинамического моделирования (ТМ), модель идеаль- 
ных растворов продуктов взаимодействия (ИРПВ) для описания конденсированной фа- 
зы, БД АСТРА.ВА8 и БД АСТРА. 0\ѴЫ [9, 7]. Исходные рабочие тела имели состав 
(масс. %): 99% Си + 1% Аг. ТМ выполнено при Р = 10 5 Па, Т = 400-3000 К с шагом по 
температуре 10-100°. В составе раствора ИРПВ учитывали атомы меди, конденсиро- 
ванные кластеры Си 2 , Си 3 и Си 6 ; в газовой фазе кроме Аг, атомы (Сщ), частицы Си 2 , 
Си 3 , Си +1 , е - газ. 

3. РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ 

Составы конденсированной меди в твердом и жидком состоянии, как уже указывалось, 
были рассчитаны с использованием модели ИРПВ (см. рис. 1, А и Б). Из рис. ІА видно, 
что при 400-1200 К твердая фаза может содержать, кроме атомов, метастабильные кла- 
стеры. При 400 К, т.е. при комнатной температуре, общее содержание кластеров равно 
0-353 мол. %, при 1200 К ~ 15 мол %. Таким образом, можно предполагать, что в тем- 
пературной области твердого (или кристаллического) состояния (до 1357 К) существу- 
ют кристаллы из атомов меди, как основы твердого раствора, и кристаллы меди из ди- 
меров, тримеров и бинарных (двойных) тримеров меди. Возможность образования кри- 
сталлов меди из этих частиц представляется реальной в рамках принятых представле- 
ний о структуре меди, т.е. плотных кубических и гексагональных упаковок [10]. Учи- 
тывая недостаточно высокие возможности идентификации методами РФА и РСА по- 
добных кристаллов, статистически равномерно распределенных в основной структуре, 
образованной только атомами меди, можно предполагать, что медь в кристаллическом 
состоянии содержит кристаллы из атомов и малых кластеров меди. 

Увеличение температуры сопровождается ростом содержания кластеров меди, 
при всех температурах наблюдается иерархия: 

х[к*Сщ] > х[к*Си 2 ] > х[к*Си 3 ] » х[к*Си 6 ]. 

Численно это выглядит следующим образом: 

х[к*Си п ], мол. % 
п 12 3 

400 К 99.647 0.318 0.0348 

1200К 84.955 10.69 4.343 


6 

1.74- ІО' 9 
0.0107. 
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Рис. 1 . Содержание (мол. %) атомов и метастабильных кластеров в конденсиро- 
ванной фазе меди в зависимости от температуры. 

А. 400-1200 К (твердая фаза). Б. 1400-2800 К(жидкая фаза). 

1-3 - х[к*Сиі - к*Сиз]; 4 - х[к*Сііб] в жидкой фазе. 
х[к*Си 6 ] в твердой фазе = 1.74-10' 9 (400 К) -0.0107 (1200 К) 
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Приблизительная оценка количеств (или чисел) плотных кубических упаковок 
из атомов, димеров, тримеров и двойных гримеров, образованных из 1 моля исходных 
атомов меди с учетом данных ТМ при 1200 К, показывает, что на 1 ячейку из двойных 
тримеров образуется -810 ячеек из тримеров, - 3000 ячеек из димеров и - 47700 ячеек 
из атомов. Таким образом, при равномерном распределении "кластерных" кристаллов 
по объему кристаллической фазы меди возможность идентификации "кластерных" кри- 
сталлов даже при 1200 К очень невелика. Поскольку РФА и РСА приводят, как прави- 
ло, при комнатных температурах (или близких к ним), то идентифицировать существо- 
вание кластерных кристаллов технически невозможно. Отметим также, что, исходя из 
результатов ТМ, кластерные кристаллы будут образовываться преимущественно из 
димеров и тримеров. 

В результате ТМ установлено, что температура кипения модельного расплава - 
2900 К, АН кипения — 288 кДж/моль, что удовлетворительно согласуется с данными раз- 
ных авторов (2845 К и 300,5 кДж/моль по [3]; 2816 К и 302 кДж/моль по [11]). 

Для области жидкого состояния (рис. 1Б) характер изменения х[к*Си п ] = і'(Т) и 
иерархия содержания компонентов модельного расплава подобны рассмотренным для 
области кристаллического состояния и характеризуются следующими показателями: 

х[к*Си п ], мол. % 

п 1 2 3 6 

1400 К 81.44 12.934 5.6 0.0268 

2800 К 68.564 20.767 10.355 0.313. 

Закономерности изменения х[к*Си п ] = к(Т, п) коррелируют с результатами изучения 
расплавов индивидуальных щелочных металлов (ЩМ) с учетом, кроме атомов, мета- 
стабильных кластеров ЩМ 2 - ЩМ 5 [12, 13]. 

Представление о содержании газовой фазы над модельными расплавами пока- 
зывает рис. 2. С ростом температуры от 1400 К до 2800 К содержание Си,. 3 в газовой 
фазе закономерно возрастает, М(Си,) > М(Си 2 ) > М(Си 3 ) при каждой температуре, т.е. 
коррелирует с закономерностями изменения х(к*Си м ) в конденсированной фазе меди 

и с данными о М(Іл,. 5 ) = Г(Т, п), полученными в системе расплавов лития + Аг с учетом 
кластеров с п = 2-5 [12, 13]. 
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Рис. 2. Соотношение чисел молей (МО летучих частиц Сиі (1), Си 2 (2) и Си 3 (3) в 
расчете на 1кг исходного рабочего тела в зависимости от температуры 

4. ЗАКЛЮЧЕНИЕ 

С использованием рассчитанных, по-видимому, впервые термодинамических свойств и 
функций метастабильных кластеров к*Си 2 , к*Си 3 , к*Си 6 и Си 3 выполнены равновесные 
компьютерные эксперименты в системе Си + Аг при 400-3000 К и общем давлении 
10~Па. Показана термодинамическая возможность существования в твердом и жидком 
состояниях меди, кроме атомов, конденсированных кластеров меди; рассчитанные со- 
ставы конденсированной и газовой фаз. Высказано предположение о возможности об- 
разования при 400-1357 К, наряду с кристаллами из атомов меди, "кластерных" кри- 
сталлов из димеров, тримеров и двойных тримеров на основе плотных кубических и 
гексагональных упаковок. 

Закономерности изменения состава модельных конденсированной и газовой фаз 
от температуры и числа атомов в кластерах коррелируют с результатами компьютер- 
ных экспериментов по определению состава модельных расплавов и газовой фазы с 
участием атомов и малых кластеров щелочных металлов. Значения Т и АН кипения мо- 
дельного расплава меди согласуются с известными данными. 

Работа выполнена при финансовой поддержке программы Президиума РАН 
"Фундаментатьные проблемы физики и химии наноразмерных систем и наноматериа- 
лов". 
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КОМПЛЕКСООБРАЗОВАНИЕ ХРОМА (III) С МОНО АМИННЫМИ 
КАРБОКСИМЕТИЛЬНЫМИ КОМПЛЕКСОНАМИ И САЛИЦИЛОВОЙ 
КИСЛОТОЙ В ВОДНЫХ РАСТВОРАХ 

В.И. КОРНЕВ, Г. А. МИКРЮКОВА 


Удмуртский государственный университет, Ижевск, Россия 


АННОТАЦИЯ. Спектрофотометрическим методом изучено взаимодействие 
хрома(ІІІ) с метилиминодиуксусной (Н 2 Міс1а), (3-гидроксиэтилиминодиуксусной 
(Н 2 НеісІа), нитрилотриуксусной (Нз№а) и салициловой (Н 2 3а1) кислотами. Показано, 
что комплексообразование в двойных и тройных системах зависит от концентрации 
реагирующих компонентов и кислотности среды. Обнаружены следующие комплек- 
сы: СгНЗа1 2+ , СгЗаГ, СгМібаНЗаІ, СгМЫаЗаГ, СгНеібаНЗаІ, СгНеібаЗаГ, Сг№аНЗаГ 
и СгИіаЗаІ" . Логарифмы констант устойчивости этих комплексов, рассчитанные при 
ц = 0,1 (ЪІаСЮД, соответственно равны: 19,04 ± 0,27, 17,06 ± 0,24, 29,64 ± 0,57, 28,18 
± 0,16. 29,73 ± 0.17, 26,46 ± 0,39, 31,37 ± 0,12, 27,67 ± 0,72. 


Широкое исследование металлов, в том числе и комплексонатов хрома(Ш), определя- 
ется, прежде всего, растущей практической значимостью этого класса соединений. С 
точки зрения координационной химии весьма интересными являются комплексы хро- 
ма(Ш) с комплексонами, производными карбоксиалкилированных аминов. Причем при 
образовании монокомплексонатов хрома(ІІІ) наблюдается неполное насыщение коор- 
динационной сферы металла комплексоном. Поэтому к такому комплексу можно при- 
соединить дополнительный лиганд, который явится элементом достройки структуры и 
максимального использования координационной емкости центрального атома. Иссле- 
дованию комплексообразования таких смешанолигандных комплексов и посвящена 
данная работа. В качестве первичного лиганда выступали моноаминные комплексоны: 
метилиминодиуксусная (Н 2 Міс1а), (3-гидроксиэтилиминодиуксусная (Н 2 НеісІа) и нитри- 
лотриуксусная (Нз№а) кислоты. В качестве дополнительного лиганда была 
использована салициловая кислота (Н 2 3а1). Литературные данные по гетеролигандному 
комплексообразованию в системах типа Сг(ІІІ) — ЦКотр— Н 2 3а1, где ЦКотр — молекула 
комплексона, отсутствуют. Нами не найдено также каких-либо сведений и о комплек- 
сах хрома(ІІІ) с салициловой кислотой. Поэтому нам пришлось изучить и комплексооб- 
разование в системе Сг(Ш)-Н 2 За1. Следует отметить, что характер равновесий в систе- 
ме Сг(ІІІ)-НіКотр был детально описан нами в работе [1]. 
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Исследования проводили спектрофотометрическим методом (прибор СФ-26, 
кварцевые кюветы длиной 5 см). Концентрацию ионов водорода измеряли на иономере 
И-130.2 М (стеклянный электрод ЭСЛ-43-07, хлорсеребряный электрод ЭВЛ-1М3.1). 
Иономер калибровали и проверяли на стенде УПКП-1, а также при помощи стандарт- 
ных буферных растворов, приготовленных из фиксаналов. Необходимую кислотность 
среды поддерживали растворами НС10 4 и ИаОН. Постоянство ионной силы (р = 0,1) 
создавали при помощи ИаСІСД марки “ч.д.а.”. Рабочие растворы комплексонов и сали- 
циловой кислоты готовили растворением точных навесок химически чистых препара- 
тов в воде. Раствор перхлората хрома получали растворением точной навески химиче- 
ски чистого металла в хлорной кислоте. Для увеличения скорости реакций приготов- 
ленные смеси кипятили в течение 15 минут. Математические расчеты проводили по 
программе СРЕ88Р [2]. 

Исследование комплексообразования в двойных и тройных системах было осно- 
вано на изменении оптической плотности гидратированных ионов хрома(ІІІ) в присут- 
ствии различных лигандов. Добавление лигандов к раствору соли металла приводит к 
гиперхромному эффекту и гипсохромному сдвигу. 

Для изучения процессов, протекающих в исследуемых системах, была выбрана 
оптимальная длина волны 540 нм. 

По характеру кривой А = ДрН) (рис. 1, кривая 2) можно заключить, что в систе- 
ме Сг(ІІІ)— НгЗаІ протекает процесс комплексообразования, зависящий от кислотности 
раствора. Известно, что в растворах комплексонатов хрома в кислой среде образуется 
по одному комплексу [1]. 



Рис. 1 . Зависимость оптической плотности растворов Сг(СЮ 4 ) 3 (1) и его смесей с 
Н 2 §а1 (2), Н 2 Міс1а (3), Н 2 Неісіа (4), Н 3 №а (5), Н 2 Міс1а + Н 2 8а1 (6), Н 2 НеЫа + 
Н 2 8а1 (7), Н 3 №а + Н 2 3а1 (8) состава 1:2 (3), 1:3 (2, 4), 1:5 (5) и 1:3:3 (6, 7, 8) 
от рН; С Сг = 8,48 • 1СГ 4 моль/л, А. = 540 нм 
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Мольное соотношение компонентов в двойных системах в области значений рН 
образования и существования комплексов, равное 1:1, устанавливали методами насы- 
щения и изомолярных серий. В качестве примера на рис. 2 приведена диаграмма насы- 
щения в системе Сг(ІІІ)— НіЗаІ (кривая 1). 

Анализ кривых А = ДрН), спектров поглощения реагирующих компонентов и их 
смесей показывает, что в тройных системах Сг(ІІІ)-НіКотр-Н 2 8а1 протекает смешано- 
лигандное комплексообразование. Действительно кривые 6, 7 и 8 (рис. 1), характери- 
зующие тройные смеси при рН > 2,0, заметно отличаются от соответствующих кривых 
3, 4 и 5 для двойных систем. Они не совпадают с ними и не являются их полусуммой. 

Образование гетеролигандных комплексов в системах Сг(Ш)-Н|Котр-Н 2 8а1 
можно подтвердить также кривыми насыщения. Так насыщение смеси Сг(ІІІ)— НгЗаІ 
(1:3) метилиминодиуксусной кислотой (рис. 2, кривая 2) приводит к увеличению опти- 
ческой плотности до некоторой постоянной величины, соответствующей поглощению 
тройной системы Сг(ІІІ)-Н 2 МісІа-Н 2 8а1 при мольном соотношении 1:1:1. С другой сто- 
роны, добавление раствора салициловой кислоты к метилиминодиацетату хрома соста- 
ва (1:3) также указывает на образование гетеролигандного комплекса, поскольку опти- 
ческая плотность уменьшается до постоянной величины (рис. 2, кривая 3). При этом 
мольное отношение компонентов в смешанном комплексе равно 1:1:1. Аналогичные 
равновесия прослеживаются и в других тройных системах. 


Г /Г 
Н 2 8аГ с Сг 



( 3 ) 


Рис. 2. Кривые насыщения в системе Сг 3+ -Н 2 8а1 (1), а также смесей 
Сг 3+ + Н 2 8а1 (1:3) (2) и Сг 3+ + Н 2 Міс1а (1:3) (3) растворами Н 2 МЫа (2) и 
Н 2 8а1 (3) при рН = 3,1 3,3; Сс г = 8,48 ■ ІО -4 моль/л, X = 540 нм 


ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. Том 7, № 1 


73 


В. И. КОРНЕВ, Г. А. МИКРЮКОВА 


Математическое описание равновесий, как в двойных, так и в тройных системах, 
содержащих в растворе несколько гидролизованных ионов хрома(ІІІ), а также различ- 
ные протонированные формы полидентатных лигандов, представляет собой весьма 
сложную задачу. Учитывая, что величины определяемых констант устойчивости ком- 
плексов зависят от того, каким набором комплексных частиц описывается та или иная 

система, то выбор правильного набора комплексных форм имеет принципиальное зна- 
чение. 


В программе СРЕ88Р, которая использовалась в данной работе, построение мо- 
дели равновесной системы сводится к последовательному поиску адекватной модели в 
форме закона действующих масс и на основании физико-химических измерений, в на- 
шем случае величины оптической плотности. В первом приближении формируется со- 
вокупность наиболее вероятных комплексов, а затем производится последовательное 
усложнение модели путем включения в нее только тех комплексов, которые вносят 
ощутимый вклад в среднеквадратичное отклонение и усредняют его. 

Эта программа с наибольшей надежностью позволяет оценить параметры равно- 
весной системы, а также количество реакций, их стехиометрию и термодинамические 
характеристики химических процессов. 

Вопрос о целесообразности выбора той или иной формы диссоциации гидрокси- 
да хрома(ІІІ), комплексона или салициловой кислоты решался на основе минимизации 
значения критерия Фишера. 


Каждая модель содержала обязательные компоненты двойных или тройных сис- 
тем: Н п 8аГ 2 , ЬфКотр'Л продукты гидролиза хрома(ІІІ) Сг(ОЩ н и наиболее вероят- 
ные комплексные ионы состава: Сг р Н п 8а1 ч 3р+п ~ 24 , Сг р НіКошр 2 3р+і_тг и 

Сг р НіКотр 2 Н п За1 ч 3р гт ‘ ч , где т - основность комплексона. Учет гидроксоформ ком- 
плексов не приводит к уменьшению значения критерия Фишера. Наилучшие результа- 
ты получаются при р = 2 = ч = 1, что согласуется с экспериментальными данными. 

Ступенчатые константы диссоциации комплексонов, салициловой кислоты и 
константы устойчивости монокомплексонатов хрома(ІІІ), используемые в данной рабо- 
те, приведены в таблице 1 . 


Таблица 1 . Константы диссоциации лигандов [3 - 9] и константы устойчивости моно 


комплексонатов хрома(Ш) [1] для ц = 0,1 


Лиганд 

рК диссоциации 

Литература 

Комплекс 

1§Р 

рКі 

рК, 

рк 2 

Н 2 8а1 


3,00 

13,60 

[3] 



Н 2 Міс1а 

1,62 

2,12 

9,65 

И, 5] 

СгМіёа" 

11,08 

Н 2 Неі<За 

1,60 

2,20 

8,73 

[6] 

СгНеЫа" 

10,68 

Нз№а 

ТТттст НТТчНо 

1,10 

V , — О 73 

1,84 

2,49 

[7-9] 

Сг№а 

12,37 


74 


ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. Том 7, № 1 


КОМПЛЕКСООБРАЗОВАНИЕ ХРОМА (III) С МОНОАМИННЫМИ КАРБОКСИМЕТИЛЬНЫМИ 
КОМПЛЕКСОНАМИ И САЛИЦИЛОВОЙ КИСЛОТОЙ В ВОДНЫХ РАСТВОРАХ 


Моноаминные комплексоны, применяемые в работе, можно представить сле- 
дующей общей формулой 


+ 


КНК(СН 9 СООН) 9 , 


где К - СНз (Н 2 Мі<іа), СН 2 СН 2 ОН (Н 2 Неіёа) и СН 2 СООН (Н 3 №а). 

Образование салицилатных комплексов хрома(ІІІ) удается описать с помощью 
математической модели только в том случае, если предположить, что связывание Сг 3 ^ в 
растворе при 0,5 < рН < 4,0 происходит преимущественно за счет реакций: 

Сг 3+ + Н 2 8а1 = СгН8а1 2+ + Н + 

СгН8а1 2+ = Сг8аГ + Н + 

Константы устойчивости этих комплексов соответственно равны 19,04 ± 0,27 и 17,06 ± 
0,24. Характер равновесий в системе Сг(ІІІ)— Н 2 8а1 наглядно отображает диаграмма до- 
левого распределения частиц в зависимости от рН. Как следует из рис. 3 (а) максимум 
образования комплекса СгН8а1 2+ достигается при рН 1,2 (78 %). При рН 2 долевое со- 
держание комплексов выравнивается (по 50 %), а при рН 4 наблюдается стопроцентный 
выход комплекса Сг8а1 + . 




Рис. 3. Диаграммы долевого распределения частиц в системе Сг(ІІІ)— Н 2 8а1 (а) и 
Сг(ІІІ)-Н 2 Неіс1а-Н 2 8а1 (б) в зависимости от рН для: Сг(СЮ 4 ) 3 (1), СгН8а1 2+ 
(2), Сг8аГ (3), СгНеЫаН8а1 (4) и СгНеМа8а1 (5); С Сг = 8,48 • Ю^ 4 моль/л. 
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При описании гетеролигандных комплексов мольного состава 1:1:1 наилучшие 
результаты получаются только с учетом протонированной формы салициловой кисло- 
ты. Учет этой структуры значительно понижает значение критерия Фишера. Установ- 
лено, что во всех исследуемых системах образуется по два комплекса (табл. 2). Сравне- 
ние устойчивости соответствующих комплексов показывает, что она изменяется парал- 
лельно с изменением основности аминогруппы комплексона, а также с увеличением 
основности донорного атома кислорода фенольного гидроксила за счет влияния орто- 
расположенного карбоксилат-иона молекулы салициловой кислоты. 

Необходимо учитывать, что кислород недиссоциированной ароматической гид- 
роксогруппы может принимать участие в координации ионом металла. Для характери- 
стики устойчивости протонированных смешанных комплексов удобнее использовать 
не полные константы устойчивости, а частные константы равновесия реакций. Значе- 
ния констант равновесия частных реакций были рассчитаны исходя из констант устой- 
чивости комплексов и соответствующей константы диссоциации салициловой кислоты, 
поэтому доверительные интервалы для 1§ К в табл. 2 не приводятся. 

Таблица 2. Константы равновесия реакций и константы устойчивости гетеролиганд- 
ных комплексов хрома(ІІІ) для р = 0,1 (ЛаСЮД 


Уравнение диссоциации комплекса 

1§К 

1§Р 

СгМіс1аН8а1 = Сг 3+ + Місіа 2 ' + 8а1 2 ' +Н + 

- 

29,64 ±0,57 

СгМіс1аН8а1 = Сг 3+ ± Мі<3а 2 ~ + Н8аГ 

16,04 

- 

СгМЫаЗаГ = Сг 3+ + МіОа 2 ' + 8а1 2 ' 

- 

28,18 ± 0,16 

СгНеИаНЗаІ = Сг 3+ + Неісіа 2 ' + 8а1 2 ' + Н + 

Н 

- 

29,73 ±0,17 

СгНеісІаНЗаІ = Сг 3+ + Неісіа 2 ' + Н8аГ 

16,13 

- 

СгНеИаЗаГ = Сг 3+ + Неісіа 2 ' + Заі 2 ' 

- 

26,46 ±0,39 

СгМаНЗаГ = Сг 3+ + №а 3 ~ + 8а1 2 ' + Н + 

- 

31,37 ± 0,12 

Сг№аН8аГ= Сг 3+ + №а 3 ~ + Н8аГ 

17,77 

- 

Сг№а8а1 2 " = Сг 3+ + №а 3 ~ + 8а1 2 ' 



27,67 ±0,72 
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Как видно из табл. 2 протонированные гетеролигандные комплексы с участием 
НгМісіа и НгНеісІа близки по устойчивости. Комплекс же Сг№аН8а1 более чем на поря- 
док прочнее первых двух. Этот факт можно объяснить тем, что в комплексе Сг№аН8а1 
нитрилотриацетат, возможно, ведет себя как тетрадентатный лиганд, образуя с ионом 
хрома(ІІІ) дополнительный пятичленный цикл за счет группы -СН 2 СООН. Однако, этот 
цикл должен быть стерически неблагоприятным и, вероятно, легко раскрывается в слу- 
чае среднего гетеролигандного комплекса. Действительно, устойчивость комплекса 
Сг№а8аг близка к устойчивости комплекса СгМісІаЗаІ, в котором образование допол- 
нительного цикла невозможно, поскольку НгМіба является трехдентатным лигандом. 
Раскрытие дополнительного цикла в комплексе Сг№а8а1 2 ~, очевидно, связано с пере- 
распределением электронной плотности, которое возникает в результате отщепления 
протона от фенольного гидроксида, связанного с сопряженной системой бензольного 
ядра. 

Структуры смешанолигандных комплексов с участием Н 3 №а можно представить 
следующими скелетными схемами: 





ОоСНіСМ Сг о 



Сг1МіаН8аГ 


Сг№а8а1 2 ~ 


При таком представлении структуры среднего смешанного комплекса у хрома(ІІІ) ос- 
тается еще одно координационное место, которое в водном растворе может быть занято 
молекулой воды. Хотя с другой стороны кажется маловероятным, что молекула воды 
занимает у металла координационное место, а энергетически выгодная метилкарбок- 
сильная группа остается свободной от координации. 

Обращает на себя внимание и тот факт, что устойчивость комплекса СгНеісіа8аГ 
несколько ниже, чем устойчивость других средних комплексов. Очевидно, это связано 
со значительным отрицательным индукционным эффектом оксиалкильной группиров- 
ки в молекуле НгНеісіа. Однако этот эффект мало сказывается в случае протонирован- 
ного комплекса. Возможно, что на характер связей в гетеролигандных комплексах хро- 
ма(ІІІ) влияют еще какие-то процессы, природа которых остается пока невыясненной. 

Для всех гетеролигандных комплексов были построены диаграммы долевого 
распределения частиц в зависимости от рН. В качестве примера на рис. 3 (б) представ- 
лено долевое распределение в системе Сг(ІІІ)-Н 2 Неіёа-Н 28 а 1 . Как следует из рис. 3 (б) 
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при рН 1,6 в растворе присутствуют в основном комплексы СгНЗа1 2+ и СгНеібаНЗаІ, 
доля их накопления составляет по 40 %. Понижение концентрации водородных ионов 
приводит к накоплению комплекса СгНеісІаНЗаІ, долевое содержание которого дости- 
гает 80 % при рН 2,4. Дальнейшее повышение рН растворов способствует снижению 
доли этого комплекса и накоплению нормального комплекса СгНеісІаЗаГ. Так при рН 
4,0 в растворе существуют комплексные частицы СгНеіёаНЗаІ (18 %) и СгНеісІаЗаГ (80 
%). В системе Сг(ІІІ)-Н 2 Міс1а-Н28а1 при рН 2,5 накапливается нормальный комплекс 
СгМЫаЗаІ (80 %), а в системе Сг(Ш)-Н 3 №а-Н28а1, при этом же значении рН, преобла- 
дает протонированный комплекс состава Сг№аН8аГ (85 %). При рН > 4,0 в обеих сис- 
темах накапливаются только нормальные гетеролигандные комплексы. 

В заключение следует отметить, что гетеролигандные комплексы хрома(ІІІ) на 
основе моноаминных карбоксилированных комплексонов и салициловой кислоты име- 
ют октаэдрическое строение и образуются они за счет вытеснения воды из гексааква- 
комплекса хрома. Количество молекул воды во внутренней сфере гомо- и гетероли- 
гандных комплексов можно установить исходя из координационного числа хрома(ІІІ), 
равного шести и дентатности лигандов, входящих в состав того или иного комплекса. 
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ЗЦММАКѴ. Іпіегасбоп ЬеСѵѵееп сЬгошіит(ІИ) апё теіЬШтіпобіасеІіс (МША 
Н 2 Мі<За), Р-ЬуёшхіеіЬуІітіпобіасебс (НЕГОА, Н 2 Неіба), піігііоігіасейс (ХТА, Н 3 №а) 
апб аізо заіісіі (Заі, Н 2 3а1) асЫз \ѵаз зШбіеб Ьу зресСгоіЫотеігіс шеазигешет. И ѵѵаз 
Ьееп зЬоѵѵп Шаі сотріех І'огшаііоп іп биріех апб Ігіріех зузіет із берепсі оп Йіе геа§ет 
сопсепігаіюп апё рН. И \ѵаз Гоипб Нте пехі сотріехез: СгНЗа1 2+ , СгЗаГ, СгМЫаНЗаІ 
СгММаЗаГ, СгНеібаНЗаІ, СгНеісІаЗаГ, СгМаНЗаГ и Сг№аЗа1 2 -. ТЬе зСаЬіПіу сопзСат 
1о§агііЬт оГ Фезе сотріехез аі ц = 0,1 (ИаСЮД аге 19,04 ± 0,27, 17,06 ± 0,24. 29.64 ± 
0,57, 28,18 ± 0,16, 29,73 ± 0,17, 26,46 ± 0,39, 31,37 ± 0,12, 27,67 ± 0,72 гезресбуДу.’ 


78 


ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. Том 7, № 1 


УДК 541.11. 


СТАНДАРТНЫЕ ЭНТАЛЬПИИ ОБРАЗОВАНИЯ (СЭО) 
КОНДЕНСИРОВАННЫХ МЕТАСТАБИЛЬНЫХ БИНАРНЫХ 
САМОАССОЦИАТОВ (КЛАСТЕРОВ) РАЗЛИЧНЫХ ЭЛЕМЕНТОВ 

Г.К. МОИСЕЕВ 

Российская академия наук, Уральское отделение, ГУ Институт металлур- 
гии, Россия, Екатеринбург 

АННОТАЦИЯ. Для 25 бинарных конденсированных кластеров элементов ІА, ПА и 
ІНВ-ѴІВ подгрупп Системы исследованы зависимости известных величин ДН° 98 от 
размеров атомов, положительных ионов и величин первых потенциалов ионизации 
для элементов каждой подгруппы. Предложены эмпирические уравнения для оцен- 
ки величин ДН° 98 . 


1. ВВЕДЕНИЕ 

В работах [1-3] предложены различные методы оценки величин ДН° 98 (стандартных 
энтальпий образования - СЭО) метастабильных конденсированных самоассоциатов 
("малых" кластеров, в том числе бинарных) для отдельных элементов (Эл) системы: 

СЭО(к*Эл п )/СЭО(Эл п ) ~ СЭО(к*С п )/СЭО(С п ), ( 1 ) 

СЭО(к*Эл п ) ~ 5940. 1 пЕ в (к*Эл п ), Дж/моль, (2) 

СЭО(к*Эл п ) ~ 5940. 1 п-Е в (к*С п )ОЭО(Эл)/ОЭО(С), Дж/моль, (3) 

іде к*Эл п и Эл п — соответственно кластеры в конденсированном и летучем состояниях; 
СЭО(к*С п ) и СЭО(Сп) - АН° 98 кластеров углерода по данным [4, 5]; Е в (к*Эл п ) - энер- 
і ия в расчете на 1 атом (эВ/атом) системы из п атомов или средняя энергия связи между 
атомами в изолированном п- атомном кластере [6]; ОЭО (Эл) и ОЭО (С) - относитель- 
ные электроотрицательности Эл и углерода по [7]. Применение (1) — (3) приводит к со- 
гласованным величинам СЭО (к*Эл п ) для элементов І-ІѴ групп Системы [2]. Исходные 
данные Е в (к*Эл п ) взяты по данным молекулярной динамики (например, [8, 9]), 
СЭО(Эл п ) - из БД АСТРА. ВА8 [10], СЭО(к*С п ) и СЭО(С п ) - из [4, 5]. 
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К настоящему времени оценены СЭО для нескольких десятков конденсирован- 
ных малых кластеров с числом атомов от 2 до 8 [1, 3], однако наиболее представи- 
тельными являются сведения о бинарных кластерах (см. табл. 1), что дает возможность 
проанализировать зависимости СЭО кластеров в каждой из подгрупп Системы от раз- 
личных физико-химических характеристик атомов, составляющих самоассоциаты. Ли- 
тературных данных о систематических исследованиях подобного рода не обнаружено. 

Целями настоящей работы являлись: - выявление закономерностей изменения 
СЭО бинарных кластеров из элементов ІА, ПА и ІІІВ-ѴІВ подгрупп Системы в зависи- 
мости от размеров атомов, размеров положительных ионов, величин первого потенциа- 
ла ионизации (ППИ); - конструирование и проверка эмпирических уравнений для оцен- 
ки величин СЭО бинарных кластеров. 

2. МЕТОДИКА ИССЛЕДОВАНИЙ 

Для определения зависимостей СЭО(к*Эл 2 ) - Г(г ат ., г+ иона , ППИ) использованы приве- 
денные значения ЛН° 98 ь равные 

ДН2 98 (к*Эл2)і/Ыі, кДж/(моль-Мі), (4) 

где Ыі - номер і - того элемента в Системе, эквивалентный величине положительного 
заряда ядра атома і - того элемента Ъ ѵ Предварительный анализ показал, что рацио- 
нальнее использовать величину І§ ^ для определения численных зависимостей от 
Га ‘ (і)’ г +ион(і) и ППИ(і), т.е. при этом выявляются более простые математические урав- 
нения. Исходные данные приведены в табл. 1 . Сведения о СЭО кластеров взяты из [1-3] 
и для ряда элементов (к*Ва 2 , к*Р 2 , к*8 2 ) рассчитаны по (1) с использованием данных 

[10] . Характеристики атомов различных элементов (г ат , г +ионов , ППИ) взяты по данным 

[ 11 ] . 

3. РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ 

На рис. 1-3 показаны результаты анализа данных табл. 1, т.е. численные значения Т, 
= Дг ат ., г+иона, ППИ)і для 25 изучаемых бинарных кластеров элементов каждой из под- 
групп системы. Из рисунков можно сделать следующие общие заключения: 

- во всех подгруппах наблюдаются зависимости 1§ Т,- от радиусов атомов, ионов и пер- 
вых потенциалов ионизации (ППИ), причем с увеличением радиусов атомов и положи- 
тельных ионов У = 1§ *Р уменьшается; с увеличением ППИ - увеличивается; с ростом 

номера элемента (или заряда ядра атома) наблюдается уменьшение величин У в под- 
группе; 
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Таблица 1. Исходные данные для определения численных закономерностей изменения 
СЭО бинарных кластеров от радиусов атома, иона и первого потенциала 
ионизации 


N элемен- 
та 

Конденс. 

бинарный 

кластер 

ДН“ 98 , 

кДж/моль 

1§Т 

Гатома.э^ 

Г+ІІОНз А 

1 потенц. 

ионизации 

атома, еВ 

Подгруппы ІА и ІВ 

1 "> 

2) 

ІЛ 2 

9.35 

0.4937 

1.55 

0.68 

5.390 

11 

N32 

6.38 

-0.2365 

1.89 

0.98 

5.138 

се: 

К 2 

5.76 

-0.5186 

2.36 

1.33 

4.339 

37 

НЬ 2 

5.08 

-0.8623 

2.48 

1.49 

4.176 

55 

Сз 2 

4.90 

-1.0500 

2.67 

1.65 

3.893 

29 

Си 2 

24.75 

-0.069 

1.28 

0.93 

7.726 

Подгруппа НА 

I 4 

Ве 2 

27.24 

0.8331 

1.13 

0.34 

9.320 

П ІТ~ 

С— 

М§ 2 

15.10 

0.0998 

1.60 

0.74 

7.640 

ро^ 

Са 2 

15.78 

-0.1029 

1.97 

1.04 

6.111 

1 38 

Зг 2 

14.30 

-0.4244 

2.15 

1.20 

5.692 

: 56 

Ва 2 

16.05 

-0.5435 

2.21 

1.38 

5.210 

Г 

одгруппа III 

В 


|5 

В 2 

38.55 

0.8870 

0.91 

0.23 

8.298 

13 

АІ2 

21.58 

0.2200 

1.43 

0.57 

5.986 

Г 

Іодгруппа ІѴ1 

3 


I 6 

1 

с 2 

38.30 

0.805 

0.77 

0.15 

11.261 

! 14 

ЗІ2 

26.00 

0.269 

1.18 

0.41 

8.149 

32 

Ое 2 

22.30 

-0.1096 

1.39 

0.53 

7.880 

50 

Зп 2 * 

19.10 

-0.3587 

1.58 

0.71 • 

7.322 

82 

РЬ 2 

14.10 

-0.6421 

1.75 

0.84 

7.415 


Подгруппа ѴВ 



тт - 

Р, 

5.67 

-0.4225 

1.10 

0.35 

-10.5 

33 

Аз 2 

8.80 

-0.5735 

1.21 

0.47 

9.810 

51 

зь 2 

10.875 

-0.6710 

1.41 

0.62 

8.640 

83 

Ві 2 

10.14 

-0.9130 

1.46 

0.74 

7.237 


п 

одгруппа ѴІЕ 

1 


16 

§2 

5.93 

-0.4312 

1.04 

1.84 

10.357 

34 

8е 2 

6.65 

-0.7087 

1.16 

1.98 

9.750 

52 

Те 2 

7.52 

-0.8398 

1.35 

2.21 

9.010 
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Рис. 1. Зависимости 1§ У = Г(г ат .) для кластеров элементов ІА, НА, ІІІВ-ѴІВ под- 
групп системы (1-6). Названия элементов приведены на рисунке 
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Рис. 2. Зависимости 1§ Т = Г(г +ион .) для кластеров элементов ІА, НА, ІІІВ-ѴІВ под- 
групп системы (1-6). Названия элементов приведены на рисунке 
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Рис. 3. Зависимости 1§ 'Р - Г(ППИ ) для кластеров элементов ІА, ПА, ІІІВ-ѴІВ под- 
групп системы (1-6). Названия элементов приведены на рисунке 

- для кластеров элементов рассмотренных подгрупп, за исключением ІѴВ подгруппы 
при Хі = ППИі, проявляются отчетливые тенденции к образованию простых линейных 
зависимостей У = Г(х). Не исключено, что ППИ для элементов ІѴВ подгруппы рассчи- 
таны недостаточно корректно (в [11] ППИ этих элементов приняты по данным, опуб- 
ликованным в 1974 г. [12]) и требуют уточнения; 

- 1§ Т(к*Си 2 ), т.е. элемента ІВ подгруппы, в зависимости от аргумента "согласуется" 
или может "включен" в число величин У для кластеров ІѴВ подгруппы (х = г ат и х = 
ППИ) или ІА подгруппы (х - г иона ). Поскольку для других элементов ІВ подгруппы Ае 
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и Аи величины У неизвестны, то вопрос о "принадлежности" 1§ 'Р этого кластера к ІА 
или ІѴВ подгруппам остается открытым; 

- для кластеров ШВ подгруппы известны только две величины 1§ Ѵ Р (для В 2 и Аіг). По- 
этому можно только указывать тенденции их изменения в зависимости от разных аргу- 
ментов. 

Для подгрупп ІА, ПА, ІѴВ-ѴІВ, несмотря на ограниченное число точек (5, 5, 5, 4 
и 3, соответственно), рассчитаны уравнения У = Г(х), среднеквадратичные отклонения 
(8) и коэффициенты корреляции (г), а также средние отклонения расчетных СЭО от 
исходных (см. табл. 2). 

Из табл. 2 видно, что наблюдаются удовлетворительные коэффициенты корре- 
ляции для всех линейных уравнений, за исключением полученных при х = ППИ для ІА 
и ІѴВ подгрупп. Расчетные значения СЭО отличаются от исходных в пределах ±7 - 
±22%. Совокупность полученных данных дает основания полагать, что выбранный 
путь (или идеология) построения численных уравнений с учетом N=2, г ат ., г +ИОН ов и ППИ 
является правомерным и может быть использован. Целесообразно в первую очередь 
для кластеров каждой группы использовать линейные зависимости У = Дг ат и г +иона ), 
поскольку, как уже отмечалось, величины первых потенциалов ионизации могут быть 
недостаточно корректными. 

К сожалению, количество известных СЭО бинарных кластеров в подгруппе ѴІВ 
не дает возможность считать зависимости У = Г(х) в этой подгруппе статистически дос- 
товерными. Тем не менее, мы сочли необходимым оценить СЭО (к*Ро 2 ) в этой под- 
группе, а также СЭО (к*Ргг) и СЭО (к*Яаг) (ІА и ПА подгруппы). Для этих кластеров 
ДН“ 98 (в кДж/моль) соответственно равны 4.98±1.2; 4.81±0.66 и 17.9±1.6. 

4. ЗАКЛЮЧЕНИЕ 

Для 25 бинарных метастабильных конденсированных кластеров ІА, НА, ІѴВ-ѴІВ под- 
групп Системы впервые выполнен анализ зависимостей 1§(СЭО/ІЧ) = Т (г ат , г +ион , ППИ) в 
каждой из подгрупп. Выявлены численные линейные уравнения У = Г (г ат ., г +ио „, ППИ), 
расчеты по которым во всех подгруппах описывают известные величины СЭО класте- 
ров в пределах ±(7-55.5)%, в то время как уравнения У = р (г ат ., г +ИО н) в пределах 
±(7-5-23)%. При анализе зависимостей У = Г (ППИ) в ІА и ІѴВ подгруппах расхождение 
расчетных и известных величин СЭО составляет ±(35—55.5)%. Высказано мнение, что 
для ряда элементов этих подгрупп величины ППИ нуждаются в уточнении. 

С использованием полученных уравнений впервые оценены СЭО конденсиро- 
ванных кластеров Рг 2 , Наг и Рог. 

Представляется, что предложенная методология анализа может быть использо- 
вана при систематическом исследовании других малых кластеров с числом атомов п > 
2 . 
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Таблица 2. Численные уравнения зависимостей 1 8 У = Г(г ах , г +иона , ППИ), рассчитанные 
по данным табл. 1 ; з и г - среднеквадратичные отклонения и коэффициенты 
корреляции 


X* 

СЭО 

У = 1§У=1§ ** 

N 

3 

г 

Отклонения ве- 
личин СЭО рас- 
четных от СЭО 

исходных, % 

ІА 

1 

Гат. 

2.19784-1.21389-х 

0.1182 

0.9816 

±22.9 

Г+иона 

1.29400-1.40604-х 

0.1169 

0.9800 

±18.95 . 

ППИ 

-4.5 1 862+0. 89028-х 

0.2263 

0.9463 

±35.4 


ПА 


Гат. 

2.19784-1.21389-х 

0.1182 

0.9816 

±15.5 

Г+иона 

1.29400-1.40604-х 

0.1169 

0.9800 

±18.95 

ППИ 

-2. 18750+0. 31788-х 

0.1235 

0.9805 

±18.56 


ІѴВ 



Гат. 

1.96872-1.48126-х 

0.0353 

0.9985 

±16.4 

Г+иона 

1.09889-2.09502-х 

0.0636 

0.9952 

±17 

ППИ 

-2.62092+0.3 1095-х 

0.2847 

0.8994 

±55.5 

ППИ 

-17.31317+3. 55559-х- 

0.2064 

0.9659 

±36.3 


0.17285-х 2 





ѴВ 

— 

Гат. 

0.81584-1.12806-х 

0.095 

0.9262 

±20.1 

Г+иона 

-0.0022371-1.17938-х 

0.0533 

0.9774 

±7.1 

ППИ 

-1.93306+0. 1423 8-х 

0.0371 

0.9891 

±6.3 


ѴІВ 

] 

Гат. 

0.83238-1.26108-х 

0.0966 

0.9449 

±12.3 

Г+иона 

1.45362-1.05150-х 

0.0994 

0.9416 

±12.6 

ППИ 

-3.55988+0. 29879-х 

0.0761 

0.9661 

±9.7 


■ г ахома и г+ иона в А; ППИ — в еВ . 


- СЭО ДН 298 в кДж/моль; N — номер элемента в Системе. 


Работа выполнена при финансовой поддержке программы исследований Прези- 
диума РАН "Фундаментальные проблемы физики и химии наноразмерных систем и на- 
номатериалов"; РФФИ (код проекта N04-03-33109). 
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ИССЛЕДОВАНИЕ ВНУТРЕННЕЙ ДИНАМИКИ ТРЕХАТОМНЫХ 
КЛАСТЕРОВ БЛАГОРОДНЫХ ГАЗОВ МЕТОДОМ ВЫДЕЛЕНИЯ 
ЭФФЕКТИВНЫХ МОД ДВИЖЕНИЯ 

А. А. РЫБАКОВ 

Московский государственный университет имени М.В. Ломоносова, 
Москва, Россия 

АННОТАЦИЯ. В работе предлагается метод определения мод движения и эффек- 

тивной размерности фазового пространства в задаче анализа внутренней динамики 
нелинейной многочастичной системы. Метод протестирован на модели гармониче- 
ского осциллятора. Представлены результаты исследования внутренней динамики 
трехатомных кластеров благородных газов с помощью данного метода. 


ВВЕДЕНИЕ 

Малые кластеры благородных газов представляют собой слабосвязанные системы, 
взаимодействие между атомами, в которых носит Ван-дер-ваальсовый характер и хо- 
рошо описывается потенциалом Леннарда-Джонса [1]. Взаимодействие атомов в кла- 
стере характеризуется сильной нелинейностью, вызванной большой величиной оттал- 
кивания на малых расстояниях. В силу такой ангармоничности потенциала взаимодей- 
ствия, даже в тримерах благородных газов возможно хаотическое движение, наблю- 
даемое даже при энергиях, незначительно превышающих значение энергии нулевых 

колебаний, при этом объем хаотической компоненты возрастает с ростом полной энер- 
гии системы. 

В рамках классического подхода к описанию внутренней динамики малых кла- 
стеров используются следующие методы численного эксперимента: метод Монте- 
Карло и метод молекулярной динамики. Так, интересным явлением, обнаруженным с 
помощью этих методов, является квазиплавление и квазиотвердение в кластерах. Пока- 
зано, что квазитвердая фаза характеризуется регулярным движением составляющих 
кластер частиц около положений равновесия, а квазижидкая — хаотическим движени- 
ем. 

Традиционно во многих работах изучается структура потенциальной энергии 
кластеров с разной степенью агрегации. Для многочастичных кластеров, состоящих из 
достаточно большого числа частиц, поверхность потенциальной энергии принадлежит 
многомерному пространству, и ее наглядное исследование затруднительно. Однако на- 
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хождение минимумов и седловых точек на этой поверхности представляет собой впол- 
не осуществимую и уже решенную для многих конфигураций задачу даже для класте- 
ров с достаточно большим числом частиц. Возможно также нахождение путей изоме- 
ризации кластера как энергетически наиболее выгодных переходов от одного миниму- 
ма к другому через связывающее их седло [2]. 

Трехатомные кластеры атомов благородных газов имеют следующую поверх- 
ность потенциальной энергии: существуют глобальный минимум, соответствующий 
конфигурации равностороннего треугольника, и три седловых точки, соответствующие 
линейным конфигурациям, различающимся перестановкой частиц [3]. Особенностью 
трехатомных кластеров является то, что переход между различными треугольными 
конфигурациями может быть осуществлен как с переходом через линейную конфигу- 
рацию, так и простым поворотом без прохождения линейной конфигурации. 

Большинство теоретических работ по исследованию внутренней динамики кла- 
стеров выполнено для невращающихся систем. Однако, для моделей, описывающих 
внутреннюю динамику кластеров благородных газов, наряду с полной энергией систе- 
мы, угловой момент также является интегралом движения. Вращение кластеров явля- 
ется нежестким, и оказывает существенное влияние на их внутреннюю динамику. Так, 
движение кластеров с нулевым угловым моментом стохастично при любых начальных 
условиях и, видимо, эргодично - существует только одна стохастическая компонента, 
заполняющая всю энергетическую поверхность. При учете даже небольшого вращения 
ситуация меняется - в фазовом пространстве появляются области регулярного движе- 
ния. [4,5] Установлено, что зависимость доли регулярных траекторий от значения угло- 
вого момента носит немонотонный характер. Влиянию направления суммарного угло- 
вого момента на характер движения атомов в кластере посвящена работа [3]. Выявлены 
качественные и количественные различия в зависимостях максимального показателя 
Ляпунова от направлений вектора углового момента. 

1. ПОНЯТИЕ О МЕТОДЕ ВЫДЕЛЕНИЯ МОД И ОПРЕДЕЛЕНИЯ 
ЭФФЕКТИВНОЙ РАЗМЕРНОСТИ 

Для описания внутренней динамики трехатомного кластера с наложенным условием 
отсутствия поступательного движения необходимо двенадцать координат (шесть коор- 
динат и шесть импульсов). При такой размерности фазового пространства достаточно 
трудно выявить особенности внутренней динамики системы. Для решения этой про- 
блемы в данной работе предлагается использовать разложение движения системы в п- 
мерном фазовом пространстве на ортогональные составляющие, которые мы будем на- 
зывать модами движения, обладающие экстремальными свойствами: суперпозиция из 
т мод наиболее точно аппроксимирует заданные особенности эволюции системы среди 
всех т-мерных составляющих фазовой траектории, т=\,2,...,п. Для многих встречаю- 
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щихся на практике динамических систем их фазовая траектория в многомерном фазо- 
вом пространстве располагается в пространстве меньшего числа координат. Предла- 
гаемая методика позволяет построение такой системы координат, упорядоченной по 
вкладу в эволюцию системы, что при отбрасывании п-т координат мы всегда получаем 
наиболее точную аппроксимацию исходных координат. 

Эволюция изучаемых в работе систем будет рассматриваться в дискретные мо- 
менты времени е [О, Г), N > п ; в этом случае фазовая траектория задается в виде 

матрицы <5 ѵ - - 1,...,ц,у = 1,...,УѴ размера пх N , п - число степеней свободы, ко- 
торую будем называть траекторной матрицей динамической системы и обозначать я 
Решением поставленной выше задачи будет система координат, получаемая в результа- 
те сингулярного разложения траекторной матрицы [ 6 ]. По теореме о сингулярном раз- 
ложении, любая матрица вида может быть записана в виде 


п 

= I X К ММ/’ і = 


к=] 


(1) 


где числа Л 1 ,...,Л Ы , \ > Л 1 > ... > Л п > 0 называются сингулярными числами матрицы 
бч, а ортонормированные наборы векторов е* и к = \,...,п называются сингуляр- 
ными базисами. При этом вектор е* задает “направление” к -ой моды движения, точ- 
нее, направления и относительную амплитуду движения составляющих кластер частиц, 
а произведение Л к {§ к } / - амплитуду, зависящую от времени і . , / = 1,...,7Ѵ^ , то есть пе- 
редает динамику к -ой моды. 

Величина 



А=1 


( 2 ) 


задает вклад А: -ой моды в аппроксимацию матрицы 04 на интервале [О, Г). Если сум- 
ма вкладов первых т мод достаточно велика, то мы можем пренебречь вкладом ос- 
тальных мод и, таким образом, выделить подпространство размерности т<п, наилуч- 
шим образом аппроксимирующим исходное пространство. При этом зависимость раз- 
мерности подпространства т от точности приближения мы будем называть эффектив- 
ной размерностью фазовой траектории. 

Так как разложение справедливо не только для исходной матрицы д , но и для 
любой матрицы вида 0 $ , мы можем раздельно рассматривать преобразования коорди- 
натного и импульсного пространства, задавая оператор 0 , как оператор, оставляющий 
от траекторной матрицы только координационную или только импульсную часть. 
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Рассматривая координатное пространство, отметим, что средний по времени 
квадрат нормы вектора Од характеризует средний размер области пространства, где 
происходит движение материальной точки. В наибольшей степени эта область вытяну- 
та вдоль вектора е, , соответствующего первой моде в нашем разложении, проекция на 
который вектора Ос ] в среднем по времени имеет наибольший квадрат нормы. Анало- 
гично, двумерная плоскость е х е 2 , проекция Од , на которую в среднем по времени име- 
ет максимальный квадрат нормы, описывает ориентацию этой области в двух направ- 
лениях. Проекция “пространственной” части фазовой траектории, на эту плоскость дает 
двумерное описание координатного пространства. 

Если мы будем рассматривать часть траекторной матрицы, содержащую им- 
пульсы составляющих кластер частиц, мы сможем выделить направления в фазовом 
пространстве, в которых сосредоточена большая часть кинетической энергии системы. 

9 

Здесь важно отметить то, что сумма 2] ^ квадратов сингулярных чисел равна средне- 

ш 

му по времени значению кинетической энергии системы, а квадраты соответствующих 
сингулярных чисел Я, 2 - долям средней кинетической энергии в каждой из мод. Проек- 
ция фазового пространства на две первые моды с максимальной точностью дает дву- 
мерное описание пространства импульсов. 

2. АПРОБАЦИЯ МЕТОДА ВЫДЕЛЕНИЯ МОД И ОПРЕДЕЛЕНИЯ 
ЭФФЕКТИВНОЙ РАЗМЕРНОСТИ 


В качестве модельной задачи для апробации метода сингулярного разложения траек- 
торных матриц рассмотрим модель трехмерного гармонического осциллятора. 

Уравнение движения материальной точки во внешнем центральном поле запи- 
сывается в следующем виде [7]: 


г 


дЦ г_ 
дг г 


( 3 ) 


Простейшее (гармоническое) приближение для О (г) будет следующим: 


II (г) = г - (4) 

Отсюда получаем, что, положив частоты колебаний в направлениях всех трех 
осей равными единице, движение описывается следующей системой уравнений: 
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X = -X 

' У = -У (5) 

Было произведено численное интегрирование системы уравнений (5) следую- 
щих, произвольно выбранных начальных условиях: х 0 = 0.25, у 0 = 0.5, = 0.75: 

По = - 1 . Но = 0, ѵ. 0 - -1 . 

После выделения мод в импульсном пространстве были получены следующие 
результаты: 

1 . Подтверждено, что полусумма квадратов сингулярных чисел равна средней по вре- 
мени кинетической энергии системы (квадраты сингулярных чисел имеют значения 
1.302 ПО 3 , 195.254, -9.971 • ІО" 14 соответственно). 

2. Экстремумы кинетической энергии системы совпадают с экстремумами квадратов 
преобразованных координат (рис. 1). Экстремумы зависимостей исходных коорди- 
нат от времени не совпадают с экстремумами зависимости кинетической энергии от 
времени (рис. 1). 



Рис.1. Кинетическая энергия, импульс и проекции импульсного подпространства 
на первые две моды. 

® — кинетическая энергия, проекция на первую моду, — проекция на 
вторую моду, - ■ - - проекция импульса на одну из осей 
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3. В рассматриваемой системе вся кинетическая энергия полностью сконцентрирована 
в двух направлениях. Вклад координаты, соответствующей третьему сингулярному 
числу, практически равен нулю. Это хорошо согласуется с тем, что траектория дви- 
жения частицы в центральном поле лежит целиком в плоскости, перпендикулярной 
вектору углового момента. 

4. Две первые оси нового базиса образуют плоскость, полностью совпадающую с 
плоскостью движения частицы. Преобразованные координаты совпадают с импуль- 
сами в системе координат, в которой плоскость V — V совпадает с плоскостью, в 
которой лежат импульсы частицы. 

3. ПРИМЕНЕНИЕ МЕТОДА ВЫДЕЛЕНИЯ МОД И ОПРЕДЕЛЕНИЯ 

ЭФФЕКТИВНОЙ РАЗМЕРНОСТИ ДЛЯ ИССЛЕДОВАНИЯ ДИНАМИКИ 
ТРЕХАТОМНОГО КЛАСТЕРА 


Кластер из трех одинаковых атомов благородного газа рассматривается как связанное 
или долгоживущее квазисвязанное состояние системы с гамильтонианом 


Н=і 

/= 1 


— • 1 

_д_ 

2т 


+ ХЩг/Д 


( 6 ) 


где Ги = 


я -я, 


Р, и я, 


векторы импульса и координаты і -го атома, т - его мас- 


са. 

Принимая во внимание только попарное взаимодействие между атомами в класте- 
ре, для численного описания выбирается модельный потенциал, хорошо описывающий 
ван-дер-ваальсовые взаимодействия — потенциал Деннарда- Джонса 


ад = 4-с/ 0 . 


г V 2 

а 


г 

V» У 


а 


к г >и 


(7) 


где І1 0 - глубина потенциальной ямы и а - характерная ширина, соответствующая по- 
ложению нуля потенциала. Значения глубины потенциальной ямы и ее характерной 
ширины для различных атомов хорошо известны. Далее мы будем пользоваться естест- 
венной системой единиц, в которой масштабами служат т , 1/ 0 и а. Для кластеров, со- 
стоящих их молекул аргона, т = 39.945 а. е.м. ІІ 0 = 99.55см~' и а = 3.757/2 |/б А . Вы- 
раженная в естественных единицах постоянная Планка Ті = ~ г — . достаточно мала: 

^тіі^а 2 
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Ь 2.93-ю 2 , 1.55 -10 2 и 9.9 -КГ 3 для аргона, криптона и ксенона соответст- 
венно. Малая величина П может оправдывать классический подход к рассмотрению 
динамики кластеров [4]. 

Основное состояние кластера с энергией Е = -3 и нулевым угловым моментом 
соответствует неподвижным атомам, размещенным в вершинах равностороннего тре- 

угольника с длиной стороны а’ =аіІ2. Движение системы в состоянии с максималь- 
ным для данной энергии угловым моментом представляет собой вращение вокруг 

центра симметрии равностороннего треугольника, лежащего в плоскости, перпендику- 
лярной вектору углового момента. 

Структура фазового пространства кластера исследовалась методом молекуляр- 
ной динамики. Схема численного интегрирования уравнений движения модели деталь- 
но описана в [4]- в Р е мя расчета траектории составляло 2500 естественных единиц вре- 
мени. Начальные условия при заданных полной энергии и суммарном угловом моменте 
выбирались случайным образом. Отобранные для исследования траектории относятся к 
областям хаотического и регулярного движения кластера. При этом для определения 
типа движения (регулярного или хаотического) для каждой траектории вычислялся 

максимальный показатель Ляпунова по схеме Бенеттина - Галгани - Стрелсина опи- 
санной в [8]. 

Рассмотрим результаты, полученные для трехатомных кластеров с помощью из- 
ложенной выше методики. 

Сингулярное разложение пространственной части траекторной матрицы позво- 
ляет отметить следующее: 

а) Для невращающегося кластера эффективная размерность регулярной компо- 
ненты значительно меньше, чем хаотической. Регулярное движение атомов кластера 
может быть аппроксимировано 2-мя модами с погрешностью в 2-3 % (рис.2). 

б) Для кластера с высоким угловым моментом (М = 0.8 М тт ) эффективные раз- 
мерности регулярной и хаотической компонент практически совпадают и движение 
системы можно аппроксимировать двумя модами с погрешностью в 5-7 % (рис 2) 

Сингулярное разложение импульсной части траекторной матрицы позволяет от- 
метить следующее: 

а) в случае нулевого углового момента и регулярного движения вся кинетиче- 
ская энергия с точностью 5-7% сконцентрирована в трех направлениях. Хаотическое 
движение системы характеризуется большим числом мод, необходимых для описания 
кинетической энергии системы с той же точностью (рис.З). 

б) при большом угловом моменте (М = 0.8М тах ) кинетическая энергия атомов в 

кластере сконцентрирована, в основном, в двух модах как для регулярной, так и для 

хаотической компонент. Погрешность такого двумерного приближения составляет 12- 
15% (рис.З). 
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Рис. 2. Эффективная размерность координатного подпространства в зависимости 
от относительной погрешности аппроксимации. Полная энергия -2.5. а- 
невращающийся кластер, б - относительный угловой момент М=0.8М тах . ■ 
- регулярное движение; о - хаотическое движение 



Рис. 3. Распределения кинетической энергии по модам импульсного подпро- 
странства. Полная энергия -2.5. а - невращающийся кластер, б - относи- 
тельный угловой момент М=0.8М тах . и - регулярное движение; о - хаоти- 
ческое движение 
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в) для невращающегося кластера двумерные проекции траекторий на направле- 
ния, характеризуемые наибольшей кинетической энергией атомов в кластере, значи- 
тельно различаются для регулярного и хаотического движений: при регулярном движе- 
нии существуют ограничения на взаимные направления импульсов атомов кластера 
(рис.4). 

г) для кластера с большим значением углового момента кластера двумерные 
проекции траекторий на направления, характеризуемые наибольшей кинетической 
энергией атомов в кластере, практически совпадают, и имеют вид окружностей (рис. 5). 



Рис.4. Проекция импульсного подпространства на плоскость еі-е 2 . Полная энер- 
гия —2.5, невращающийся кластер, а — регулярная траектория, б — хаоти- 
ческая траектория 



Рис. 5. Проекция импульсного подпространства на плоскость еі-е 2 . Полная энер- 
гия -2.5, относительный угловой момент 0.8М тах . а - регулярная траекто- 
рия, б - хаотическая траектория 
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Известно, что, помимо ангармоничности потенциала взаимодействия, еще одним 
источником хаоса является кинетическая энергия, сосредоточенная в колебательных 
степенях свободы. Так как кинетическая энергия состоит из колебательной и враща- 
тельной, то представляет интерес выявление различий поведения регулярной и хаоти- 
ческой компонент с ростом углового момента системы. Общепринятой точкой зрения 
здесь является то, что наличие вращения регуляризует динамику [9]. Однако, нами бы- 
ло установлено, что существует немонотонный переход от движения с нулевым угло- 
вым моментом к движению с большим угловым моментом. Это можно объяснить тем, 
что с ростом углового момента сначала число мод увеличивается, так как увеличивает- 
ся и средняя кинетическая энергия, а затем эффективное число мод стремится к двум, 
так как наибольшая часть кинетической энергии сосредоточена во вращении. Такой 
вывод основывается на зависимостях погрешности в кинетической энергии при ап- 
проксимации эволюции системы двумя модами от величины углового момента (рис. 6). 



Угловой момент 

Рис. 6. Зависимость ошибки аппроксимации эволюции системы двумя модами в 
зависимости от величины относительного углового момента. Полная 
энергия кластера -2.5. ■ - регулярное движение; о - хаотическое движе- 
ние 


ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. Том 7, № 1 


97 


А. А. РЫБАКОВ 


ВЫВОДЫ 

Как показали исследования внутренней динамики трехатомного кластера аргона, выде- 
ление главных мод движения позволяет дать детальное описание поведения достаточно 
сложных нелинейных динамических систем. Было показано, что нелинейная динамика 
слабосвязанной системы может быть с хорошей точностью описана небольшим числом 
мод. С помощью методики выделения мод движения обнаружены некоторые новые 
особенности внутренней динамики кластера, приближающие нас к пониманию приро- 
ды хаоса и влиянию нежесткого вращения на внутреннюю динамику. 
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АННОТАЦИЯ. Рассчитаны термодинамические свойства и функции метастабиль- 
ного конденсированного кластера Рег (к*Ре 2 )и оценены ДН^ для к*Рез и к*Ре 4 . С 
использованием методологии термодинамического моделирования (ТМ) определен 
состав жидкого железа при Р = ІО 5 Па и Т = 1900-3100 К с учетом атомов (к*Реі ) и 
кластера к*Рег; оценены содержания к*Рез и к*Ре 4 . 

1. ВВЕДЕНИЕ 

Из обзора [1] следует, что при 298-3145 К в газовой фазе над железом образование кла- 
стеров не наблюдается. Термохимические характеристики летучих и конденсирован- 
ных метастабильных кластеров железа в доступной литературе отсутствуют. Складыва- 
ется впечатление, что существование малых кластеров железа маловероятно. Для про- 
верки этого предположения целесообразно использовать данные работы [2], где для 
кластеров к*Рег-4 методами молекулярной динамики с использованием различных 
представлений о структуре частиц рассчитаны энергии связи в кластерах (Б, кДж/моль) 
(см. табл. 1) . 

Цели настоящей работы были следующими: - рассчитать термодинамические 
характеристики конденсированного кластера Ре 2 , - выполнить равновесные компью- 
терные эксперименты с использованием данных об атомах (к*Реі ) и кластера к*Ре 2 по 
определению состава расплава железа в широком интервале температур и оценить со- 
держание в расплаве к*Рез и к*Ре 4 . 


Таблица 1. Энергии связей атомов кластеров (Э) и оценочные величины ДН ^ кла- 
стеров железа 


Связи в кластерах 

Энергии связи (Э), 
кДж/моль, [2] 

Оценка ДН 2 98 , 
кДж/моль 

Ре - Ре 

189 

~ 189 (к*Ре 2 ) 

Ре - Ре 

192 

~ 253 (к*Ре 3 ) 

Ре - Ре 

248 

-315 (к*Ре 4 ) 
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2. МЕТОДЫ ИССЛЕДОВАНИЯ 

Для оценки Дн!| 98 (СЭО) кластеров использовали величины Э из табл. 1 . При этом счи- 
тали, что СЭО (к*Ре 2 ) — О (Ре-Ре), поскольку СЭО в этом случае тождественна энергии 
связи между атомами в расчете на 1 моль. Для к*Ре 3 и к*Ре 4 использовали соотношения 

СЭО (к*Ре 3 ) ~ СЭО (к*Ре 2 ) + 0.333 [О(Ре-Ре)] и ( 1 ) 

СЭО (к*Ре 4 ) ~ СЭО (к*Ре 3 ) + 0.25[О(Ре-Ре)], (2) 

где численные коэффициенты 0.333 и 0.25 отражают доли, приходящиеся на связи Ре- 
Ре в каждом из кластеров. Результаты оценки СЭО (к*Ре п ) приведены в табл. 1 . Из 
табл. 1 видно, что для образования кластеров требуется большая энергия, значительно 
превышающая величины СЭО (к*Ме 2 . 5 ), рассчитанные для метастабильных кластеров 
Ьі 2 - 5 , Си 2 , Ве 2 , М§ 2 , Са 2 , 8г 2 , Ва 2 , А1 2 , 8і 2 , 8і 2 , Ое 2 , 8п 2 , РЬ 2 [3], т.е. в среднем СЭО 
(к*Ре 2 . 4 ) больше в ~ 10-15 раз. Поэтому для количественной оценки состава модельного 

расплава железа, кроме атомов, был выбран кластер к*Ре 2 с наименьшей величиной 
СЭО. 

Рассчитанные по методам [4, 5] термохимические характеристики к*Ре 2 приве- 
дены в табл. 2 вместе с температурной зависимостью приведенной энергии Гиббса. 

Компьютерные эксперименты выполнены с использованием методологии тер- 
модинамического моделирования (ТМ), пакета АСТРА.4 с БД АСТРА. ВА8 [6] при 
1900-3200 К с шагом 10-100» и Р - ю'Па, Модельный расплав железа содержал атомы 
и кластеры Ре 2 , газовая фаза - Аг, Ре, Ре * и е — газ. 

3. РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ 

Содержание кластера в модельном расплаве железа по данным ТМ показано на рис. Из 
него видно, что с ростом температуры от 1900 К до 3100 К мол. % [к*Ре 2 ] изменяется 
от 3.12-10 до 0.026, т. е. содержание метастабильной частицы увеличивается, как это 
оыло ранее установлено для конденсированных кластеров в расплавах щелочных ме- 
таллов [7], кремния [8], алюминия [9] и меди [10]. Исходя из логичного предположе- 
ния, основанного на данных [7-10], что при одинаковых температурах х[(к*Ре п )] будет 

обратно пропорционально величинам СЭО (к*Ре п ) (см. табл. 1), оценены содержания 
кластеров Ре 3 и Ре 4 при 1900 и 3100 К: 

х[к*Ре 3 ], мол. % х[к*Ре 4 ], мол. % 

1900 К ~ 2.3Т0 4 ~ 7 3- 10" 4 

3100 К • — 0.019 ~ппіа 
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С учетом этих оценок можно считать, что расплав железа содержит даже при 
3100 К > 99.95 мол. % атомов, т. е. практически является моноатомным. Это приводит в 
результате испарения к моноатомному составу газовой фазы над расплавом железа, что 
хорошо согласуется с данными экспериментальных исследований [1]. 

Для модельного расплава железа рассчитаны Т и ДН кипения по методу, приве- 
денному в [6]. Эти характеристики согласуются с известной экспериментальной ин- 
формацией: 



по ТМ 

по [11] 

Ткип, К 

-3130 

3143 

ДНкип, кДж/моль 

~ 351 

350 


4, ЗАКЛЮЧЕНИЕ 

В результате проведенного исследования показано, что содержание конденсированных 
малых кластеров Ре 2 , Ре 3 и Ре 4 в жидком железе в температурном интервале существо- 
вания жидкой фазы не превышает 0.05 мол. %. Поэтому можно считать, что расплавы 
железа практически состоят из атомов. Полученные данные согласуются с известными 
экспериментальными наблюдениями. 
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АННОТАЦИЯ. Для закладки выработанных пространств комбинатом «Нориль- 
ский Никель» используются составы на основе портландцемента с добавлением ме- 
стного техногенного отхода - ангидрита. Эти закладочные смеси отличаются высо- 
кой стоимостью и низкой долговечностью. В этой работе рассматривается замена 
портландцемента на каустический магнезит. Составы на основе каустического маг- 
незита имеют более высокие прочностные характеристики и являются более устой- 
чивыми к сульфатной коррозии. Закладочные смеси на основе магнезиального вя- 
жущего позволяют максимально использовать техногенные отходы предприятия та- 
кие как ангидрит, металлургические шлаки и карфосидерит, что снижает не только 
стоимость закладочных смесей, но и снижает экологическую нагрузку на регион. 

1. ВВЕДЕНИЕ 

К настоящему времени из недр земли извлекается до 15 млрд, т горной породы в год, из 
которой в хозяйственный оборот вовлекается только 5 млрд, т, остальные две трети ос- 
таются в отвалах [1]. Необходима перестройка традиционных технологических про- 
цессов и методов, которая приведет к созданию и внедрению производств с высокими 
технико-экономическими показателями, исключающих накопление техногенных мате- 
риалов и снижения уровня антропогенного воздействия на окружающую среду. 

К составам, изготовленным на основе техногенного сырья, предъявляется спектр 
повышенных требований, которые должны обеспечивать долговечность и безопас- 
ность материалов. Для решения этой проблемы необходимо установление оптимальных 
интервалов соотношений компонентов, важных для обеспечения технологических ас- 
пектов приготовления, транспортировки и получения требуемых физико-механических 
характеристик и свойств материалов. При этом техногенные материалы должны ис- 
ключать вторичные процессы, приводящие к деструкции затвердевшего материала. 
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Десятки лет металлургическая и химическая промышленность в г. Норильске 
наряду с полезной продукцией вырабатывают техногенные материалы, включающие 
металлургические шлаки и шламы гальванических производств, вскрышные породы, 
такие как ангидрит, каустический магнезит. Одновременно при разработке полезных 
ископаемых образуются выработанные пространства, которые требуют заполнения за- 
кладочными смесями. Производство закладочных работ призвано обеспечить безопас- 
ную и эффективную эксплуатацию месторождения и предотвратить самопроизвольное 
обрушение налегающих толщ и земной поверхности. Закладка выработанного про- 
странства является неотъемлемой технологической операцией в общем процессе добы- 
чи руды в горно-металлургической компании «Норильский никель». 

Для закладки выработанных пространств в настоящее время компанией «Но- 
рильский никель» используется в качестве вяжущего портландцемент. Несмотря на 
низкие требуемые прочностные характеристики закладочных смесей, идет большой 
расход цемента для обеспечения требуемых реологических свойств. Существующая 
технологическая схема приготовления закладочной смеси приведена на рис. 1. 

Закладочная смесь транспортируется в выработанное пространство по закладоч- 
ным трубопроводам в самотечном режиме. Как видно из технологической схемы в су- 
ществующем составе присутствуют портландцемент и ангидрит, которые согласно [2] 
способствуют образованию эттрингита. В качестве наполнителя закладочной смеси ис- 
пользуются шлаки - гранулированные отходы металлургического производства. Необ- 
ходимо отметить, что в шлаках, используемых для приготовления закладочных смесей, 
присутствуют сернистые соединения железа (табл. 2). В работе [3] отмечается, что при 
наличии этих соединений железа в техногенных отходах вследствие протекающих хи- 
мических реакций образуются гипс, который в сочетании с портландцементом также 
приводит к образованию эттрингита. Известно [4], что эттрингит, образующийся в за- 
твердевшем цементном камне, создает внутренние напряжения и нарушает целостность 
в композиционном материале. 



Рис. 1 . Технологическая схема производства закладочных смесей 
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Таким ооразом, замена портландцемента на другой вид вяжущего целесообразна 
не только с экономической точки зрения, но и для исключения деструктивных процес- 
сов при эксплуатации затвердевшей закладочной смеси. 

Учитывая необходимость устранения из состава закладочной смеси портландце- 
мента, была разработана смесь, включающая в качестве связующей матрицы ангидрит 
и каустический магнезит промышленного производства. Необходимость введения кау- 
стического магнезита в состав связующей матрицы основана на выводах, сделанных в 
работе [5]. Автор этой работы показал, что в такой композиции достигается обеспече- 
ние аморфной структуры новообразований, которые на этапе транспортировки закла- 
дочной смеси улучшают реологические свойства закладочной смеси, а при затвердева- 
нии, как установлено в работе [6], способны подвергаться структурированию, обеспе- 
чивая повышенные механические характеристики твердеющей композиции. 

С целью исключения усадочных явлений при твердении закладочной смеси в бе- 
тонную смесь, используемую для закладки выработанных пространств, вводился гра- 
нулированный шлак. Затворение смеси производилось раствором хлорида магния. В 
качестве структурирующей добавки использовался карфосидерит. Для улучшения пла- 
стичности бетонной смеси применялся супер пластификатор С-3. 

2. МАТЕРИАЛЫ 

Для приготовления закладочных смесей использовались составляющие, обеспечиваю- 
щие вяжущие свойства - ангидритовая порода Горозубовского месторождения (табл. 1) 
и каустический магнезит. В качестве заполнителя применялись шлаки (табл. 2), кото- 
рые измельчались в мельнице при совместном помоле с другими компонентами. Для 

затворения закладочной смеси использовался раствор бишофита с плотностью 1.16 
г/см 3 . 

Таблица 1 . Химический состав ангидритовой породы 


Содержание оксидов элементов, % 


Н 2 0 

8Ю 2 

А1 2 0з 

Ре 2 0з 

СаО 

М§0 

50 3 

1.08 

6,32 

1,78 

1,35 

35,42 

2,97 

45,25 


Минералогический состав ангидритовой породы, %: СаС0 3 + М§С0 3 - 12 
СаЗСД + Са80 4 -2Н 2 0 - 78, глинистые минералы - 9,5. 

Химический состав каустического магнезита: содержание М§0 не менее 75% 
8 і 0 2 не более 4,5%, СаО не более 3,5%. Влажность - 1 ,5%. 
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Таблица 2. Химический состав шлаков 


Вид шлака 

Содержание компонентов, % 

Ре 

8 

8Ю 2 

СаО 

М§0 

А1 2 0 3 

Г ранулированный шлак Надеждинского 

металлургического завода (НМЗ) 

40,4 

1.11 

35.4 

3.1 

1.59 

3.9 

Гранулированный отвальный шлак Ни- 

кель завода (НЗ) 

33,1 

1.38 

37.5 

3.73 

2.37 

5.76 


Минералогический состав (рис. 2а, б) гранулированных шлаков НМЗ, НЗ пред- 
ставлен фаялитом (железистым оливином). Шлак НЗ (рис. 26) наряду с кристалличе- 
ской фазой имеет аморфную составляющую. Дисперсность гранул шлаков НМЗ со- 
ставляет 2-5 мм, шлаков НЗ - 1-3 мм. 

Для обеспечения структурирования связующей матрицы применялся дисперс- 
ный отход гальванического производства сагГозМегйе (НзО)Ре 3 (8С>4) 2 -(ОН)б с размером 
частиц меньше 6 мкм. Порошок карфосидерита нерастворим в воде, состоит из кри- 
сталлических частиц (рис. За. б). Химический состав карфосидерита приведен в табл. 3. 


а 


б 




Рис. 2. Рентгенофазовый анализ металлургических шлаков: (а) - гранулированного 
шлака НМЗ, (б) - гранулированного отвального шлака НЗ 

Таблица 3. Химический состав карфосидерита 


Химические компоненты 

Ре 2 0з 

Н 2 0 

80 3 

Содержание, % 

49,83 

16,86 

33,31 
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а 


б 




Рис. 3. Рентгенофазовый анализ карфосидерита (а), микроструктура карфосидери- 
та (б) 

Основные требования к закладочному составу - удовлетворительная кинетика 
твердения закладочного бетона и реологические свойства, обеспечивающие норматив- 
ные прочности и транспортировку смесей по подземным трубопроводам и укладку в 
выработанные пространства. Подвижность закладочной смеси должна соответствовать 
осадке стандартного конуса равной 12 см, прочность бетона после твердения закладоч- 
ной смеси в возрасте 28 суток должна быть не ниже 4,5 МПа, в возрасте 180 суток - не 
менее 10 МПа. 

3. ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ 

Микроскопические исследования составов показали значительные изменения микро- 
структуры при добавлении структурирующей добавки - карфосидерита (рис. 4). При 
отсутствии структурирующей добавки ангидрит затвердевает в кристаллической форме 
в виде блоков-пластинок (рис. 4а), введение карфосидерита с суперпластификатором - 
С-3 приводит к аморфизации структуры и образованию структурированного межфазно- 
го слоя (рис. 46) на поверхности частиц шлака (рис. 4в). 

При этом отмечается существенное повышение прочности композиции, с дос- 
тижением максимального показателя прочности при сжатии 1 10 МПа в возрасте 28 су- 
ток. 
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а 


б 


в 



Рис. 4. Микроструктура закладочной смеси: (а) - при отсутствии карфосидерита, 
(б) - фрагмент структурирования связующей матрицы частицами карфоси- 
дерита, (в) - структура граничного слоя на поверхности частиц шлака 



Рис. 5. Рентгенофазовый анализ композиции в возрасте 2 месяцев с вяжущей мат- 
рицей на основе каустического магнезита (80 %) и ангидрита (20 %), затво- 
ренных хлоридом магния 
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Рентгенофазовый анализ закладочной смеси в возрасте двух месяцев показал 
образование оксихлоридов магния М§зС1(0Н)5-4Н 2 0 и наличие оксихлоридов магния, 
подвергнутых карбонизации М§ 4 (0Н) 2 СІ2(С0 3 )-6Н 2 0 (рис.5). В процессе твердения в 
возрасте девяти месяцев в спектре закладочной смеси отмечается значительное сниже- 
ние интенсивности отражений оксихлоридов магния и существенное возрастание ин- 
тенсивности карбонизированных оксихлоридов магния (рис. 6). Снижается также со- 
держание оксида магния. 

Учитывая отсутствие необходимости в такой высокой прочности закладочной 

смеси возможно увеличение содержания техногенных отходов в качестве заполните- 
ля. 

В составе закладочной смеси предполагается максимальное использование ме- 
стных техногенных материалов. В частности, в экспериментах в качестве связующей 
матрицы использовался состав, включающий 80 % молотого ангидрита и 20 % каусти- 
ческого магнезита. Для структурирования матрицы применялся карфосидерит 
(Нз0)Рез(80 4 )2-(0Н)б в количестве 1 % от общей массы, в качестве поверхностно- 
активной добавки, способствующей равномерному распределению карфосидерита в 
ооъеме матрицы и, одновременно, для повышения пластичности смеси, применялся су- 
перпластификатор С-3 в количестве 0,3 % от общей массы связующей матрицы. 



ренных хлоридом магния 


ПО 
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БЕСЦЕМЕНТНАЯ КОМПОЗИЦИЯ ДЛЯ ЗАКЛАДКИ ВЫРАБОТАННЫХ ПРОСТРАНСТВ НА 

ОСНОВЕ ТЕХНОГЕННЫХ МАТЕРИАЛОВ 


Расход раствора бишофита устанавливался, исходя из подвижности смеси, соот- 
ветствующей 22 см расплыва теста на приборе Суттарда, что обеспечивает технологи- 
ческие требования к пластичности закладочной растворной смеси. 

Таким образом, разработанная связующая матрица на основе техногенных мате- 
риалов горно-металлургической компании «Норильский никель» способна обеспечить 
необходимые физико-механические характеристики закладочной смеси. Отсутствие в 
составе матрицы алюминатной составляющей предопределяет ее устойчивость против 
сульфатной коррозии. Использование в качестве исходных компонентов только техно- 
генного сырья позволяет существенно улучшить экологическую обстановку, так как 
отходы производства вновь возвращаются в выработанные пространства в виде закла- 
дочных смесей и снижают экологическую нагрузку в районе города Норильска. 

4. ЗАКЛЮЧЕНИЕ 

На основе полученных результатов может быть решена проблема утилизации техно- 
генных отходов производственного комплекса «Норильский Никель». Замена порт- 
ландцемента на каустический магнезит целесообразна не только для исключения дест- 
руктивных процессов при эксплуатации затвердевшей закладочной смеси, но и для по- 
лучения более высоких прочностных характеристик. Использование разработанных со- 
ставов на основе магнезиального вяжущего дает значительное снижение стоимости за- 
кладочной смеси и позволяет решить проблемы защиты окружающей среды в окрест- 
ностях г. Норильска. 
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81ЛѴШАКУ. ТЬе иза§е оГ саизііс та§пе 2 ІСе іпзСеас! оГ СЬе рогСІапЬ сетепС Гог сЬе 
ргерагаСюп оГ СЬе Гі11іп§ тіхСигез із ехреЬіепС поС опіу Гог СЬе ехсіизіоп оГсЬе сІезСгисСіоп 
ргосеззез оГ СЬе ЬагЬепеё тіхСиге с1игіп§ іС’з ехріоісасіоп ЪиС Гог ІС’з есопошісаі изеГиІпезз 
аз ѵѵеіі. ТЬе еІаЬогаСеЬ іп СЬіз \ѵогк Ьіпсііп§ таСгіх іЬе таіп сотропепСз оГѵѵЬісЬ сЬеге аге 
Ше лѵазіе таіегіаіз оГ іііе тіпіп§ апсі теСа1иг§іса1 ргосеззез ас сЬе сотрапу “Могіізк 
Міскеі епаЫез СЬе епЬапсе СЬе рЬузісаІ апсі тесііапісаі сЬагасІегізСісз оГ сЬе Ьагсіепесі 
Й11іп§ тіхСиге. ТЬе аЬзепсе іп сЬе Ьіпсііп§ таСгіх оГ СЬе аІіитіпаСе сотропепС із 
сопёіСіопіп§ іС’з ёигаЬіІіСу Со СЬе зиІрЬаСе соггозіоп. ТЬе еІаЪогаСес! іп сЬіз ѵѵогк Ш1іп§ 
тіхСиге сопзшпез СЬе ''Ыогіізк МіскеГ сотрапу ѵѵазСе. ТЬегеГоге СЬе гезиІСз оГ сЬіз ѵѵогк 
аге ѵегу ітрогСапС Гог СЬе зо1ѵіп ё оГ СЬе ргоЫет оГ СЬе сіесгеазе оГ СЬе епѵігоптепС 
роІІиСіоп іп СЬе Могіізк. 
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КВАНТОВО-ХИМИЧЕСКИЕ РАСЧЕТЫ 


РКОВРЕМ8 ОР ТНЕ СКЕАТIОN ОР ^Г<ЮКЕАСТОК8 РОК ТНЕ 
8У1ЧТНЕ8І8 ОР МЕТАРРІС^ШРАКТІСРЕ8 ЕЧ САКВОN 8НЕРЬ8 

А.М. ЬІРАРЮѴ, Ѵ.І. КСЮОЬОѴ, N2/. КНОКНКІАКОѴ, А. А. ЭГОІК, 
Ѵ.Ѵ. КСЮОЬОѴ А, РТУ. 8ЕМАКША 


Вазіс Кезеагск - Ні§к екисаііопаі сепіег о! 1 скетісаі ркузісз апк тезозсору 
Ыктигі зсіепііГіс сепіег, Ька1 Оіѵізіоп, Киззіап Асакету ок Зсіепсез, ІгЬеѵзк, 
Киззіа 

е-таіі: коко1@ізіи.икт.ги 


АВ8ТКАСТ. Рог іке зупікезіз оі' теіаіііс папорагіісіез апсі папоѵѵігез іп сагЬоп зкеііз, іі 
із ехрекіепі Іо изе іке теіаі іопз ікаі сап Ье гесіисесі іп папогеасіогз Ьу Іке іпіегасііоп 
ѵѵіік іке пе§аііѵе1у скаг§ек ѵѵаііз оі' папогеасіогз. Тке зупікезіз із а гекисііоп-охікаііоп 
ргосезз, ікегекоге аготаііс кукгосагкопз ог ог§апіс роіутегз іо§еікег ѵѵіік теіаі заііз аге 
гекисііоп-охікаііоп раігз. Ѵѵіікіп іке кгатеѵѵогк оі" Напгее-Роск ікеогу ѵѵіііі іке тіпітаі 
Ьазіз 3-210, іке епег§іез о к тоіесиіез апк тоіесиіаг кга§тепіз аге саісиіаіек. Іп ікіз еп- 
ег§у токек іке ериіІіЬгіит аіотіс §еотеігу ок іке кга§тепІз апк іке іпіегасііоп епег§у 
ок іке геасііоп сепіегз аге сіеіегтіпесі. 

Тке токе1іп§ гезиііз аге геаіігек іп іке ехрегітепіз зіагііп§ \ѵіік тіхіп§ теііз ог $о- 
Іиііопз ок іке Зк шеіаі скіогікез (Сг, Мп, Ре, Со, Ыі, Си) ѵѵіік аготаііс кукгосагкопз ог 
роіуѵіпуі аісокоі зоіиііоп. 

Неаііп§ ок іке геасііоп таззез ир іо 400 С із саггіесі оиі іп а зіер-ііке таппег. Тке 
сопігоі оѵег іке ргосеззез із регкогтек изіп§ зресігаі, тісгозсоріс апк скетісаі теікокз. 
Хапоргокисіз аге зіикіек ѵѵіік іке изе ок Х-гау ркоіоеіесігоп зресігозсору, ігапзтіззіоп 
еіесігоп апк зсаппіп§ еіесігоп тісгозсору, кіккгасіотеігу апк іке теікок ок еіесігоп ті- 
сгокіккгасііоп. Оерепкіп§ оп іке сопкіііопз ок іке ргосезз, §іапі киііегепез апк іикиіепез 
іпсіизіѵе ок теіаіз апк теіа!-сопіаіпіп§ сотроипкз аге когтек. 


ШТКОБРСТЮХ 

Оигіп§ Іазі уеагз, скетізігу іп папогеасіогз із еѵо1ѵіп§ гарікіу. Аі ргезепі, іке іегт 
"папогеасіог’ із ге§агкек іп а ѵѵікег ітріісаііоп аз іке Гіеік ѵѵіік опе ок іке Ііпеаг рагатеіегз 
Іезз ікап 1 тст ѵѵкеге іке зупікезіз ок папорагіісіез ог ікеіг геасііопз сап оссиг. Ііпкег ікіз 
кеРтіііоп, іке іегт 'папогеасіог’ тау іпсіике іке кекесі ге§іопз ок теіаі заіі роіусгузіаіз; 
іпіегіауег ге§іоп іп теііз апк зикзіапсез ѵѵіік а Іауег зігисіиге; ехіепкек саѵіііез когтек Ьу 
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РІ.Ѵ. 5ЕМАКША 

тасго-тоіесиіез іп §е!з ог зоіиііопз оР роіутегз. ТЬе гедиігетепі Іо сіеѵеіор есо!о§іса11у 

сіеап ргосіисііопз \ѵі11 §іѵе ап оррогіипііу оР ѵѵісіе арріісаііоп оР папогеасіогз іп сЬетізІгу 
апсі теіаі 1иг§у. Ноѵѵеѵег, іі із сііРРісиІІ Іо ітріетепі сіігесіесі зупіЬезез іп папогеасіогз ѵѵііЬ- 
ои! ргеіітіпагу саісиіаііоп ехрегітепі ргесіісііпо *е ЬеЬаѵіог оРсЬетісаі рагіісіез апсі іЬеіг 
ІгапзРогтаІіопз іп папогеасіогз. Рог оЪіаіпіп§ теіаіііс папорагіісіез апсі папоѵѵігез іп сагЬоп 
зЬеІІз, іі із ехресііепі Іо еѵаіиаіе ІЬе роззіЬіІіііез оР гесіисііоп-охісіаііоп соиріесі геасііопз 
\ѵіШ рагіісіраііоп оР іЬе іопз оР теіаіз апсі ог§апіс сотроипсіз; іЬеп ІЬе гезиііз оР тосіе1іп§ 
зЬоиЫ Ье ехрегітепіаііу сопРтпесі. ТЬе аготаііс зисЬ аз ЬусігосагЬопз парЫЬаІепе, апіЬга- 
еепе, рЬепапіЬгепе, апсі роіушегіс Рипсііопаі сотроипсіз Ике роіуѵіпуі аІсоЬоі аге оРРегесі 
іог а зоигсе оР сагЬоп. Рог ІЬе таігісез ѵѵііЬ папогеасіогз (ѵѵііЬ саѵіііез ог іпіеііауег ге§іопз 
ѵѵЬеге ІЬе Рогтаііоп оР папорагіісіез Іакез ріасе) ѵѵе Ьаѵе сЬозеп еиіесііс теііз оР Зсі теіаі 

сЫогісіез (Сг, Мп, Ре, Со, N1, Си), §е1з оР роіуѵіпуі аісоЬоі апсі роіуѵапасііит асісі апсі ііз 
сіегіѵаіез іпіегсаіаіесі ѵѵііЬ теіаіз. 

МАТНЕМАТІСАЬ МООЕЬШС ОР ^ІМОКЕАСТОК РКОСЕ88Е8 

МаіЬетаІісаі тосіе1іп ё ѵѵаз сопсіисіесі іп зеѵегаі зіерз. Іп соппесііоп ѵѵііЬ іЬе іазк Іо 
сіеіегтіпе іЬе роззіЬііііу оР ІЬе Рогтаііоп оР §1оЬиіаг ог суііпсіег папорагіісіез, аі іЬе Рігзі 
зіер, ІЬе Ргаотепіз оР іп1егас1іп§ рагіісіез апсі асііѵе сепіегз оР іЬе папогеасіог ѵѵаііз ѵѵеге 
зіп§1есі оиі. Іп сазе оР аготаііс ЬусігосагЬопз, іі ѵѵаз зи§§ез1есі іЬаі іп ІЬе теіаі сЫогісіе 
таігіх, ІЬе тоЬШіу оі ІЬе теіаі іопз іп іЬе саѵіііез оР папогеасіогз ѵѵііЬ іЬе ѵѵаііз розіііѵеіу 
сЬаг§есі із зиРРісіепІ Рог ІЬе Рогтаііоп оР іЬе сотріехез ѵѵііЬ аготаііс писіеі оР сотроипсіз, 
“сігаѵѵп” іпіо іЬезе саѵіііез. Іп іЬе саісиіаііопз, іЬе тоіесиіез оР Ьепгепе апсі парЫЬаіепе 
ѵѵеге изесі Рог іЬе тоіесиіез ітііаііп» іЬе аготаііс ЬусігосагЬопз изесі іп іЬе ехрегітепі. ТЬе 

тесЬапізт оРіЬе тіегасііоп апсі іЬе еѵоіиііоп оРіЬе сІеЬусігороіусопсіепзаІіоп аге сіізріауесі 
іп Рі§. 1. 

Іп іЬіз тоёеі егозіоп’ оРіЬе папогеасіог ѵѵаііз із зи§§ез1есі ассотрапіесі Ьу іЬе НЬ- 
егаііоп оРЬус1го§еп сЫогісіе апсі ііз рагііаі аЬзогрІіоп іп іЬе таігіх. 


а 



Ь 






с 


Еі§. 1. ЗсЬете оР іпіегасііоп оі Ьепгепе ѵѵііЬ теіаіз: а) гіп§ зігисіиге оРЬепгепе, 

Ь) сотріех теіаі / Ьепгепе гіп§ ѵѵііЬ Си, с) сотріех теіаі / Ьепхепе гіп§ ѵѵііЬ 
Со апсі № 
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РЯ0ВЬЕМ8 ОР ТНЕ СЯЕАТЮ1Ч ОР ЫАШЯЕАСТОЯ8 РОЯ ТНЕ 8ѴПТНЕ8І8 ОР 
МЕТАЫЛСКАМ0РАЯТІСЕЕ8 Ш САЯВОЫ 8НЕЬЬ8 


Тііе тоёе1іп§ оЕ (Не ргосеззез оссиггіп§ іп (Не §е1з оЕ роіуѵапаёіит асіё апсі і(з ёе- 
гіѵаііѵез ёигіп§ (Не репе(га(іоп оЕ (Не роіуѵіпуі аІсоНоІ тоіесиіез іп(о (Неіг іпіегіауег зігис- 
Шге \ѵаз сопёис(её изіп§ іЬе Ега§теп( оЕ роіуѵапаёіит асісі апё еіЬапоІ ог Нехапе(гіо1 1, 3, 
5. Аз (Ніз (оок ріасе, (Не ргосеззез оЕ ёеНуёгаёоп апё ёеНуёго§епа(іоп \ѵеге еѵа1иа(её; і( 
§аѵе (Не Ьазіз Еог (.Не Еогтаёоп оЕ роіуепе сНаіпз апё сагЬопігеё ІиЬиІепез (На( аге 1а(ег 
сопёепзеё. 

Рог (Не циап(ит-сНетіса1 са1си1а(іопз оЕ (Не ргосеззез іп (Не роіуѵіпуі аІсоНоІ §е1з, 
(Не іп(егас(іоп оЕ (Не те(а1 іопз Егот (Не ге1а(её сЫогіёез ѵѵаз зи§§ез(её; (Неп (Неу репе- 
(га(её іп(о пе§а(іѵе1у сНаг§её саѵі(іез оЕ (Не папогеас(огз Еогтеё іп (Не роіуѵіпуі аІсоНоІ 
§е1з соп(аіпіп§ аІсоНоІ Ега§теп(з; Нехапе(гіо1 1, 3, 5 \ѵаз зе1ес(её Еог (Не аІсоНоІ. 

АЕ(ег зе1ес(іоп апё зиЬз(ап(іа(іоп оЕ (Не Ега§теп(з Еог (Не тоёе1іп§, (Не зеті- 
етрігісаі яиап(итсНетіса1 са1си1а(іопз ѵѵеге сопёис(её. А( (Ніз з(ер, (Не шіпішаі Ьазіз 
МІМ ог 3-210 \ѵаз изеё, апё (Не зеагсН оЕ ециіІіЬгіит §еоте(гіса1 сопГі§ига(іоп оЕ (Не 
Еогтеё сотріехез \ѵаз регЕогтеё Еог апу ргоіес(іоп оЕ (Не епег§у §гаёіеп(. 

А( (Не пех( з(ер, (Не са1си1а(іопз \ѵеге сопёис(её ѵѵё(Н (Не Неір оЕ (Не ех(епёеё Ьазіз 
Т2Ѵ апё аізо (Не те(Ноё 2ШЭ0/1 изіп§ (Не зоЕ(\ѵаге ргоёис( НурегсНет апё (Не сотри(ег 
(есНпо1о§у оЕ (Не Віо-ТесНпо1о§у Іпз(і(и(е, Огапаёа ЕНііѵегзку. 

А( (Не Еоиг(Н з(ер, (Не Еігз(ргіпсір1ез саісиіайопз \ѵі(Н (Не изе оЕ НаіЕгее-Роск те(Ноё 
апё (Не зойѵѵаге сотріех ОАМЕ88 \ѵеге саггіеё ои(. А( (Не зате (іте, (Не ас(іѵа(іоп епег§у 
апё (Не зресіЕіс га(е оЕ (Не ргосезз ѵѵеге еѵа1иа(её; (Неу ѵѵеге изеё Еог ёе(егтіпа(іоп оЕ (Не 
ргеЕегеп(іа1 рго§гезз оЕ (Не геас(іопз. ТНе роззіЬіІкіез оЕ (Не Еогтаёоп оЕ Еиііегепез апё (и- 
Ьиіепз ѵѵеге ёеГтеё іп (Не Егате\ѵогк оЕ(Не тоіесиіаг ёупатісз тоёеі ММ+ [1]. 

Рог (Не Рігз( сазе, (Не С-Н Ьіпёіп§ епег§іез ѵѵеге са1си1а(её аз (Не ёіЕЕегепсе оЕ (Не еп- 
ег§у оЕ (Не сотріех іп (Не ециіІіЬгіит сопГі§ига(іоп апё і(з епег§у іп (Не аЬзепсе оЕ Нуёго- 
§еп ((Не епег§у оЕ Н \ѵаз (акеп еципі (о хего): 

Ес-н = Е(С 6 Н 6 М) п+ - Е(С 6 Н 6 М) <п ' 1)+ ; 

апё Еог (Не ез(іта(іоп оЕ (Не те(а1 са(іоп іп(егас(іоп (Не Ео11о\ѵіп§ Еогтиіа \ѵаз изеё: 

Ем-сбНб = [Е(СбН 6 ) + Е (М п+ )] -Е(С 6 Н 6 М) п+ . 

Аз і( Еоііоѵѵз Егот ТаЫез 1 апё 2, ѵѵНеге (Не Ьіпёіп§ епег§іез оЕ С-Н апё оЕ (Не 
те(а!-агота(іс гіп§ аге §іѵеп, (Не Ьез( з(іти!ап( Еог (Не Нуёго§еп ёе(асНтеп( Егот (Не гіп§. 
саизіп§ (Не ёіз(игЬапсе оЕ і(з сор1апагі(у, із ап іоп оЕ Ьіѵа1еп( піскеі. Іп (Ніз сазе, (Не гіп§ 
асцшгез (Не Еогт зітііаг (о (Не сопЕогтаёоп оЕ (Не сусІоНехапе “Ьа(Н”. 

ТНе гезиКз оЕ (Не са1си1а(іопз аге §іѵеп іп ѵѵогкз [2 - 4]. ТНе тоёе1іп§ гезиПз аі- 
Іоѵѵеё (о саггу ои( (Не ехрегітеп(з Еог оЬ(аіпіп§ те(а1 папорагёсіез апё папо\ѵігез іп (Не 
сагЬоп зНеІІз \ѵі(Н а сотрагайѵеіу НШе сопзитр(іоп оЕ епег§у. 

Рог (Не тоёе1іп§ оЕ (Не папораііісіез зуп(Незіз Егот роіуѵіпуі аІсоНоІ іп(о роіуѵа- 
паёіс асіё ёегіѵаёѵез аз папогеасГогз (Не Ега§теп(з оЕ роіуѵапаёіс асіё апё е(Напо1 аге изеё. 
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ТаЫе 1. ТЬе ЪіпШп§ епегду оР С-Н (Найгее) 


ТЬе зсЬете оР Ше Ьус1го§сп ОеіасЬтепі Епег§у 
С б Н 6 С 6 Н 5 ‘ + Н + 0,674485 

(С 6 Н 6 Си) + С 6 Н 5 Си + Н + 0,50 1 245 

(С 6 Н 6 Со) + С 6 Н 5 Со + Н + 0,498053 

(С 6 Н 6 Мп) + С 6 Н 5 Мп + Н + 0,498 1 26 

(С 6 Н 6 №) + С 6 Н 5 № + Н + 0,305434 


ТаЫе 2. ТЪе епег§у оТ Ше теіаіаготаііс гіп§ (Найгее) 


ТЬе зсЬете оТ Ше іпіегасііоп 
(С 6 Н 6 М) П+ С 6 Н 6 + М п+ 
(СбНбСи) С 6 Нб + Си + 
(СбНбСо) + С 6 Н 6 + Со + 
(СбНбМп) + С 6 Н 6 + Мп + 
(СбНбНі) 2 СбНб + № 2+ 


Епег§у (Т2Ѵ) 
Т2Ѵ 
-0,05210 
-0,05588 
-0,05392 
-0,22156 


ТЬе ригрозе оТ Шіз іпѵезп 8 аііоп із Іо сіеіегтіпе Ше аІсоЬо! ЬеЬуШаііоп роззіЪІІІіу ѵѵіШ Ше 
Ро11оѵѵіп§ ёеЬуШо§епа1іоп апё роіутегігаііоп Шігіпё Шеіг іпіегасііоп ѵѵіШ Ше Іауегз оТ ѵа- 
пасііит охіОе зузіетз, оп Ше ехатріе оР ѵапаШит охіёе апё ѵапасііп асіёз агкЗ еШапоІ. 

ТЬе Шгесі тос!е1іп§ оР сЬетісаІ іпіегасііоп оР еШапоІ тоіесиіе ѵѵіШ сотріехез соп- 
Іаіпіп§ Іѵѵо Ѵ205 ругашіёз із аізо саггіесі оиі. ТЬе сазе ѵѵЬеп аіі оху§еп аіошз аге Ьоипё 
\ѵіШ Ьу<Зго§еп аіошз іп Ше сотріех із сопзіёегеё (Рі§. 2), аз \ѵе11 аз Ше сазе ѵѵЬеп еасЬ 
ругатісі Ьаз опе зіорреО сЬетісаІ Ьопсі оп оху§еп аіот ѵѵіШ зузіет сЬаг§е —2 (Рі§. 3). 

Оигіп§ Ше тосіеі ехрегітепі зЬоѵѵп іп Рі§. 3 Іѵѵо ЬуШо§еп аіотз Ргот еШуІ аІсоЬо! 
тоіесиіе ІгапзТег Іо оху§еп аіот оР опе оР ѵапаШит ругатіОз апсі іо оху§еп Ьгіс1§е аіот. 

ТЬіз ргоѵез Ше роззіЬіІііу оР еШуІ аІсоЬоІ 0еЬуШо§епаІіоп іп Ше іпіегасііоп ргосезз ѵѵіШ 
Ріа! Іауег оР ѵапаШит охісіе. 


ТНЕ РОКМАТІОХ ОР МЕТАЬ ХАІѴОРАКТІСЕЕ8 АХБ І\АІЧО\ѴІКЕ8 IX ТНЕ 
САКВОХ 8НЕЬЕ8 РКОМ ТНЕ 80ЫГП0Х8 АХБ МЕЬТ8 ОР ТНЕ МІХТІЖЕ8 ОР 
МЕТАЬ СНЕОКШЕ8 АШ ОКСАМС СОМРОШН8 

Рог опе оР Ше ѵагіапіз оР оЬіаіпіп§ папозігисіигез Ргот аготаііс ЬуШосагЬопз, Ше тіхіиге 
оР апіЬгасепе ѵѵаз изеО іп Ше еиіесііс тек соп!аіпіп§ аіитіпит сЫогіОе, зосііит (Ьеііег. 
ІііЬіит) сЫогМе апё сЫогіОе оР опе оР Ше Зсі теіаіз (Си, №, Со ог Мп). ТЬе Іетрегаіиге оР 
Ше теіі ѵѵаз кері іп Ше гап§е оР 300 -500 С 0ерепШп§ оп Ше ігапзіііоп теіаі изеШ ТЬе ііте 
оРЬо1Шп§ ѵѵаз езІаЫізЬес! Ьазесі оп Ше сЬап^е оР Ше оріісаі зресіга, Ше аЬзепсе оРШе НС1 


116 


ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. Том 7, № 1 


РКОВЬЕМЗ ОР ТНЕ СКЕАТЮИ ОР ИАИОКЕАСТОК5 РОК ТНЕ 5УИТНЕ5ІЗ ОР 
МЕТАиЛСИАИОРАКПСЬЕЗ Ш САКВОЙ ЗНЕЬЬЗ 



Гі§. 2. ЗсЬете ой еіЬапоІ іпіегасііоп ѵѵііЬ Ѵапасііит Охійе йга§тепі іп \-ѵЬісЬ Оху§еп 
аіогпз аге соппесіей ѵѵііЬ Нуйго§еп 




Гі§. 3. ЗсЬете ой еіЬапоІ іпіегасііоп ѵѵііЬ Ѵапасііит Охійе Йга§тепі Ьаѵіп§ іѵѵо поп 
сотрепзаіесі Оху§еп аіотз апй іЬе сЬаг§е -2 

ІіЬегаііоп апй іЬе арреагапсе ой Іііе ргасіісаііу зіаЫе зигйасе §еотеігу. АЙег іЬе епсі оі' іЬе 
ргосезз іЬе геасііоп тазз ѵѵаз ѵѵазЬей ѵѵііЬ ЬуйгосЫогіс асісі. Ьоі ѵѵаіег, асеіопе апсі огіЬо 
хуіепе йог іЬе зерагаііоп о Г Ше папоргойисі йот іііе тайіх. ТЬе папоргойисі оЬіаіпесІ ѵѵаз 
сіізрегзесі іп асеіопе іп іЬе иЙгазоипсі йеісі апсі Йіеп ѵѵаз іпѵез1і§аіей Ьу теапз оі - ійе ігапз- 
тіззіоп еіесігоп апсі зсаппіп§ еіесігоп тісгозсорез, Йіе теіЬойз оі - еіесігоп тісгойіййгасііоп 
апсі Х-гау рЬоіоеІесігоп зресйозсору. ТЬезе іпѵезІі§аІіопз геѵеаіесі папозігисіигез іп Йіе 
&гт ой ІиЬиІепез, йіііегепез апй папосгузіаіз. ТЬе оиіриі ой сагЬоп тиііі-ѵѵаіі папоіиЬез 
тсгеазез ѵѵЬеп іЬе ргосезз із саггіесі оиі оп іЬе ргосеззей теіаі йоііз. Іп Рі§. 4, опе сап зее 
тісго§гарЬз ой папоіиЬез оЬіаіпесІ оп іЬе апойігей йоіі. 

Рог іЬе зупіЬезіз ой папорагіісіез апй папоѵѵігез йот ійе тіхіиге ой теіаі заііз апй 
роіуѵіпуі аІсоЬоІ (РѴА), іЬе а^иеоиз зоіиііопз ой заііз ѵѵеге тіхей іп а сегіаіп гаііо ѵѵііЬ ійе 
а^иеоиз зоіиііоп оі РѴА. ТЬе теап тоіаг гаііо ой РѴА іп іЬе тіхіиге ѵѵаз 5. ТЬе ехрегі- 
тепіз ѵѵеге саггіей оиі оп іЬе §1азз зиЬзігаіез; аЙег іЬе оЬіаіпей тіхіигез Ьай Ьееп йгіей, 
іЬеу йогтей соіогей ігапзрагепі йітз. Оп зоте затріез, іЬе йітз ѵѵеге Ьгокеп йие іо а 1аг§е 
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Гі§. 4. Місго§гарЕз оР тиііі-ѵѵаіі папоіиЬез оп Ае апоАгес! аіишіпиш Роіі: 
а) - Ае ТЕМ тісго§гарЕ; Ь) - Ае 8ЕМ шісго§гарЬ 

зигіасе іепзіоп. ТЕе йігпз ѵѵеге Ееаіесі аі і = 250 С ипііі Аеіг соіог, сотрозіііоп апсі тог- 
рЕо1о§у сЕап§еА Рог Ае сопігоі оѵег Ае ргосезз, а сотріех оР теАосІз ѵѵаз изесі, і.е. рЕо- 
іосоіогітеігу, оріісаі шісгозсору, Х-гау рЕоіоеіесігоп зресігозсору апсі аіотіс роѵѵег ші- 
сгозсору. 

\Ѵ1іеп РѴА \ѵаз асісіесі іо Ае росѵсіегз оР теіаі сЫогіёез, Ае соіог оР Ае тіхАге 
сЕап§еА Ае тіхіиге оР соррег сЫогісІе Ьесате уе11о\ѵ§гееп, Ае соЪаІІ сЫогіёе тіхАге - 
Ыие, апсі піскеі сЫогіёе - ра1е§гееп. ОЬзегѵіп§ Ае соіог сЕап§ез, опе сап Аа\ѵ а сопсіизюп 
Ааі ѵѵЕеп РѴА іпіегасіз \ѵіА теіаі сЫогісіез, Ае Рогтайоп оР сошріех сотроипёз іакез 
ріасе (Рі§. 5). 

Атоп§ Ае аЬоѵе-Азсиззес! шеіаіз, ігоп із тозі асііѵе. Вгоѵѵпгесі іпсіизіопз оп Ае 
рЕоіо§гарЕ еѵіёепсе Ае Рогтайоп оР Ае сотріех ігоп сотроипйз. Іп асЫійоп, оп аіі Ае 
рЕоіо§гарЕз сіерісііпр Ае тіхАгез сопіаіпіпе теіаі сЫогісІез, опе сап зее а пеі оР ѵѵеаѵез, 
ѵАісЬ аге тозі Іікеіу Ае гейесііопз оР папозігисіигез. 



Гі§. 5. ТЕе рЕоіо§гарЕз оРАе затріез сопіаіпіп§ РѴА апсі соррег сЕІогіёе (а), соЬак 
сЕіогіёе (Ь), апсі піскеі сЕІогісіе (с) 
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РЯОВЬЕМЗ ОР ТНЕ СКЕАТІСЖ ОР КАШЯЕАСТОЯЗ РОЯ ТНЕ 5УЖНЕ5ІЗ ОР 
МЕТАЫЛСКАРТОРАЯТІСЬЕЗ ІИ САЯВОЫ ЗНЕЬЬЗ 


Іп огёег Іо сотраге іЬезе зішсіигез, Йіе іпѵезІі§аІіопз оР Іііе шогр1іо1о§у оР іЬе йішз 
сЬап§іп§ оѵег а сегіаіп гап§е оР Іетрегаіигез \ѵеге саггіесі оиі ѵѵііЬ іЬе Ьеір оР аіотіс роѵѵег 
тісгозсору (Рі§. 6). 

\Ѵ1іеп Іііе папоргосіисі рісіигез оЬіаіпесі Ьу аіотіс роѵѵег тісгозсору апсі оріісаі ті- 
сгозсору аге сотрагесі ѵѵіііі Іііе ТЕМ тісго§гарЬ оР Іііе папоргосіисі Ігеаіесі іЬегтаІІу апсі 
ѵѵіііі ас]иеоиз зоіиііоп Тог Іііе таігіх гетоѵаі, опе сап поіісе зоте соггезропсіепсе Ъеіѵѵееп 
іЬет. ТЬе папоргосіисі; гергезепіз іпІегѵѵеаѵіп§ ІиЪиіепз сопІаіпіп§ Си(І), Си(ІІ). Іп Рі§. 7, 
гііеге аге аізо оріісаі еіТесІз іпсіісаІіп§ 1і§Ы роіагігаііоп аі 1і§Ьі Ігапзтіззіоп 11ігои§1і ІІіе 
Тіітз оѵѵіп§ Іо іЬе сіеРесІз арреагіпа ёигіп§ ІІіе Рогтаііоп оР іЬе сотріех сотроипсіз аі іЬе 
іпіііаі з1а§е оР ІІіе ргосезз. 

Оие Іо іЬе Расі, іЬаі теіаі іопз аге асііѵе, іп ІІіе роіутег тесііит іЬеу іттесііаіеіу 
арреаг іп ІЬе епѵігоптепі оР Іііе РѴА тоіесиіез апсі і'огт Ьопсіз ѵѵіііі ІІіе Ьусігохуі §гоирз оР 
1Ыз роіутег. Роіуѵіпуі аісоЬоі геріісаіез іЬе зігисіиге оР іЬе рагіісіе іЬаі іі зиггоипсіз; Ьоѵѵ- 



Гі§. 6. ТЬе Місго§гарііз оі' ІІіе зигРасе §еотеІгу оР іЬе РѴА Рііт ѵѵіііі піскеі сЫогісіе 



Рі§. 7. Місго§гарЬз оР ІиЬиіепез 
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еѵег, ёие іо іке іепсіепсу ок іке тоіесиіез ок іке теіаі заііз ог оікег теіаі сотроипёз іо 
сотЬіпе, РѴА аз ік епѵеіорз іке рохѵсіег рагіісіез, апсі ікегекоге іке когтз о Г іке оЬіаіпесі 
папозігисіигез сап Ье сііккегепі. Тке оріісаі тісгозсору теікосі аііоѵѵз Іо сіеіегтіпе іке 
зігисіиге о к іке папозігисіигез аі іке еагіу зіа§е. 

\ѴЬеп іке затріез аге кеаіеё, ёекусігаііоп оссигз, апсі аз а гезиіі, теіаі-сопіаіпіп§ 
папоіиЬез і'огт. Ткезе ргосеззез аге ікогои§Ыу ёезсгіЬесі іп ѵѵогкз [5-6]. Эекусігаііоп Іеасіз 
іо іке сіагкепіпе оі іке к 1т. АЙег іке затріез Ьаѵе Ьееп кеаіеХ оп іке ркоіо§гарк іке ге- 
таіпіп§ пеі ок лѵеаѵез сап Ье зееп, і.е. іке зігисіиге тогрко1о§у каз гетаіпеХ То зоте ех- 
іепі, ікіз касі іпсіісаіез ікаі іке іпіііаііу когтесі зігисіиге о к таігісез із іпкегііесі. Тке теікокз 
оі оріісаі зресігозсору апсі Х-гау ркоіоеіесігоп зресігозсору а11о\ѵ іо сіеіегтіпе Іке епегеу 
оі іке іпіегасііоп ок Іке скетісаі рагіісіез іп іке папогеасіогз ѵѵіік іке асііѵе сепіегз ок іке 
папогеасіог ѵѵаііз, ѵѵкіск зіітиіаіе гесіисііоп-охісіаііоп ргосеззез. 

Оерепсііп§ оп іке паіиге ок іке теіаі заіі апсі іке еіесігоскетісаі роіепііаі ок іке 
теіаі, сііккегепі теіаі гесіисііоп папоргокисіз іп іке сагкоп зкеііз ёіккегіп§ іп зкаре аге 
іоітесі. Вазек оп ікіз гезик \ѵе тау зреак аЬоиі а пе\ѵ зсіепііііс Ьгапск — папотеіа11иг§у. 

ССЖСЫШСЖ 

Тке роззіЫе лѵауз ког оЫаіпіп§ теіаіііс папозігисіигез іп сагЬоп зкеііз каѵе Ьееп сіеіег- 
тіпек. Тке іпѵезіі§аііоп гезиііз аііосѵ іо зреак аЬоиі іке роззіЬііііу окіке ізоіаііоп ок теіаі- 
ііс апсі теіаі-сопіаіпіп§ папо-рагіісіез іп іке сагЬоп зкеііз скккегіп§ іп зкаре апсі зігисіиге. 
Ноѵѵеѵег, ікеге аге зіііі ргоЫетз геіаіек іо іке саісиіаііоп апсі ехрегітепі Ьесаизе изіп§ іке 
ехізііп§ іпѵезіі§аііоп теікокз іі із кійісиіі ипатЬі§иоизіу іо езіітаіе іке §еотеігу апсі еп- 
ег§у рагатеіегз ок папогеасіогз ипсіег іке сопсііііоп ок ‘егозіоп’ ок ікеіг ѵѵаііз сіиппе іке 
когтаііоп ок теіаіііс папозігисіигез іп ікет. 
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ПЕРСОНАЛИИ 


К 70 - ЛЕТИЮ АЛЕКСЕЯ МАТВЕЕВИЧА ЛИПАНОВА 



Академик РАН 

Алексей Матвеевич ЛИПАНОВ 

3 марта 2005 г. исполняется 70 лет выдающемуся ученому, академику Российской 
академии наук Алексею Матвеевичу Липанову. 

А.М. Липанов родился в поселке Усть-Баргузин в семье рабочих. В 1954 г. он 
поступил на физический факультет Томского государственного университета, по 
окончании которого в 1959 г. направлен на работу в Научно-исследовательский 
химико-технологический институт. В это время началось создание твердотопливных 
ракетных двигателей, для которых не были исследованы многие проблемы внутренней 
баллистики, особенно начального периода работы двигателей. За решение этой 
проблемы и взялся А.М. Липанов. Она послужила стимулом для его исследований в 
области прикладной математики. По данной тематике в 1964 г. Алексеем Матвеевичем 
была защищена кандидатская диссертация «Газодинамические уравнения внутренней 
баллистики РДТТ и их решение методом характеристик». Проблема нестационарной 
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газодинамики затем получила развитие в его последующих работах и в работах его 
учеников. 

Развитие ракетной техники потребовало более глубокого понимания физико- 
химических процессов, протекающих в двигателях и топливе. К этому периоду 
относятся работы А.М. Липанова по расчету двухфазных течений в камерах и соплах 
РДТТ, процессам горения твердых топлив, созданию систем автоматизированного 
проектирования РДТТ, использованию статистических методов для оценки разбросов 
внутрибаллистических параметров двигателей. Обобщением этого цикла работ А.М. 
Липанова стали докторская диссертация «Некоторые проблемы внутренней баллистики 
РДТТ» и монография «Нестационарные и квазистационарные режимы работы РДТТ». 
В них впервые изложены подходы к математическому описанию сложных физико- 
химических процессов, протекающих при горении современных твердых топлив и 
разложении полимерных материалов. 

Следует отметить, что теоретические разработки А.М. Липанова отражали потребности 
развития современных ракетных систем на твердом топливе различного назначения, в 
создании которых он принимал непосредственное участие от стадии проектирования до 
летных испытаний. Среди двигателей, разработанных с его участием, были и 
двигатели, использованные в первом полете на Марс в 1971 г. За работы в области 
ракетостроения А.М. Липанову в 1985 г. была присуждена Государственная премия 
СССР. 

С 1976 г. начинается плодотворная деятельность Алексея Матвеевича в 
должности ректора Ижевского механического института. Он многое создал для 
развития материально-технической базы кафедр и факультетов института. Им открыта 
кафедра проектирования и производства двигателей летательных аппаратов, создано 
конструкторско-технологическое бюро «Восход», создан ряд других новых 
коллективов и специальностей. В этот период он также продолжает интенсивно 
заниматься научной деятельностью. Разработанная Липатовым А.М. универсальная 
автоматизированная система формирования областей, интегрирования при расчете 
параметров механических и термогазодинамических процессов позволила на порядок 
сократить затраты труда на составление программ при расчете процессов заполнения 
внутреннего объема камеры двигателя полимерными вязко-упругими составами, при 
исследовании напряженно-деформированного состояния элементов конструкций РДТТ. 
горения или течения продуктов сгорания. В этот же период им разработан новый класс 
двигателей для МГД-генераторов и газодинамических лазерных установок, так 
называемых двигателей перезарядки, позволяющих выполнять многократные запуски в 
течение заданного времени. 

Работая заместителем министра высшего и среднего специального образования 
РСФСР, А.М. Липанов явился инициатором создания мелкосерийных производств при 
вузах, открытия ряда НИИ и КБ, установления межгосударственных научных 
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контактов и укрепления материально-технической базы вузов. К этому периоду 
относятся его первые работы по мезоскопическим системам. В 1987 г. он избирается 
членом-корреспондентом Академии наук СССР. В г. Ижевске им создан Институт 
прикладной механики и Удмуртский научный центр Уральского отделения РАН, в 
рамках которого сегодня формируются коллективы будущих научно- 
исследовательских организаций. В этот период он активно занимается проблемой 
образования кластеров в процессах горения металлических частиц при 
переконденсации и дроблении материалов. Липановым А.М. создана теория 
молекулярно-динамического эволюционирования мезочастиц с использованием 
современных видов потенциалов межмолекулярного взаимодействия. Новым 
направлением в его деятельности становится проектирование пожаробезопасных 
материалов с использованием вспенивающихся материалов и ингибиторов процесса 
горения в сочетании с явлением образования графитоподобных веществ. Это позволило 
создать и экспериментально отработать образцы материалов, оригинальных по 
цветовой гамме и выдерживающих высокие тепловые нагрузки. 

В 2000 г. А.М. Липанов избирается действительным членом Российской 
академии наук. Среди достижений Алексея Матвеевича и его школы в 2000 - 2004 г. 
необходимо отметить следующие работы: 

- исследование конвективного и лучистого потоков и их влияние на нагрев и 
зажигание топливного заряда в РДТТ с застойными зонами; 

разработка моделей и исследование течений вязко-упругой полимеризующейся 
жидкости в технологических процессах ракетной техники; 

- разработка моделей и исследование напряженно-деформированного состояния 
полимерных композиционных материалов сложной пространственной конструкции; 
создание методов высокого порядка точности для прямого численного 
моделирования турбулентных течений сжимаемого газа и применение их для 
решения задач о турбулентных течениях газа; 

- разработка математических моделей технологических процессов в нефтяной и 
газовой отраслях и др. 

Созданной А.М. Липановым научной школой подготовлено более шестидесяти 
докторов и кандидатов наук. Алексей Матвеевич ведет активную научно- 
общественную работу, участвует в различных международных научных конференциях, 
в том числе в качестве председателя и члена оргкомитета, организовывая российские и 
международные встречи. За успехи в научной, научно-педагогической и 
организационной деятельности А.М. Липанов награжден многими государственными 
премиями и наградами. 

Редколлегия журнала сердечно поздравляет Алексея Матвеевича со славным 
юбилеем и желает ему больших творческих успехов в его дальнейшей деятельности. 
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РК0ГЕ880К СЕШАБѴ ЕЕКЕМОѴІСН 2АІКОѴ: 
НАЬЕ А СЕМШѴ Ш 8СIЕNСЕ 


РгоР. Еіі М. Реагсе 

Негтап Р. Магк Роіутег КезеагсЬ Іпзіііиіе РоІуІесЬпіс Шіѵегзііу 333 Зау 81., 
8іх МеІгоІесЬ Сепіег Вгоокіуп, ЫУС, 178 А 


Рапиагу 7, 2005 \ѵі11 Ье іЬе 70' 1 ' ЬігіЬсІау оРРгоР. О. Е. 2аікоѵ апсі 50 уеагз оРЬіз зсіепІіРіс 
асііѵііу/ 2аікоѵ ѵѵаз Ьогп іп Отзк, 8іЬегіа ѵѵЬеге Ье §гасіиаІе<і Ргот іЬеіг ргітагу, пжісПе. 
апО Ьі§Ь зсЬооІ. Не аізо ^гасіиаіесі Ргот а тизісаі зсЬооІ ѵѵЬеге Ье зіисііесі ѵіоііп. Ноѵѵеѵег, 
Ьіз рагепіз ЕРгет апсі Маігепа сіесісіесі іЬаІ іі ті§Ьі Ье Ьеііег Рог ІЬеіг зоп іо сопііпие Ьіз 
есіисаііоп Ьу Ро11оѵѵіп§ іп іЬе Рооізіерз оР Ьіз тоіЬег - а сЬетізІгу ІеасЬег іп Ні§Ь ЗсЬооІ 
апО Отзк’з МеОісаІ Іпзіііиіе (Ьіз РаіЬег ѵѵаз а таіЬешаІісіап апсі ІапО-зигѵеуог). ТЬегеРоге, 
іп 1952 Оеппасіу шоѵеё Іо Мозсо\ѵ ѵѵЬеге Ье епіегеё ІЬе Мозсоѵѵ Зіаіе Ыпіѵегзііу (М81І), 
апсі Ье §гас1иа1ес1 ѵѵііЬ а сЬешізІгу сіе§гее іп ОесетЬег 1957. Ніз ЬасЬеІог’з <3е§гее сіеаіі \ѵііЬ 
ІЬе ргоЫет оР зерагаііп§ Ьіб апсі Іл7 ізоіорез. АРЙег іЬіз Ье ]оіпесі іЬе Іпзіііиіе ѵѵаз зрііі 
іпіо Іѵѵо рагіз: N. N. Зетепоѵ Іпзіііиіе оР СЬетісаІ РЬузісз (ІСР) апсі N. М. Етапиеі 
Іпзіііиіе оР ВіосЬетісаІ РЬузісз. Аі іЬе ргезепі Ііте РгоР. О. Е. 2аікоѵ із ѵѵогкіп§ аі ІЬе N. 
М. Етапиеі Іпзіііиіе оР ВіосЬетісаІ РЬузісз (ІВР). 

Оеппасіу ѵѵаз огі ё іпа11у іпѵііесі Іо ІСР Ьу РгоРеззог Nіко1аі МагкоѵісЬ Етапиеі. 
ШОег Ьіз §иЫапсе, О.Е.2аікоѵ ёеРепсіеё іп 1963 Ьіз РЬ.Б. іЬезіз епііііеё “Сотрагізоп оР 
ІЬе Кіпеіісз апсі МесЬапізт оР Охісіаііоп оР іЬе Ог§апіс Сотроипёз іп Оазеоиз апсі Ьщшс! 
РЬазез іп 1963. Іп 1968, Ье ОеРепсІесІ а Носіог оР Зсіепсе іЬезіз епііііеё “ТЬе Коіе оРМесІіа 
іп КаОісаІ-СЬаіп Охіёаііоп Кеасііопз”. Іп 1970 Ье Ьесате а Риіі ргоРеззог. 

Іп 1966, Оеппасіу Ье§ап Іо Ьесоте іпѵоіѵесі ѵѵііЬ роіутег зсіепсе. N.М. Етапиеі 
сЬаг§е<3 2аікоѵ ѵѵііЬ іЬе огуапігаііоп оР ѵѵогк оп ргоЫетз аззосіаіеё ѵѵііЬ а§іп§ апсі 
зІаЬіІігаІіоп орроіутегз, апсі, Іаіег, ѵѵііЬ іЬе сотЬизІіоп оР роіутегіс таіегіаіз. Іп іЬе 1970з 
іЬеге ѵѵеге аЬоиІ 1000 зсіепіізіз іп іЬе 118. 8. В. ѵѵогкіп§ оп іЬезе ргоЫетз іпс1ис1іп§ 200 
зсіепіізіз Ргот ІСР ипсіег 2аікоѵ’з ІеаОегзЬір. ТЬе гезеагсЬ ѵѵаз сопсіисіеё оп аіі азресіз оР 
ІЬезе роіутег ргоЫетз іЬегтаІ ёе§гас1аІіоп, охісіаііоп, огопоіузіз, рЬоЮс!е§гас1аІіоп апсі 
гасііаііоп <Зе§гаёаІіоп, Ьусігоіузіз, Ьіо<іе§гас1аІіоп, тесЬапісаІ 0е§гас1а1іоп, ругоіузіз, апсі 
ЙаттаЬіІіІу. 

Айег “регезігоіка апсі ёе§гас!аІіоп” оР іЬе ІА8.8.К. іп 1991, іЬе пеѵѵ Киззіап 
§оѵегптепІ ОесгеазеО іЬе Ртапсіаі зиррогі оР зсіепсе зі§пійсап11у. 8о, О.Е.2аікоѵ Ьаз поѵѵ 
ѵѵііЬ Ьіт іп ІЬе N. М. Етапиеі Іпзіііиіе опіу 20 соѵѵогкегз (іпзіеасі оР200 іп 1970 -1980з). 
Рогіипаіеіу, Ье Ьаз §ооё зсіепіШс соорегаііоп ѵѵііЬ 28 гезеагсЬ сепіегз оР іЬе Рогтег 
118.8.К. (е.§., Киззіа, Шіаіпе, Оеог§іа, Веіогиззіа, Агтепіа, Тафкізіап) апсі тапу зоте 
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гезеагсЬ сепіегз оР Еигоре (Рогіи§а1, 8раіп, Ве1§іит, СгесЬ КериЫіс, Роіапсі, Ви1§агіа, 
Поташа). 

О.Е.2аікоѵ із ап оиізІапс1іп§ зсіепіізі ѵѵііЬ ехрегіізе іп \ѵіёе агеаз оР сЬетізІгу: 
сЬетісаІ апсі Ьіо1о§іса1 кіпеіісз, сЬетізІгу апсі рЬузісз оР роіушегз, Ьізіогу оР сЬетізІгу, 
ЪіосЬетізігу. 

Іп асШіііоп Іо Ьіз розіііоп аі іЬе N. М. Етапиеі Іпзіііиіе, Ье із а Іесіигег аі іЬе 
Мозсош Іпзіііиіе оРРіпе СЬетісаІ ТесЬпо1о§у. Не іаи§Ьі Ьіз зіисіепіз Ргогп Ьіз о\ѵп Ьоокз: 
Ое§гасіаііоп апсі ЗіаЪіІігаііоп оР Роіушегз, РЬузісаІ МеіЬоёз іп СЬетізІгу, апсі Асісі Каіпз 
апсі Епѵігоптепіаі РгоЫетз. О.Е.2аікоѵ Ьаз ѵѵгіііеп аЪоиі 1800 огі§іпа1 агіісіез, 130 
шопоёгарЬз (20 іп Киззіап апсі ПО іп Еп§1ізЬ), апсі 250 сЬаріегз іп 55 ѵоіишез. И із 
аррагепі Ргош іЬіз ѵѵогк іЬаі Ье Ьаз тасіе ѵаІиаЫе сопігіЪиііопз - а§іп§ апсі сіеѵеіортепі оР 
пе\ѵ зіаЪіІігегз Рог роіушегз, ог§апІ 2 аііоп оР іЬеіг іпсіизігіаі ргоёисііоп, ІіРе-ііте ргесіісііопз 
Рог изе апО зіога§е, апсі іЬе тесЬапізтз оР охісіаііоп, огопоіузіз, Ьуёгоіузіз, ЪіосІе§гасІаііоп, 
апсі 0есгеазіп§ оР роіутег ЙаттаЬіЫу. Ыеѵѵ теіЬоёз оР роіутег тосіісіісаііоп изіп§ іЬе 
ргосеззез оР сіе§гасіаііоп сѵеге іпігосіисесі іпіо ргасіісе Ьу 2аікоѵ. ТЬезе теіЬосіз аііоѵѵ іЬе 
ргосіисііоп оР пеѵѵ роіушегіс таіегіаіз \ѵііЬ ішргоѵесі ргорегііез. Мозі гесепііу, Ье із аізо 
ѵегу асііѵе іп іЬе Ріеіё оР зетісопсіисіогз апсі еіесігосопсіисііѵе роіушегз, роіутег Ыепёз, 
апсі роіутег сотрозііез іпсіиёіпё папосотрозііег. 

О.Е.2аікоѵ із а тетЬег оР тапу еёііогіаі Ьоагёз оР ритаіз риЫізЬеё іп Киззіа, 
Роіапсі, Ви1§агіа, іЬе Ы.8.А. апсі Еп§1апс1. Веіоѵ/ 1 \ѵои1с! Ііке іо §іѵе а Іізі оРЬіз асііѵііу іп 
іЬіз Ріеіё. 
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РЕФЕРАТЫ 


УДК 536.7 + 541.12 

Применение функций Лагранжа и пространственно-энергетических параметров 
для определения вероятности процессов. Липанов А.М., Кодолов В. И., Кораблев Г.А. 
ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. 2005. Том 7, № 1. С. 6-22. 

Предложен метод оценки вероятности протекания процессов с использованием 
функций Лагранжа. Обосновано введение пространственно - энергетического парамет- 
ра, полученноі о на основе анализа поверхностной и объемной энергий систем и вывода 
приведенной энергии. Приведены примеры применения предложенного метода для 
расчета энергий активации химических реакций. Табл. 5. Библиогр. 3. 

УДК 532.519.6 

Математическое моделирование течения жидкости с переменной структурой. Ли- 
панов А.М., Вахрушев А.А., Вахрушев А.В. ХИМИЧЕСКАЯ ФИЗИКА И 
МЕЗОСКОПИЯ. 2005. Т7. №1. С. 23-30. 

Сформулирована задача расчета течения жидкости, в которой при изменении 
макропараметров течения изменяются её микроструктурные свойства и реологическое 
поведение. Выполнены численные исследования течения жидкостей с температурным 
пределом изменения вязкости. Ил. 5. Библиогр. 8. 

УДК 541.123+621.791.052 

Расчет энергии активации объемной диффузии и самодиффузии элементов в твер- 
дых телах. Соловьев С.Д., Кораблев Г.А., Кодолов В. И. ХИМИЧЕСКАЯ ФИЗИКА И 
МЕЗОСКОПИЯ. 2005. Т7. №1. С. 31-40. 

В работе приведен метод расчета энергии активации самодиффузии и объемной 
диффузии атомов в твердых телах, основанный на эффективной энергии межатомного 
взаимодействия. Учитывается температурный фактор процесса диффузии атомов и его 
влияние на энергию межатомного взаимодействия. Показано, что результаты прове- 
денных многочисленных расчетов удовлетворительно согласуются с эксперименталь- 
ными данными по диффузии элементов в твердых телах. Табл. 3. Библиогр. 1 7. 
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УДК 541.1+536.4 


Влияние концентрации железа и давления в камере на параметры парогазовой 
смеси плазмогазодинамической установки. Бесогонов А.П., Кисаров Ю.Ф., Шиш- 
кина Л.В. ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. 2005. Том 7, № 1. С. 41-52. 

Предложена одномерная математическая модель расчета процесса неравновес- 
ной конденсации металла в сверхзвуковом сопле Лаваля. Установлено, что давление в 
камере плазмогазодинамической установки на количество конденсата влияет очень 
слабо, а увеличение содержания паров металла в исходной парогазовой смеси приводит 
к уменьшению концентрации конденсата в фиксированных сечениях сопла. Ил. 6. Биб- 
лиогр. 6. 

УДК 539.3 

Расчет потенциала парного заимодействия наночастиц. Вахрушев А. В., Липа- 
нов А.М. ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ.2005. Т7. №1. С. 53-62. 

Предложен метод расчета потенциала парного взаимодействия наночастиц на 
основе аппроксимации результатов численных расчетов метода молекулярной динами- 
ки. На основе потенциала парного взаимодействия наночастиц рассчитана зависимость 
предела прочности монодисперсного порошкового нанокомпозита от размеров состав- 
ляющих его наночастиц. Ил. 7 . Библиогр. 6. 

УДК 541.11 

Метастабильные самоассоциаты (малые кластеры) в системе медь+аргон. 

Моисеев Г.К. ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. 2005. Т.7. №1. С. 63-70. 

Рассчитаны термодинамические свойства и функции метастабильных самоассо- 
циатов Сиг, Си 3 , Си 6 . Методами термодинамического моделирования с использованием 
модели идеальных растворов продуктов взаимодействия рассчитаны составы конден- 
сированной и газовой фаз при Р = 10 5 Па и Т = 400-3000 К с учетом существования кла- 
стеров. Табл. 1. Ил. 2. Библиогр. 13. 

УДК 541.49+541.121 

Комплексообразование хрома(ІІІ) с моноаминными карбоксиметильными комп- 
лексонами и салициловой кислотой в водных растворах. Корнев В. И., Микрю- 
кова Г. А. ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. 2005. Том 7, № 1. С. 71-78. 

Исследованы взаимодействия в системах с участием моноаминных комплексонов 
(МИДА, ГЭИДА, НТА) и салициловой кислоты. Обнаружено формирование двух одно- 
роднолигандных комплексных соединений для салициловой кислоты различного про- 
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тонного состава, а также образование смешанолигандных комплексов с участием грех 
изученных комплексонов, в состав которых входит разное число протонов. В качестве 
основного использовался спектрофотометрический метод исследования, как вспомога- 
тельный применялся потенциометрический способ контроля. Показано, что процессы 
комплексообразования в тройных системах осложнены гидролизом и сильно зависят от 
кислотности среды. Выявлены механизмы процессов комплексообразования в тройных 
системах; установлены оптимальные области значений рН комплексообразования, со- 
став образующихся комплексов. Рассчитаны константы равновесия процессов комплек- 
сообразования, а также константы устойчивости образующихся комплексных частиц. 
Схематически представлено образование смешанолигандных комплексов хрома(ІІІ). 
Табл. 2. Ил. 3. Библиогр. 9. 

УДК 541.11. 

Стандартные энтальпии образования (СЭО) конденсированных метастабильных 
бинарных самоассоциатов (кластеров) различных элементов. Моисеев ГК. 
ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. 2005. Т.7. №1. С. 79-87. 

Для 25 бинарных конденсированных кластеров ІА, ПА и ІІІВ-ѴІВ подгрупп 
Системы исследованы зависимости известных величин ДН^ от размеров атомов, по- 
ложительных ионов и величин первых потенциалов ионизации для элементов каждой 
подгруппы. Сконструированы эмпирические уравнения для оценки величин ДН° 98 би- 
нарных кластеров. Табл. 2. Ил. 3. Библиогр. 12. 

УДК 593.139 

Исследование внутренней динамики трехатомных кластеров благородных газов 
методом выделения эффективных мод движения. Рыбаков А.А. ХИМИЧЕСКАЯ 
ФИЗИКА И МЕЗОСКОПИЯ. 2005. Т.7. №1. С. 88-98. 

В работе предлагается метод определения мод движения и эффективной размер- 
ности фазового пространства в задаче анализа внутренней динамики нелинейной мно- 
і очастичной системы. Метод протестирован на модели гармонического осциллятора. 
Представлены результаты исследования внутренней динамики трехатомных кластеров 
благородных газов с помощью данного метода. Ил. 6. Библиогр. 9. 
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УДК 541.11 

О метастабильных малых кластерах в жидком железе. Моисеев Г. К. 
ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. 2005. Т.7. №1. С. 99-103. 

Рассчитаны термодинамические характеристики метастабильного конденсиро- 
ванного кластера Рег и оценены ЛН^ для Рез и Ре 4 . С использованием методологии 
термодинамического моделирования определен состав жидкого железа при Р = 1 0 ' Па и 
Т = 1900-3100 К с учетом атомов и кластера Рег; оценено содержание в железе класте- 
ров Рез и Ре 4 . Табл. 2. Ил. 1. Библиогр. 11. 

УДК 622.273.2 

Бесцементная композиция для закладки выработанных пространств на основе 
техногенных материалов. Плеханова Т. А. ХИМИЧЕСКАЯ ФИЗИКА И 
МЕЗОСКОПИЯ. 2005. Т. 7. №1. С. 104-112. 

Для закладки выработанных пространств комбинатом «Норильский Никель» ис- 
пользуются составы на основе портландцемента с добавлением местного техногенного 
отхода - ангидрита. Эти закладочные смеси отличаются высокой стоимостью и низкой 
долговечностью. В этой работе рассматривается замена портландцемента на каустиче- 
ский магнезит. Составы на основе каустического магнезита имеют более высокие 
прочностные характеристики и являются более устойчивыми к сульфатной коррозии. 
Закладочные смеси на основе магнезиального вяжущего позволяют максимально ис- 
пользовать техногенные отходы предприятия такие как ангидрит, металлургические 
шлаки и карфосидерит, что снижает не только стоимость закладочных смесей, но и 
снижает экологическую нагрузку на регион. Табл. 3. Ил. 6. Библиогр. 6. 

РгоЫеіт оГ 1Ье сгеаііоп оГ папогеасіог» Гог іЬе 8уп*Ье8І8 оГ іпеЫІіспапорагІісІе» іп 
сагЬоп «ЬеІБ. Ьірапоѵ А.М., Кобоіоѵ Ѵ.І., КЬокЬгіакоѵ Ы.Ѵ., Оісіік А.А., Кобоіоѵа Ѵ.Ѵ., 
Зетакіпа МѴ. ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. 2005. Т7. №1. С. 113-120. 

Для получения металлических наночастиц, включая нанопроволоки, в углерод- 
ных оболочках целесообразно оценить возможности протекания сопряженных окисли- 
тельно-восстановительных реакций с участием ионов металла и органических соедине- 
ний с последующим экспериментальным подтверждением результатов моделирования. 
В качестве источников углерода предложены такие углеводородные соединения, как 
ароматические углеводороды, например, нафталин, антрацен или фенантрен, и поли- 
мерные функциональные соединения, например, поливиниловый спирт. Матрицами, 
которые имеют нанореакторы (полости или межслойные пространства, где протекает 
образование наночастиц), выбраны эвтектические расплавы хлоридов 36 металлов (Сг, 
Мп, Ре, Со, №, Си), гели поливинилового спирта или поливанадиевой кислоты и ее ин- 
теркалированных металлами производных. Табл. 2. Ил. 7. Библиогр. 6. 
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